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ABSTRACT 


The  hydrodesulfurization  of  benzothiophene  over  a 
Ni -Mo-y-Al 20 3  catalyst  was  studied  in  a  fixed  bed  microreac¬ 
tor  at  350  to  450°C  and  atmospheric  pressure.  Before  use, 
the  catalyst  was  first  reduced  in  situ  at  250°C  for  4  hours 
in  a  stream  of  hydrogen  and  then  presulfided  at  250°C  in  a 
10%  mixture  of  hydrogen  sulfide  in  hydrogen  for  another  4 
hours.  The  catalyst  was  kept  on  stream  for  up  to  40  hours 
for  a  given  set  of  reaction  conditions,  after  which  point  it 
was  discarded. 

Ethylbenzene  was  the  major  product;  dihydrobenzothio- 
phene  and  styrene  were  also  found  in  the  product  stream  in 
much  smaller  quantities  and  were  identified  as  intermedi¬ 
ates.  The  overall  reaction  rate  increased  from  350  to  400°C, 
but  was  virtually  identical  at  400  and  450°.  It  was  hypothe¬ 
sized  that  as  temperature  increases  beyond  400°,  the  major 
reaction  path  shifts  from  initial  hydrogenation  of  benzo¬ 
thiophene  to  dihydrobenzothiophene  to  initial  desulfuriza¬ 
tion  to  styrene,  with  the  decrease  in  the  rate  of  benzothio¬ 
phene  hydrogenation  caused  either  by  adsorption  limitations 
on  one  of  the  reactants  or  by  a  shift  in  the  equilibrium 
between  benzothiophene  and  dihydrobenzothiophene. 

Catalyst  activity  was  also  found  to  decrease  linearly 
with  time  on  stream,  at  a  rate  which  ranged  from  0.0  to  0.2% 
per  hour.  This  rate  of  deactivation  was  found  to  be  indepen¬ 
dent  of  time  on  stream  and  to  depend  only  on  the  reaction 
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conditions,  decreasing  with  an  increase  in  conversion 
benzothiophene  and  increasing  with  temperature.  It  is 
gested  that  the  loss  of  activity  was  due  to  a  loss  of 
from  the  catalyst,  and  that  this  was  inhibited  by  the 
hydrogen  sulfide  partial  pressures  in  the  reactor  at 
fractional  conversions. 
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1.  INTRODUCTION 


It  has  long  been  recognized  that  the  hydrodesulfuriza¬ 
tion  of  petroleum  fractions  has  a  number  of  benefits  beyond 
the  reduction  of  sulfur  oxide  emissions  from  combustion. 
Poisoning  of  precious-metal  reforming  catalysts  is  less¬ 
ened1,  the  ratio  of  middle  distillate  to  heavy  fuel  oil 
obtained  from  the  crude  is  increased2,  and  the  odor  and 
color  of  the  final  product  are  improved. 

Until  about  fifteen  years  ago,  application  of  hydrode¬ 
sulfurization  to  heavy  feeds  had  been  deterred  by  the  unique 
characteristics  of  these  feeds  and  the  associated  costly 
processing  requirements.  With  the  increased  development  in 
recent  years  of  alternative  fuel  sources,  in  particular  the 
high-sulfur  and  high-nitrogen  liquids  derived  from  tar 
sands,  coal,  and  oil  shale,  the  process  of  hydrodesulfuriza¬ 
tion  has  taken  on  a  new  importance. 

The  design  of  large  scale  catalytic  reactors  usually  is 
improved  by  an  understanding  of  the  reaction  mechanisms 
involved  as  well  as  by  the  availability  of  kinetic  data. 
Mechanistic  data  are  scarce  in  the  case  of  hydrodesulfuriza¬ 
tion  reactions,  and  most  of  the  available  data  have 
described  the  behavior  of  one  of  the  simpler  sulfur  com¬ 
pounds,  thiophene3"11.  Detailed  kinetic  data  are  even  more 
difficult  to  obtain,  and  the  available  literature  has  again 
focused  on  thiophene. 
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A  significant  amount  of  the  sulfur  present  in  crude 
fossil  fuels  is  present  as  benzo ( b) thiophenes ,  henceforth 
referred  to  as  benzothiophenes ,  however.  Martin  and  Grant12 
found  that  two-ring  thiophenes  accounted  for  22%  of  the  sul¬ 
fur  compounds  and  33%  of  the  thiophenic  compounds  in  a  sam¬ 
ple  of  Middle  East  crude,  and  14-25%  of  the  sulfur  compounds 
and  23-35%  of  the  thiophenic  compounds  in  a  number  of  sam¬ 
ples  of  American  crudes.  Clugston  et  al.13  also  found  sig¬ 
nificant  quantities  of  benzothiophenes  in  samples  of  bitumen 
from  the  Athabasca  and  Cold  Lake  districts,  and  in  a  sample 
of  heavy  oil  from  the  Lloydminster  area. 

Only  a  few  desulfurization  studies1 4-1 7  have  been  per¬ 
formed  using  benzothiophene  as  a  feed,  and  these  have  tended 
to  be  more  qualitative  than  quantitative.  In  addition,  there 
is  little  agreement  among  researchers  concerning  the  rela¬ 
tive  rates  of  various  reaction  pathways,  or  the  detailed 
mechanisms  of  even  the  major  paths.  For  this  reason,  it  was 
decided  to  investigate  the  benzothiophene  hydrodesulfuriza¬ 
tion  reaction  more  closely. 

Until  recently,  the  standard  catalyst  for  hydrodesul¬ 
furization  has  been  sulfided  cobalt-molybdenum  supported  on 
y-alumina.  Nickel-molybdenum  catalysts,  however,  are  coming 
under  much  closer  scrutiny,  partly  for  economic  reasons,  and 
partly  because  of  recent  data  showing  them  to  be  more  effec¬ 
tive  for  treating  feedstocks  high  in  nitrogen  and  oxygen 
compounds  than  had  previously  been  supposed. 
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The  mode  of  action  of  the  nickel  and  the  cobalt  cata¬ 
lysts  is  supposed  to  be  very  similar18  and  the  hydrodesul¬ 
furization  activities  also  appear  comparable19"20.  Greater 
activities  for  the  nickel  catalyst  have  been  reported  for 
hydrodenitrogenation2 0 ' 2 1  and  for  hydrogenation20,  and 
although  these  results  have  been  challenged19,  this  has  been 
another  major  reason  for  the  reevaluation  of  nickel-promoted 
catalysts . 

Because  no  detailed  studies  of  the  hydrodesulfurization 
of  benzothiophene  had  been  performed  using  a  nickel  cata¬ 
lyst,  and  because  this  catalyst  may  become  as  industrially 
important  in  future  as  the  more  heavily-studied  cobalt  cata¬ 
lysts,  it  was  decided  to  use  a  nickel-molybdenum  catalyst 
supported  on  alumina  in  the  present  study. 


1 


■ 


2.  LITERATURE  REVIEW 


2.1  Hydrodesulfurization  of  Benzothiophene 

Early  work  on  the  hydrodesulfurization  of  benzothio¬ 
phene  was  more  concerned  with  predicting  the  main  reaction 
products  or  with  proving  the  effectiveness  of  a  given  cata¬ 
lyst  than  with  attempting  to  discover  a  reaction  network. 

In  I960,  Thompson  and  his  coworkers22,  in  attempting  to 
devise  a  method  of  using  the  products  of  hydrodesulfuriza¬ 
tion  as  a  means  for  identifying  sulfur  compounds,  found  that 
only  ethylcyclohexane  was  produced  when  the  reaction  took 
place  over  an  alumina-supported  palladium  catalyst  at  155°C. 
That  same  year,  Hoog  et  al.23,  working  with  a  cobalt-molyb¬ 
denum  catalyst,  and  Yamada24,  working  with  a  molybdenum  tri¬ 
sulfide  catalyst,  both  found  only  ethylbenzene  as  a  product. 
Later,  Hopkins  et  al.25,  desulfurizing  benzothiophene  in 
boiling  methanol  with  Raney  nickel,  again  found  only  ethyl¬ 
benzene  produced. 

The  first  attempt  to  study  the  mechanism  of  benzothio¬ 
phene  hydrodesulfurization  in  detail  was  that  of  Givens  and 
Venuto16  in  1970.  Using  a  series  of  alumina-supported  cata¬ 
lysts  (including  two  Co-Mo  versions)  in  a  cont i nuous- f low 
reactor  at  atmospheric  pressure  and  temperatures  ranging 
from  300  to  400°C,  these  authors  found  both  ethylbenzene  and 
dihydrobenzothiophene  as  products.  Because  the  latter  only 
appeared  at  higher  space  velocities  and  because  it  appeared 
in  increasing  amounts  at  lower  temperatures  with  no  new 
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products,  it  was  hypothesized  that  the  reaction  was  a  two- 
step  conversion  of  benzothiophene  through  the  intermediate 
dihydrobenzothiophene . 

The  mechanism  postulated  was  that  of  a  rapidly  occur¬ 
ring  equilibrium  between  benzothiophene  and  dihydrobenzo¬ 
thiophene,  with  a  slow,  irreversible  desulfurization  step 
draining  this  equilibrium.  This  was  supported  by  data  from 
tests  of  the  desulfurization  rates  of  mixtures  of  benzothio¬ 
phene  and  dihydrobenzothiophene  in  varying  proportions, 
which  indicated  that  an  equilibrium  between  the  two  was 
reached  quickly,  and  for  this  reason  the  two  appeared  to 
form  ethylbenzene  at  comparable  rates.  It  was  suggested  that 
if  dihydrobenzothiophene  were  the  source  of  ethylbenzene, 
desulfurization  might  be  a  multistep  process  involving  first 
a  cleavage  of  one  carbon-sulfur  bond  to  give  a  discrete 
mercaptan  intermediate,  followed  by  a  second  C-S  bond  scis¬ 
sion  to  give  the  hydrocarbon  plus  hydrogen  sulfide.  If 
benzothiophene  were  the  source  of  ethylbenzene,  sulfur  would 
be  lost  in  a  one-step  process  involving  the  synchronous 
breaking  of  two  C-S  bonds,  producing  the  hydrocarbon  without 
a  mercaptan  intermediate. 

These  results  were  later  corroborated  by  De  Beer  et 
al.26,  who  used  benzothiophene  hydrodesulfurization  to  test 
the  relative  activities  of  a  series  of  alumina-supported 
molybdenum  and  tungsten  catalysts  promoted  by  nickel  and 
cobalt.  Using  a  fixed  bed  reactor  operating  at  atmospheric 
pressure,  these  authors  also  found  ethylbenzene  and 
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dihydrobenzothiophene  as  the  only  products,  and  agreed  with 
the  hypothesis  that  dihydrobenzothiophene  was  an  intermedi¬ 
ate  in  the  reaction. 

The  quantitative  results  of  this  earlier  investigation 
may  be  questioned,  however,  since  the  products  of  the  reac¬ 
tion  were  condensed  prior  to  analysis  (presumably  at  room 
temperature),  rather  than  using  an  on-line  gas  sampling 
technique.  Consequently,  it  is  possible  that  some  of  the 
more  volatile  compounds  which  could  have  been  produced  in 
this  reaction  may  have  been  eluted  in  the  gas  stream.  This 
would  affect  the  relative  amounts  of  the  products  analyzed, 
and  could  have  resulted  in  a  failure  to  detect  other  inter¬ 
mediates  or  side  reaction  products. 

Also,  the  catalyst  was  used  in  a  highly  reduced  state, 
having  been  reduced  in  H2  at  over  500°C  prior  to  use.  High 
reduction  temperatures  have  been  reported  to  seriously 
degrade  the  activity  of  hydrodesulfurization  catalysts27, 
and  it  thus  seems  possible  that  many  of  the  catalyst  sites 
may  have  been  destroyed  by  the  time  that  it  was  used.  If  so, 
little  (if  any)  sulfiding  of  the  catalyst  by  the  H2S  pro¬ 
duced  in  the  reaction  may  have  occurred.  (The  effects  of 
various  methods  of  catalyst  pretreatment  are  reviewed  in 
more  detail  in  Section  2.5.)  This  is  supported  by  the  fact 
that  the  authors  found  no  effect  of  presulfiding  on  catalyst 
activity,  and  as  a  result,  their  catalyst  may  have  borne 
little  resemblance  to  most  current  industrial  hydrotreating 
catalysts . 
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The  next  study  of  this  reaction  was  performed  in  1974, 
by  Bartsch  and  Tanielan15,  who  also  used  a  Co-Mo-y-Al 203 
catalyst  in  a  continuous  flow  reactor  at  atmospheric  pres¬ 
sure,  with  temperatures  in  the  same  range  as  had  been  used 
by  Givens  and  Venuto.  In  this  study,  ethylbenzene  was  again 
the  only  product  found.  Bartsch  and  Tanielan  used  their  data 
to  derive  a  pseudo-first  order  rate  expression  and  an  appar¬ 
ent  activation  energy  (4.9  kcal/mol)  for  the  overall  reac¬ 
tion. 

This  investigation  unfortunately  included  most  of  the 
deficiencies  of  the  previous  one.  The  reaction  products  were 
again  condensed  at  room  temperature  prior  to  analysis,  and 
while  the  prereduction  and  presulfiding  conditions  were  less 
severe  (400°C),  it  is  still  possible  that  the  catalyst  may 
have  been  adversely  affected. 

Styrene  was  discovered  in  the  reaction  products  along 
with  ethylbenzene  by  Furimsky  and  Amberg14  in  1976.  These 
authors  used  a  series  of  unsupported  molybdenum  disulfide 
catalysts  as  well  as  a  commercial  Co-Mo-y-Al 203  catalyst  and 
a  chromia  catalyst  in  a  pulse  microreactor  at  atmospheric 
pressure  and  temperatures  ranging  between  250-400°C.  Both 
benzothiophene  and  dihydrobenzothiophene  were  employed  as 
feeds  in  this  study. 

Apparent  activation  energies  based  on  Arrhenius  plots 
were  calculated  for  the  desulfurization  of  benzothiophene 
and  dihydrobenzothiophene  (12  kcal/mol  and  7  kcal/mol, 
respectively).  In  addition,  the  temperature  dependence  for 
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the  yield  of  benzothiophene  from  dihydrobenzothiophene  was 
used  to  estimate  the  activation  energy  of  the  dehydrogena¬ 
tion  of  dihydrobenzothiophene  (14.8  kcal/mol).  The  fact  that 
a  higher  activation  energy  was  found  for  dehydrogenation 
than  for  desulfurization  implied  that  with  increasing  tem¬ 
perature,  the  rate  of  dehydrogenation  would  become  increas¬ 
ingly  important. 

Based  on  these  findings,  Furimsky  and  Amberg  proposed 
several  reaction  routes  for  the  overall  desulfurization,  as 
shown  in  Figure  1 . 

Reaction  route  a  suggests  participation  of  the  /?-hydro- 
gen  of  the  thiophene  ring  in  carbon-sulfur  bond  cleavage,  a 
hypothesis  supported  by  the  earlier  finding  of  Givens  and 
Venuto  that  a  marked  decrease  in  the  desulfurization  rate  of 
benzothiophene  was  found  when  this  hydrogen  was  substituted 
with  a  methyl  group.  Route  b  involves  C-S  bond  cleavage  with 
the  participation  of  hydrogen  as  the  initial  step,  though 
2-mercaptophenylacetylene  is  too  unstable  and  2-mercaptosty- 
rene  too  quickly  desulfurized  for  the  authors  to  have 
expected  to  have  detected  them  in  their  products.  Route  c 
represents  the  saturation  of  the  thiophene  ring  followed  by 
either  migration  of  one  of  the  /2-hydrogens  with  subsequent 
C-S  bond  cleavage  or  by  double  C-S  bond  scission  and  satura¬ 
tion  of  the  aliphatic  carbon. 

The  presence  of  styrene  in  the  product  stream  was 
verified  by  Bartsch  and  Tanielan28  in  1977,  while  studying 
the  hydrodesulfurization  of  benzothiophene  over  prereduced 


'  1 


9 


Figure  1  -  Furimsky  and  Rmberg  Reaction  Scheme 
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Fe203  catalysts  at  atmospheric  pressure  and  300-400°C.  These 
authors  supported  the  view  that  C-S  bond  cleavage  was  the 
primary  initial  step  in  the  overall  desulfurization  reac¬ 
tion,  and  hypothesized  two  possible  main  reaction  schemes 
similar  to  routes  b  and  c  of  Furimsky  and  Amberg.  Bartsch 
and  Tanielan  also  allowed  cleavage  of  the  S-C£  bond  as  an 
alternative  to  the  S-C,  bond,  with  2-pheny lethanethi ol  and 
2-pheny lethylenethi ol  the  intermediates  produced.  A  pseudo- 
first  order  rate  equation  was  derived  for  the  overall  reac¬ 
tion,  with  apparent  activation  energies  of  between  12.5  and 
14.0  kcal/mol  found  for  the  catalyst  reduced  at  400°C. 

No  styrene  was  found  when  the  catalysts  were  presul¬ 
fided,  however,  and  Bartsch  and  Tanielan  suggested  that  the 
sulfurized  sites  aided  in  hydrogenation.  This  was  supported 
by  data  showing  a  feed  consisting  of  pure  styrene  to  be 
completely  hydrogenated  when  subjected  to  their  reaction 
conditions . 

It  was  also  found  in  this  study  that  the  amount  of  sty¬ 
rene  produced  using  prereduced  catalysts  decreased  with  time 
on  stream,  and  these  authors  concluded  that  these  sulfurized 
hydrogenation  sites  were  formed  during  the  first  few  hours 
on  stream. 

The  early  results  of  Thompson  found  some  support  from 
Rollman29  in  1977,  who  also  found  traces  of  ethylcyclohexane 
while  hydrotreating  mixtures  of  sulfur  and  nitrogen  com¬ 
pounds  in  the  liquid  phase  at  2.1-10.3  MPa  and  300-450°C 
over  a  Co-Mo  catalyst.  Ethylbenzene  was  the  only  other 
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product  found.  No  mechanism  was  postulated  in  this  study; 
the  main  concern  of  this  investigation  was  the  relative 
reactivities  of  the  various  compounds  tested. 

The  next  researcher  to  investigate  benzothiophene 
hydrodesulfurization  was  Daly17  in  1978.  This  author  used  a 
liquid  phase  st i r red-batch  minireactor  charged  with  a  Co-Mo- 
y-Al  2 0 3  catalyst,  operating  at  200-400°C  and  8.6  MPa.  An 
apparent  activation  energy  of  16.1  kcal/mole  was  calculated 
for  the  overall  reaction. 

Traces  of  1-  and  2-phenylethanethiol  were  detected  in 
the  products  along  with  ethylbenzene  and  styrene.  It  was 
suggested  that  these  thiols  resulted  from  a  back-reaction 
between  hydrogen  sulfide  and  an  unsaturated  hydrocarbon, 
probably  styrene.  This  was  supported  by  the  result  that  an 
increased  hydrogen  sulfide  concentration  produced  an 
increase  in  the  concentrations  of  these  back-reaction  prod¬ 
ucts,  and  that  in  the  presence  of  elemental  sulfur  (known  to 
catalyze  the  reaction  of  hydrogen  sulfide  with  olefins30), 
there  was  also  an  increase  in  these  products. 

Daly  proposed  a  reaction  scheme  similar  to  that  of 
Furimsky  and  Amberg,  but  incorporated  steps  in  which  styrene 
reacts  to  form  1-  and  2-phenylethanethiol,  and  dihydrobenzo- 
thiophene  reacts  to  form  2-phenylethanethiol,  which  could 
subsequently  be  desulfurized  to  form  ethylbenzene. 

The  pretreatment  in  this  study  was  also  at  a  high  tem¬ 
perature,  and  in  addition,  mass  balances  were  not  computed 
for  the  products,  which  may  limit  the  validity  of  the 
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quantitative  results. 

In  1978,  Kilanowski  et  al.31  used  a  Co-Mo-y-Al 203  cata¬ 
lyst  in  a  pulse  microreactor  to  test  the  relative  reactivi¬ 
ties  of  various  thiophenic  compounds  at  atmospheric  pressure 
and  temperatures  of  300-400°C.  Styrene  and  ethylbenzene  were 
the  only  products  found  from  benzothiophene  hydrodesulfuri¬ 
zation,  and  it  was  concluded  that  the  reaction  network  could 
be  approximated  by  a  simple  sulfur  extrusion.  The  absence  of 
sulfides  and  mercaptans  in  the  products  was  taken  as  evi¬ 
dence  of  the  carbon-sulfur  bond  scission  being  a  slow  reac¬ 
tion  step. 

In  this  study,  the  pretreatment  consisted  of  heating 
the  catalyst  to  400°C  in  hydrogen,  then  presulfiding  at  this 
temperature  with  a  10%  mixture  of  K2S  in  H2.  This  presulfid¬ 
ing  may  have  had  only  a  short-lived  effect  on  the  catalyst, 
though,  since  pure  hydrogen  was  passed  over  the  catalyst 
constantly,  with  0.5  fil  pulses  injected  only  every  35  min¬ 
utes  to  1  hour.  Although  little  data  are  available  on  cata¬ 
lysts  which  have  been  presulfided  and  then  subjected  to 
severe  reducing  conditions,  it  seems  likely  that  a  signifi¬ 
cant  amount  of  sulfur  would  be  lost  from  the  catalyst  in 
such  a  case. 

These  authors  also  speculated  on  the  nature  of  the 
adsorption  of  benzothiophenes  on  the  catalyst  surface.  It 
was  suggested  that  single  anion  vacancies  could  accommodate 
the  sulfur  atom  in  a  "single-point"  adsorption,  as  had  been 
proposed  for  thiophene,  whereas  an  adjacent  pair  of 
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vacancies  might  allow  for  a  flat  adsorption,  involving  the  n 
electrons  of  the  thiophenic  ring. 

A  second  study  of  benzothiophene  hydrodesulfurization 
by  Kilanowski  and  Gates32  was  published  the  next  year.  The 
same  catalyst  was  used  as  previously,  although  this  time  the 
reaction  took  place  in  a  flow  microreactor  at  atmospheric 
pressure  and  temperatures  of  252,  302,  and  332°C.  In  this 
study,  ethylbenzene  was  the  only  reaction  product  found  in 
greater  than  trace  amounts. 

Only  one  of  a  number  of  empirical  rate  equations  tried 
was  found  to  even  roughly  fit  the  data  at  all  three  tempera¬ 
tures,  and  this  equation  fit  the  data  much  more  poorly  at 
332°C  than  at  the  two  lower  temperatures. 

It  was  suggested  that  the  difficulty  of  finding  a  rate 
equation  applicable  to  the  full  range  of  temperatures  stud¬ 
ied  was  due  to  a  complicated  mechanism  and  changes  in  the 
kinetics  with  temperature.  From  the  data,  an  apparent  acti¬ 
vation  energy  of  20  ±3  kcal/mol  was  calculated;  the  respec¬ 
tive  heats  of  adsorption  of  benzothiophene  and  H2S  were  cal¬ 
culated  as  -15  ±10  and  -6  ±6  kcal/mol,  with  the  large 
uncertainties  in  these  latter  figures  mainly  due  to  scatter 
in  the  data  at  252°C. 

In  1979,  Nag  et  al.33  desulfurized  benzothiophene  in  a 
high  pressure  (7.2  MPa)  batch  reactor  charged  with  a  commer¬ 
cial  Co-Mo-y-Al 203  catalyst  at  300°C.  Pretreatment  consisted 
of  presulfiding  the  catalyst  for  2  hours  at  400°C  in  a 
mixture  of  10%  H2S  in  H2.  Ethylbenzene  and  styrene  were  the 
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only  products  found,  in  the  approximate  ratio  of  5:1.  A 
pseudo-first  order  rate  constant  for  the  disappearance  of 
benzothiophene  of  2920  cm3/g  catalyst-hr  was  calculated. 

Recently,  Rwart,  Schuit,  and  Gates34  proposed  a  model 
for  the  adsorption  of  benzothiophene  onto  the  surface  of 
sulfided  catalysts,  which  they  claimed  was  much  superior  to 
both  the  "one-point"  and  the  "flat"  adsorption  models  which 
had  generally  been  used  to  describe  the  state  of  thiophenic 
molecules  prior  to  desulfurization.  This  new  model  was  then 
used  to  make  specific  predictions  about  the  paths  that  the 
desulfurization  reaction  should  take. 

These  authors  claimed  that  the  expected  bonding  charac¬ 
teristics  of  thiophenic  molecules  would  lead  to  adsorption 
at  the  C,-C2  bond,  rather  than  at  the  sulfur  atom:  the 
unshared  nonbonded  electron  pairs  on  sulfur  should  be  so 
involved  with  the  resonance  of  the  thiophene  nucleus  that 
the  sulfur  atom  should  have  little  residual  electron 
availability  for  direct  coordination  by  catalyst.  It  was 
suggested  instead  that  thiophenic  molecules  may  be  chemi¬ 
sorbed  such  that  the  C,-C2  bond  is  coordinated  at  an  anion 
vacancy  on  the  catalyst,  with  the  adjacent  sulfur  center 
then  interacting  with  a  neighboring  sulfur  atom  on  the  cata¬ 
lyst  surface.  When  the  thiophene  ring  becomes  bonded  via  the 
C(-C2  position  at  the  anion  vacancy,  the  result  is  expected 
to  be  a  change  in  the  electron  distribution  of  the  ring, 
making  the  sulfur  atom  electron-deficient,  and  consequently 
promoting  its  bonding  with  another  sulfur  atom  of  high 
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electron  density  on  the  catalyst  surface.  These  changes 
should  facilitate  hydrogen  addition  to  the  coordinated  0,-02 
centers  and  formation  of  a  dihydrogenated  intermediate.  This 
adsorbed  species  is  depicted  in  Figure  2. 

A  benzothiophene  molecule  adsorbed  at  the  C,-C2  bond 
could  provide  several  alternative  steps  to  simple  hydrogena¬ 
tion  and  desorption  to  form  dihydrobenzothiophene .  If  par¬ 
tial  hydrogenation  occurs  at  the  C,  atom,  electron  delocali¬ 
zation  could  lead  to  rupture  of  the  S-C,  bond,  which  would 
produce  either  2-mercaptostyrene  if  the  sulfur  atom  were 
hydrogenated,  or  styrene,  if  the  sulfur  atom  became  fully 
bonded  to  the  catalyst  surface  with  subsequent  C-S  bond 
scission.  Alternatively,  the  partially  bonded  protons  in  the 
structure  of  Figure  2  could  become  involved  in  covalent 
bonds  to  the  C,  and  C2  atoms  by  a  concerted  process,  with 
subsequent  ^-elimination  leading  to  the  same  products. 

These  authors  agreed  with  the  earlier  suggestion  of 
Kilanowski31  that  a  one-point  site  might  better  account  for 
the  high  reactivity  of  hydrogenated  compounds  lacking  the  n 
electrons  necessary  for  this  "three-point"  model.  They  con¬ 
sidered  their  model  superior  to  the  flat  adsorption  model  in 
that  it  did  not  require  the  interaction  of  two  separate  <r 
complexes . 

In  summary,  the  mechanism  by  which  benzothiophene  is 
desulfurized  is  still  unclear.  This  reaction  has  been  stud¬ 
ied  both  in  the  gas  phase  at  atmospheric  pressure  and  in  the 
liquid  phase  at  high  pressures.  Only  Co-Mo-y-Al 203  and 


' 


16 


Figure  2  -  Kwart,  Schuit,  and  Gates  Rdsorption  Model. 
The  dashed  lines  indicate  a  weaker  interaction  between 
centers  than  solid  lines;  the  wavy  line  indicates  a 
longer  or  more  remote  bond  between  centers. 
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Fe203-y-Al 203  have  been  employed  as  catalysts,  and  both  have 
been  pretreated  by  various  methods.  It  is  generally  agreed 
that  ethylbenzene  is  the  major  product,  but  both  styrene  and 
dihydrobenzothiophene  have  been  found  in  small  quantities  in 
the  product  stream,  and  it  appears  that  either  or  both  may 
be  intermediates  in  the  overall  desulfurization  reaction. 

The  kinetics  appear  to  be  complicated,  and  may  vary  signifi¬ 
cantly  with  temperature,  type  of  catalyst  used,  method  of 
pretreatment,  and  length  of  time  that  the  catalyst  has  been 
on  stream. 


2.2  Chemistry  of  Sulfided  Hydrodesulfurization  Catalysts 

Attempts  to  identify  the  catalyt ically  active  species 
present  in  hydrodesulfurization  catalysts  started  with 
Richardson35  in  1964.  By  measuring  the  magnetic  susceptibil¬ 
ity  of  various  cobalt  species  in  a  series  of  Co-Mo-y-Al 203 
catalysts  and  correlating  this  with  desulfurization  activ¬ 
ity,  Richardson  characterized  the  catalyst  as  a  mixture  of 
Al203,  CoA1204  (both  inactive),  Co9S8  (very  slightly 
active),  MoS2  (moderately  active  itself,  but  requiring  pro¬ 
motion  by  "active  cobalt"  for  maximum  activity),  and  possi¬ 
bly  MoO 2  (mildly  active).  This  "active  cobalt"  was  consi¬ 
dered  to  be  a  non-reducible  form  of  the  metal,  probably  in 
an  octahedral  coordination,  while  the  inactive  cobalt  was 
likely  confined  to  tetrahedral  sites. 
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Richardson  also  found  that  catalysts  with  initial  Co:Mo 
ratios  of  less  than  0.3  had  the  same  final  composition  when 
calcined  in  the  range  538-650  °C,  but  at  higher  initial 
Co:Mo  ratios,  the  final  composition  was  dependent  on  both 
cobalt  content  and  heat  treatment  method.  The  catalyst  found 
to  have  the  greatest  activity  had  a  Co:Mo  ratio  of  0.2,  and 
an  active  Co:Mo  ratio  of  0.18. 

Two  models  appeared  in  the  early  1970's  to  explain  the 
surface  chemistry  of  tungsten  and  molybdenum  catalysts  pro¬ 
moted  with  nickel  and  cobalt.  Both  of  these  models  described 
the  effect  of  the  promoters  as  an  increase  in  the  reduction 
of  the  catalyst  ions,  but  the  mechanisms  by  which  this  was 
achieved  were  significantly  different. 

The  first  of  these  two  models  was  the  intercalation 
model,  proposed  by  Voorhoeve  and  Stuiver36'38  and  Farragher 
and  Cossee39.  It  had  been  found  that  MoS2  and  WS2  have  a 
layer  structure  in  which  the  metal  is  surrounded  by  a 
trigonal  prismatic  coordination  of  six  sulfur  atoms.  The 
stacking  of  the  sulfur  layers  can  be  either  hexagonal  or 
rhombic;  in  practice  both  are  found,  with  frequent  defect 
structures.  Although  intercalation  is  impossible  for  perfect 
MoS 2  and  WS2  crystals,  it  was  postulated  that  nickel  and 
cobalt  can  intercalate  appreciably  if  the  crystals  are 
small . 

Based  on  kinetic  studies  of  the  hydrogenation  of  cyclo¬ 
hexene  and  benzene  over  a  sulfided  Ni-W  catalyst,  Voorhoeve 
and  Stuiver  concluded  that  the  active  sites  were  sulfur- 
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deficient  anion  vacancies  in  the  metal.  This  was  later  sup¬ 
ported  by  Wentrcek  and  Wise40,  who  found  a  positive  rela¬ 
tionship  between  hydrodesulfurization  activity  and  elec¬ 
tronic  "holes"  in  the  surface  of  a  single  MoS2  crystal,  and 
by  Aoshima  and  Wise41,  who  found  that  a  small  number  of  S2" 
vacancies  in  this  crystal  was  beneficial  to  hydrodesulfuri¬ 
zation  activity. 

From  the  relationships  between  an  ESR  resonance 
ascribed  to  W3+  and  H2S/H2  ratio  and  catalyst  nickel  con¬ 
tent,  and  from  the  fact  that  increasing  H2S  concentration 
decreased  the  activity  of  the  catalyst,  Voorhoeve  concluded 
that  W3  +  was  the  catalyt ically  active  site.  This  was  sup¬ 
ported  by  a  correlation  between  intensity  of  the  ESR  signal 
and  catalytic  activity  for  benzene  hydrogenation. 

The  discovery  that  the  alumina  carrier  had  no  effect  on 
the  types  of  sites  produced,  along  with  the  results  from 
optical  studies  and  X-ray  line  broadening  measurements  led 
Voorhoeve  and  Stuiver  to  conclude  that  nickel  was  located  in 
octahedral  holes  between  layers  of  WS2,  which  produced  an 
increase  in  electron  concentration  in  the  WS2  sheets  and 
thus  an  increase  in  the  number  of  active  tungsten  ions  on 
the  surface.  The  lack  of  correlation  between  the  magnitude 
of  the  W 3  +  signal  and  activity  in  cyclohexene  hydrogenation 
led  to  the  conclusion  that  this  compound  was  adsorbed  on  a 
different  type  of  site,  possibly  also  W3+  ions  but  in  a  dif¬ 
ferent  coordination. 


, 


20 


Farragher  and  Cossee  extended  this  model  with  further 
optical  and  X-ray  measurements,  showing  that  the  nickel  ions 
tended  to  be  located  on  the  edges  of  the  crystals. 

Maximum  Ni:W  ratios  based  on  this  extension  of  the 
edge- intercalat ion  model  were  predicted  to  be  in  the  range 
of  0.3  to  1;  it  was  suggested  that  varying  crystal  size 
could  explain  the  wide  variations  in  the  optimum  values  of 
this  parameter  reported  in  the  literature.  In  addition, 
nickel  and  cobalt  were  shown  to  be  both  physically  and  ther¬ 
modynamically  ideal  for  occupation  of  octahedral  holes  in 
the  crystal  structure. 

The  two  types  of  sites  hypothesized  by  Voorhoeve  and 
Stuiver  were  seen  as  a  single  W3+  ion  connected  to  an  anion 
vacancy  and  a  combination  of  two  W3+  ions  with  interconnect¬ 
ing  S2'  ions  absent. 

Although  this  model  was  considered  to  be  applicable  to 
the  Co-Mo  system  as  well,  ESR  data  for  the  corresponding 
compounds  had  not  confirmed  this  extension,  despite  attempts 
by  Lo  Jacono  et  al.42  and  others43.  Lo  Jacono  mentioned  that 
the  absence  of  an  ESR  signal  was  not  proof  of  the  absence  of 
Mo 3  +  ,  though,  due  to  the  possible  existence  of  these  ions  as 
diamagnetic  pairs. 

A  second  model  was  developed  by  Schuit  and  Gates43  and 
Hagenbach,  Courty  and  Delmon44  to  better  describe  sulfided 
Co-Mo  catalysts.  It  had  been  shown  that  although  nickel  pre¬ 
fers  an  octahedral  coordination,  cobalt  has  a  stronger 
preference  for  tetrahedral  sites45.  As  a  result,  these 
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authors  characterized  cobalt  as  incorporated  into  tetrahed¬ 
ral  sites  below  a  layer  of  ions  on  the  catalyst  surface. 

In  this  model,  y-Al203  is  assumed  to  consist  of  parti¬ 
cles  formed  by  a  one-dimensional  stacking  of  alternating  C 
and  D  layers,  where  the  D  layers  have  only  octahedral  Al , 
and  the  C  layers  have  an  equal  number  of  octahedral  and  tet¬ 
rahedral  sites.  Mo03  is  present  in  a  monolayer  on  the 
alumina  surface  and  is  arranged  with  the  opposite  symmetry 
to  the  top  alumina  layer.  To  compensate  for  the  excess 
positive  charge,  an  anion  "capping"  layer  is  assumed  to  be 
located  on  top  of  the  monolayer,  consisting  of  S2"  ions  in 
sulfided  catalysts  and  02"  ions  in  oxidic  and  partially 
reduced  catalysts.  Incorporation  of  Co2*  in  the  support  is 
assumed  to  result  in  the  location  of  these  cations  in  a  tet¬ 
rahedral  position  somewhere  below  the  monolayer.  This  is 
accompanied  by  expulsion  of  Al3  +  ions  from  the  support, 
which  then  must  find  tetrahedral  positions  in  the  monolayer, 
which  in  turn  forces  the  Mo  ions  to  find  places  in  the  octa¬ 
hedral  interstices.  As  a  result,  fewer  Mo  ions  are  incorpor¬ 
ated  and  the  ratio  of  capping  ions  to  molybdenum  ions  is 
increased.  This  was  supported  by  pulse-flow  data  of  De  Beer 
et  al.45,  who  showed  that  the  addition  of  Co  to  a  Mo-y-Al203 
catalyst  increased  the  uptake  of  H2S  and  by  the  diffuse 
reflectance  spectroscopy  of  Gajardo  et  al.46,  confirming  the 
presence  of  tetrahedral  Co2+  as  CoAl204  in  oxidic  Co-Mo- 
y-Al  2 0 3  .  Gajardo  et  al  also  found  a  significant  amount  of 
octahedral  Co3+,  however,  and  suggested  a  Co-Mo  bilayer 
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might  better  account  for  some  aspects  of  this  system. 

Additional  support  for  the  bilayer  character  of  the 
monolayer  model  was  presented  by  Delannay  et  al.47  whose  ion 
scattering  spectrometry  experiments  confirmed  that  molybde¬ 
num  was  essentially  confined  to  the  surface  layer  of  nonre- 
duced  catalysts,  and  that  cobalt  was  present  in  a  single 
layer  between  the  y-Al203  support  and  the  molybdenum  surface 
layer . 

In  the  monolayer  model,  the  stability  of  the  monolayer 
should  also  be  increased,  since  part  of  the  task  of  binding 
it  to  the  surface  layer  would  be  taken  over  by  the  alumina 
support.  Hagenbach  and  his  coworkers44  have  found  that  in 
small  quantities,  Co  does  in  fact  bring  about  a  better 
crystallinity  as  well  as  a  decrease  in  surface  area  in 
Co-Mo-y-Al 203 ,  although  at  very  low  concentrations  the 
incorporation  of  Co  seemed  harmful  rather  than  beneficial  to 
catalyst  activity. 

Another  consequence  of  the  monolayer  model  is  the  for¬ 
mation  in  the  sulfided  catalyst  of  some  Mo4+  and  Mo5+  ions, 
depending  whether  the  top  Al203  layer  is  a  C  or  a  D  layer: 
replacement  of  02~  ions  in  the  capping  layer  with  the  much 
larger  S2'  ions  results  in  extra  anion  vacancies  when  the 
latter  ions  are  removed  with  simultaneous  reduction  of  the 
underlying  Mo  ions.  This  aspect  of  the  model  received  some 
support  from  the  thiophene  hydrodesulfurization  data  and 
spectroscopic  measurements  of  Lipsch  and  Schuit48-50.  These 
authors  concluded  that  cobalt  was  almost  entirely  present  in 
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a  tetrahedral  coordination,  but  they  also  characterized  some 
of  the  active  sites  as  including  Mo4+,  Mo5  +  ,  and  Mo6+  ions 
in  an  octahedral  environment.  Support  for  the  view  that 
extra  anion  vacancies  are  produced  by  sulfiding  and  that 
these  are  catalyt ically  active  was  given  by  the  gravimetric 
data  of  Massoth51-52  in  1975  and  1977. 

Dufaux,  Che,  and  Naccache53  found  that  the  intensity  of 
the  Mo 5  +  signal  appeared  only  at  Mo03  concentrations  above 
10%,  adding  support  to  the  concept  of  an  inert  Mo03  mono- 
layer.  Lo  Jacono  et  al.42  also  found  this  effect,  but  also 
found  that  a  weak  Mo5  +  signal  disappeared  with  the  addition 
of  a  small  amount  of  Co.  They  concluded  that  the  ESR  inert¬ 
ness  of  the  first  amounts  of  Mo03  might  instead  be  ascribed 
to  a  solution  of  Mo  in  the  support.  In  other  words,  both  Co 
and  Mo  dissolve  in  the  alumina,  and  although  initially  the 
solubility  is  restricted,  once  the  two  ions  are  present 
simultaneously,  their  solubilities  increase. 

In  contrast  to  this,  Ashley  and  Mitchell54-55  found 
that  a  considerable  amount  of  Co2+  was  present  in  an  octahe¬ 
dral  environment  (as  would  be  predicted  from  the  intercala¬ 
tion  model),  but  viewed  Mo6+  as  tet rahedrally  coordinated. 
These  authors  also  found  a  synergic  effect  in  the  concurrent 
adsorption  of  Co  and  Mo,  and  a  decrease  in  the  percentage  of 
tetrahedral  Co  with  increasing  Mo  concentration,  interpreted 
as  a  competition  for  tetrahedral  sites. 

A  predominantly  tetrahedral  coordination  of  Mo  has  also 
been  found  using  diffuse  reflectance  spectroscopy  by  Chung 
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and  Massoth56  and  Gajardo  et  al.46.  Chung  and  Massoth  sug¬ 
gested  that  the  predominantly  tetrahedral  or  octahedral 
coordination  of  Mo  depends  on  calcining  temperature  as  well 
as  on  the  amounts  of  both  cobalt  and  molybdenum  in  the  cata¬ 
lyst  . 

Subsequent  spectroscopic  data  by  Mitchell  and  Trifiro57 
also  confirmed  the  presence  of  linked  Mo04  tetrahedra  and 
CoOt  octahedra  in  the  fresh  catalyst.  On  sulfiding,  it  was 
found  that  sulfur  added  to  Mo04,  and  that  no  more  than  one 
or  two  O2-  ions  per  unit  cell  were  replaced  by  S2",  probably 
those  bridging  between  Mo  and  Co.  No  evidence  was  found  for 
discrete  sulfides,  such  as  MoS2  and  Co,S8  ,  and  as  a  result, 
these  authors  concluded  that  the  active  sites  were  probably 
oxysulfides,  rather  than  sulfides.  These  catalysts  were  not 
kept  on  stream  for  a  great  length  of  time,  and  consequently, 
the  further  removal  of  02'  could  be  another  effect  associ¬ 
ated  with  the  change  of  the  catalyst  state  from  monolayer  to 
intercalated. 

A  detailed  series  of  experiments  were  performed  in  1974 
by  De  Beer  and  his  coworkers58  to  test  the  relevance  of  the 
two  models  to  the  preparation  and  operation  of  active  cata¬ 
lysts.  Since  the  intercalation  model  predicted  that  the  most 
effective  promoter  ions  should  be  those  that  prefer  octahe¬ 
dral  sites,  and  the  monolayer  model  predicted  that  they 
should  prefer  a  tetrahedral  coordination,  it  was  hypothe¬ 
sized  that  if  a  certain  mobility  were  assumed,  some 
promoters  would  become  more  active  with  use  and  some  less 
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active,  depending  on  their  initial  locations  and  inherent 
preferences . 

To  test  this,  three  methods  of  preparation  were  used:  a 
sequential  impregnation  of  the  carrier  with  Mo  and  Co  solu¬ 
tions  (to  produce  an  intrinsically  ideal  monolayer  system), 
a  slurrying  of  MoS2  with  the  carrier  before  the  promoter  was 
added  (to  yield  a  true  intercalation  system),  and  a  third 
method  designed  to  produce  a  catalyst  intermediate  between 
the  two. 

Thiophene  hydrodesulfurization  at  atmospheric  pressure 
and  400  °C  was  used  as  a  test  reaction.  The  results  indi¬ 
cated  that  the  standard  method  of  catalyst  preparation 
(sequential  impregnation)  does  result  in  a  monolayer  cata¬ 
lyst,  but  that  under  actual  run  conditions,  the  promoter 
migrates  to  the  surface,  and  within  a  matter  of  hours,  the 
catalyst  becomes  essentially  an  intercalated  system. 

This  study  also  confirmed  two  predictions  of  Farragher 
and  Cossee35.  It  was  shown  that  zinc,  although  highly  active 
initially,  was  no  promoter  for  the  intercalated  system  and 
that  nickel,  less  active  than  cobalt  initially,  becomes  more 
active  after  a  few  hours  on  stream.  This  also  explained  the 
earlier  results  of  De  Beer45,  who  had  found  zinc  and  cobalt 
to  be  much  more  effective  promoters  than  nickel.  In  addi¬ 
tion,  the  prediction  that  the  activities  of  Co-W  and  Ni-W 
catalysts  should  be  nearly  equal  was  also  verified. 

It  was  also  shown  that  catalysts  prepared  to  be 
initially  of  the  monolayer  type  showed  higher  activities 
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after  being  converted  to  intercalation  systems  than  did  cat¬ 
alysts  prepared  to  be  intercalated  initially.  It  was  sug¬ 
gested  that  the  stabilization  of  the  monolayer  by  Co2+  leads 
to  smaller  crystals,  which  are  immediately  intercalated  by 
these  ions  as  they  rise  to  the  surface. 

Schuit  and  Gates  had  suggested  that  that  the  two  models 
may  not  be  mutually  exclusive  in  either  nickel  or  cobalt 
catalysts,  and  that  perhaps  both  systems  are  present  simul¬ 
taneously  with  one  structure  more  active  for  (say)  hydrogen¬ 
ation  and  the  other  more  active  for  desulfurization. 

Behbahany  et  ai.59  have  recently  published  data  sup¬ 
porting  this  idea;  in  this  study,  a  maximum  in  thiophene 
desulfurization  activity  corresponded  to  an  intermediate 
ratio  of  tetrahedral  to  octahedral  cobalt. 

This  two-site  concept  was  not  without  precedent.  Desi- 
kan  and  Amberg5  had  concluded  from  selective  poisoning 
experiments  using  a  Co-Mo-y-Al 203  catalyst  that  two  types  of 
sites  were  present  on  the  catalyst  surface.  The  first  was 
characterized  as  being  strongly  electrophilic,  highly  sensi¬ 
tive  to  poisoning  by  the  adsorption  of  basic  nitrogen  com¬ 
pounds  such  as  pyridine,  and  mainly  responsible  for  olefinic 
hydrogenation  (and  occasionally  dehydrogenation).  The  other 
type  of  site  was  seen  as  less  electrophilic,  more  resistant 
to  poisoning,  and  more  important  for  desulfurization. 
Mitchell45  also  proposed  the  necessity  for  two  types  of 
sites,  one  more  active  for  hydrogenation  and  isomerization 
and  connected  with  the  presence  of  Co2+,  and  the  other  more 
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directly  concerned  with  desulfurization.  Ahuja  et  al.60  also 
concluded  that  there  was  a  hydrogenation/dehydrogenation 
function  as  well  as  a  different  acid  function  operating  in 
sulfided  Group  VIII-Group  Via  catalysts. 

In  1973,  Smith,  Hinckley,  and  Behbahany61  combined  a 
selective  poisoning  technique  with  the  deuter iumolysis  of 
sulfur  and  nitrogen  compounds  to  offer  further  evidence  of 
multiple  sites  on  Co-Mo-y-Al 203 .  These  authors  found  evi¬ 
dence  for  three  types  of  sites:  one  for  low  activation 
energy  multiple  exchange  and  hydrogenation  located  on  the 
y-Al203,  one  for  desulfurization,  and  one  for  a-exchange, 
the  last  two  both  located  on  the  catalyst  itself. 

Hargreaves  and  Ross62  used  the  hydrodesulfurization  of 
thiophene  at  pressures  below  atmospheric  and  temperatures  in 
the  range  200-400°C  to  investigate  the  nature  of  these  two 
major  types  of  sites  in  more  detail.  A  variety  of  Co-Mo- 
y-Al203  formulations  were  used;  it  was  found  that  the  addi¬ 
tion  of  Co  (up  to  a  maximum)  increased  the  rate  of  desulfur¬ 
ization  much  more  than  it  increased  the  rate  of  hydrogena¬ 
tion.  From  this,  it  was  concluded  that  the  desulfurization 
sites  were  dependent  on  promotion  by  Co,  while  the  hydrogen¬ 
ation  sites  probably  were  not.  Since  the  selectivity  of  the 
reaction  did  not  change  while  desulfurization  activity 
increased  when  the  degree  of  sulfiding  was  increased,  it  was 
also  concluded  that  the  hydrogenation  sites  were  probably  in 


the  sulfided  state. 
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In  1976,  Parsons  and  Ternan63  investigated  the  desul¬ 
furization,  deni trogenat ion ,  and  hydrocracking  activities  of 
various  formulations  of  unpromoted  Mo-y-Al203  catalysts  on  a 
gas  oil  in  a  high  pressure  reactor.  The  fact  that  the  extent 
of  cracking  increased  significantly  with  an  increase  in  Mo 
concentration  while  desulfurization  did  not  offered  addi¬ 
tional  evidence  that  these  two  processes  took  place  on  dif¬ 
ferent  sites,  and  that  the  desulfurization  site  was  more 
dependent  on  promotion. 

A  number  of  different  cations  were  tested  as  promoters 
for  this  system.  No  clear  correlation  was  found  between 
ionic  radius  and  activity  in  any  of  the  three  reactions; 
Parsons  and  Ternan  interpreted  this  as  evidence  for  the 
inapplicability  of  the  intercalation  model. 

The  same  year,  De  Beer  and  his  coworkers26  compared  the 
cyclohexane  isomerization  activities  of  Mo-y-Al203  and  W- 
y-Al  2 0 3  catalysts,  both  unpromoted  and  promoted  by  Ni  and 
Co.  The  two  unpromoted  catalysts  showed  a  much  higher  activ¬ 
ity  than  the  promoted  catalysts,  and  the  unpromoted  cata¬ 
lysts  on  a  more  acid  support  were  more  active  than  those  on 
a  less  acidic  base.  De  Beer  suggested  that  isomerization  was 
essentially  a  carrier  phenomenon,  even  though  unsupported 
sulfide  catalysts  had  been  shown  to  possess  some  isomeriza¬ 
tion  activity. 

In  summary,  the  complex  chemistry  of  sulfided  hydrode¬ 
sulfurization  catalysts  has  only  begun  to  be  understood 
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within  the  past  decade.  Both  the  intercalation  and  monolayer 


' 


25 


models  have  received  some  verification  from  studies  of  non- 
reduced,  reduced,  and  sulfided  catalysts  using  various 
analytical  techniques,  and  both  have  been  extended  and 
refined  as  a  result.  Each  model  has  been  successful  in 
predicting  catalyst  behavior  under  some  processing  condi¬ 
tions,  but  the  precise  conditions  under  which  either  (or 
both)  model  may  be  applicable  have  not  yet  been  defined,  and 
appear  to  depend  on  catalyst  composition,  methods  of  prepa¬ 
ration  and  pretreatment,  and  time  on  stream. 

These  studies  have  also  revealed  that  at  least  two 
functionally  different  types  of  catalyst  sites  appear  to  be 
present  on  the  catalyst  surface:  one  more  active  for  desul¬ 
furization  and  another  more  active  for  hydrogenation  and 
dehydrogenation.  It  also  appears  that  a  third  type  of  site 
responsible  for  isomerization  may  be  present.  To  date,  there 
is  no  conclusive  data  linking  a  given  structure  or  configur¬ 
ation  on  the  catalyst  with  a  particular  type  of  catalytic 
activity . 


2.3  Effect  of  Catalyst  Composition 

Most  of  the  published  hydrodesulfurization  data  have 
been  based  on  one  or  more  of  the  Group  VIII-Group  Vla- 
y-Al  20 3  catalysts.  The  Group  VIII  metal  is  usually  cobalt 
and  the  Group  Via  metal  usually  molybdenum,  but  occasionally 
either  nickel  or  tungsten  (or  both)  have  been  used  in  place 
of  these  respective  elements. 
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A  few  studies  have  attempted  to  compare  two  or  more  of 
these  alternative  formulations  under  various  processing  con¬ 
ditions.  It  is  generally  agreed  that  Ni  and  Co  are  the  best 
promoters  for  hydrogenolysis  reactions,  but  considerable 
difference  of  opinion  exists  as  to  which  is  the  better  of 
the  two  for  any  of  these  reactions. 

Beuther,  Flinn,  and  McKinley64  in  1959  were  among  the 
first  to  compare  the  relative  activities  of  Co-Mo  and  Ni-Mo 
catalysts.  Using  the  hydrodesulfurization  of  a  light  gas  oil 
at  4.1  MPa  and  315-370  °C  as  a  test  reaction,  cobalt  was 
found  to  be  only  slightly  more  active  than  nickel,  but  sub¬ 
stantially  more  nickel  than  cobalt  was  required  to  produce 
this  activity  for  a  given  amount  of  molybdenum. 

In  1970,  Ahuja,  Derrien,  and  Le  Page60  again  tried  to 
compare  the  activity  and  selectivity  of  hydrodesulfurization 
catalysts.  In  this  investigation,  Co-Mo,  Ni-Mo,  Co-W,  and 
Ni-W  catalysts  supported  on  silica  and  y-Al203  were  tested 
for  hydrogenation  and  desulfurization  activity  at  350°C  and 
atmospheric  pressure.  Nickel  was  found  to  be  the  more  active 
promoter  for  hydrogenation,  but  no  difference  was  found 
between  the  two  promoters  for  desulfurization,  or  between 
the  catalysts  for  either  hydrogenation  or  desulfurization. 
The  results  of  this  study  have  since  been  challenged  on  the 
basis  that  they  may  have  been  based  on  catalysts  of  unequal 
surface  area  4  5 . 

De  Beer  and  his  coworkers45  compared  a  series  of  Mo- 
y-Al  20 3  catalysts  promoted  by  various  amounts  of  Co,  Ni,  Zn, 
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and  Mn  at  400°C  and  atmospheric  pressure.  The  hydrodesulfur¬ 
ization  of  thiophene  and  the  hydrogenation  of  the  butenes 
produced  in  this  reaction  were  used  to  evaluate  the  perform¬ 
ance  of  the  various  catalysts.  These  authors  found  the  rate 
of  hydrogenation  to  be  relatively  insensitive  to  the  nature 
of  the  promoter  ion,  and  that  Co  and  Zn  appeared  to  have 
significantly  greater  desulfurization  activity  than  Mn  or 
Ni,  although  this  latter  effect  only  appeared  at  higher 
Group  VIII  metal  to  molybdenum  ratios.  The  finding  that  the 
maximum  activity  for  each  metal  occurred  at  a  different 
metal  to  molybdenum  ratio  implied  that  catalyst  activity  was 
not  simply  a  matter  of  filling  tetrahedral  vacancies. 

After  the  emergence  of  the  intercalation  and  monolayer 
models  of  catalyst  structure,  De  Beer58  retested  the  rela¬ 
tive  activities  of  various  catalysts  and  promoters.  A  wide 
variety  of  formulations  were  tested  for  thiophene 
hydrodesulfurization  activity,  including  Ni-Mo-y-Al 203 , 
Zn-Mo-y-Al 203 ,  Co-W-y-Al 20 3 ,  Ni-W-y-Al 203 ,  and  four  Co-Mo- 
y-Al  2 0 3  variations.  All  catalysts  were  used  in  the  reduced 
state,  although  two  Co-Mo-y-Al 203  catalysts  and  one  unpro¬ 
moted  Mo-y-Al 203  catalyst  were  run  in  both  the  reduced  and 
sulfided  states.  (Presulfiding  was  found  to  substantially 
improve  the  activity  of  the  promoted  catalysts,  but  to  sig¬ 
nificantly  impair  the  activity  of  the  unpromoted  catalyst.) 

These  authors  found  nickel  to  be  a  much  more  effective 
promoter  for  the  molybdenum  catalyst  than  cobalt,  while  zinc 
did  not  show  any  promoter  effect  at  all;  little  difference 
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was  found  between  nickel  and  cobalt  as  promoters  for  the 
tungsten  catalyst.  The  Co-Mo-y-Al 20 3  catalyst  was  also  found 
to  be  substantially  more  active  than  the  Co-W-y-Al203  cata¬ 
lyst,  although  the  surface  areas  and  cobalt  contents  of  the 
two  catalysts  were  not  matched  as  carefully  as  in  the  other 
comparisons  of  this  study.  (No  corresponding  Ni-Mo-y-Al 203 
vs.  N i -W-y-Al  2 0 3  test  was  made.) 

The  differences  between  the  results  of  this  study  and 
those  of  some  of  the  previous  studies  concerning  the  rela¬ 
tive  activities  of  the  various  promoters  was  attributed  to  a 
fundamental  change  in  the  chemistry  of  the  catalyst  surface 
during  the  first  few  hours  of  operation:  previous  data  may 
have  been  based  on  a  catalyst  configuration  which  would  be 
only  short-lived  under  actual  processing  conditions. 

In  1975,  Satterfield,  Model 1,  and  Mayer45  investigated 
the  concurrent  hydrodesulfurization  of  thiophene  and  hydro- 
denitrogenat ion  of  pyridine  over  Co-Mo-y-Al 203 ,  Ni-Mo- 
y-Al  20 3  ,  and  Ni-W-y-Al  203  at  atmospheric  pressure  and  250- 
425  °C.  Because  the  concentrations  of  the  active  ingredients 
and  the  surface  areas  per  unit  weight  of  catalyst  were  some¬ 
what  different  for  the  different  samples,  no  direct  compari¬ 
sons  could  be  made,  but  on  a  weight  basis,  the  Ni-Mo  and 
Ni-W  catalysts  were  found  to  be  significantly  more  active 
than  the  Co-Mo  catalyst  for  both  reactions. 

Parsons  and  Ternan43  recompared  the  relative  activities 
of  various  promoters  in  1976,  with  the  intention  of 
developing  catalysts  that  would  be  much  less  expensive  than 
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current  commercial  catalysts  yet  still  be  reasonably  active, 
for  use  with  feeds  with  high  metals  contents  (and  thus  high 
poisoning  rates).  A  heavy  gas  oil  derived  from  Athabasca 
bitumen  was  treated  in  a  flow  reactor  at  400°C  and  14  MPa, 
and  analyzed  for  desulfurization,  den i t rogenat ion ,  and 
hydrocracking  activity.  Co  and  Ni  were  found  to  be  the  most 
active  promoters  in  all  cases;  the  two  were  virtually  iden¬ 
tical  in  desulfurization  activity,  while  Ni  appeared 
slightly  more  active  for  denitrogenat ion  and  cracking. 

That  same  year,  De  Beer,  Dahlman,  and  Smeets26  compared 
the  high  pressure  desulfurization  and  hydrogenation  activi¬ 
ties  of  W-y-Al 203  and  Mo-y-Al203  catalysts  promoted  by  both 
Ni  and  Co.  Because  a  number  of  important  parameters  were 
varied  in  the  preparation  of  the  catalyst,  only  the  main 
effects  were  reported.  The  results  of  thiophene  desulfuriza¬ 
tion  were  in  agreement  with  De  Beer's  earlier58  data.  Also, 
the  promoter  effects  of  Co  and  Ni  for  desulfurization  were 
much  higher  for  benzothiophene  than  for  thiophene. 

In  contrast  to  the  data  of  Hargreaves  and  Ross62,  both 
Co  and  Ni  seemed  to  be  relatively  good  promoters  for  benzene 
hydrogenation.  For  this  reaction,  Ni  seemed  the  better  pro¬ 
moter  and  W  the  better  catalyst,  as  had  been  found  by  Ahuja. 
This  result  was  not  in  accord  with  the  data  of  Farragher  and 
Cossee39,  who  found  Co  to  be  a  much  more  effective  promoter 
than  Ni  for  benzene  hydrogenation  on  W  catalysts. 

Ripperger  and  Saum20  tested  Co-Mo-y-Al 203  and  Ni-Mo- 
y-Al  2 0 3  catalysts  for  desulfurization,  hydrogenation,  and 
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deni trogenat ion  activity  at  medium  pressure  in  a  trickle- 
flow  reactor  with  a  solution  of  thiophene  and  pyridine  in 
isooctane  as  the  feed.  It  was  found  that  for  a  given  degree 
of  desulfurization,  the  Ni  catalyst  possessed  a  greater 
activity  than  the  Co  catalyst  for  both  hydrogenation  and 
denitrogenat ion .  In  addition,  the  deni trogenation  activity 
of  the  Ni  catalyst  (and  to  a  lesser  degree,  the  Co  catalyst) 
could  be  significantly  increased  by  treatment  with  phos¬ 
phoric  acid.  It  was  suggested  that  the  denitrogenation 
occurs  on  the  nickel  surface,  and  that  phosphoric  acid  hin¬ 
ders  migration  of  the  nickel  ions  into  the  alumina  lattice. 
This  effect  was  also  shown  for  desulfurization,  although  the 
increase  in  activity  was  much  smaller. 

In  1978,  Furimsky,  Ranganathan,  and  Parsons1’  compared 
several  Ni-  and  Co-promoted  catalysts  for  hydrodeni trogena¬ 
tion  activity  by  treating  two  bitumen  distillates  in  a  flow 
reactor  at  14  MPa  and  320-420°C.  No  difference  between  the 
two  promoters  was  found;  the  apparent  contradiction  between 
this  result  and  those  of  previous  investigators  was 
explained  as  possibly  due  to  the  inapplicability  of  low- 
pressure  data  to  high-pressure  industrial  processes,  or  to 
differences  in  the  preparation  of  the  two  catalysts  in  the 
various  studies. 

Recently  Furimsky66  retested  the  hydrogenolysis  activi¬ 
ties  of  Co-Mo-y-Al 203  and  Ni-Mo-y-Al 203  in  a  similar  reactor 
to  the  one  used  previously.  Furimsky  again  found  little 
difference  in  desulfurization,  denitrogenation,  or 
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deoxygenation  activity  between  the  two  catalysts,  and 
reiterated  his  earlier  explanations  for  the  discrepancies 
between  his  results  and  those  of  other  researchers. 

The  only  evaluation  of  the  effectiveness  of  various 
supports  found  in  the  literature  is  that  of  De  Beer  et 
al.87.  Based  on  the  activity  of  supported  Co-Mo  for  thio¬ 
phene  desulfurization,  it  was  found  that  Si02 ,  y-Al203,  and 
77-Al  203  were  all  good  supports  as  long  as  the  specific  sur¬ 
face  areas  were  high,  but  that  A1203  was  preferred  to  Si02 
because  it  inhibited  the  formation  of  C0M0O4. 

In  summary,  the  relative  activities  of  various  catalyst 
formulations  are  still  unclear.  Cobalt  and  nickel  appear  to 
be  the  best  promoters  for  hydrodesulfurization,  molybdenum 
the  best  catalyst,  and  y-alumina  the  best  support.  The 
advantage  of  these  materials  over  others  such  as  iron, 
tungsten,  and  silica  do  not  seem  to  be  large  however,  and 
cost  may  become  the  determining  factor  for  applications 
where  the  rate  of  fouling  is  great. 

2.4  Effect  of  Catalyst  Preparation 

The  most  detailed  study  concerning  the  optimum  methods 
of  preparation  of  Ni-Mo-y-Al 203  catalysts  available  at  pre¬ 
sent  is  that  of  Laine,  Pratt,  and  Trimm.  Both  the  concentra¬ 
tion  of  the  active  metals18  and  the  method  of  preparation 8 8 


were  varied  in  this  study. 
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The  best  method  of  preparation  was  found  to  be  the 
sequential  impregnation  of  first  Ni  and  then  Mo,  a  result 
which  would  support  De  Beer's58  hypothesis  that  catalysts 
initially  approximating  the  monolayer  model  would  produce  a 
better  catalyst  even  after  that  structure  was  no  longer  in 
existence . 

The  best  starting  salt  was  found  to  be  nickel  acetate, 
despite  the  fact  that  it  was  only  partially  soluble  in  the 
liquid  used  for  impregnation.  It  was  suggested  that  a 
favorable  anion  should  interfere  with  the  adsorption  of 
ammonia,  and  should  form  an  ammonia  salt  which  would  easily 
escape  from  the  catalyst  during  the  calcining  period.  This 
was  supported  by  the  finding  that  catalysts  prepared  from 
nickel  hydroxide,  nickel  phosphate,  and  nickel  sulfate  were 
the  least  active.  While  ammonium  acetate  and  ammonium 
nitrate  would  not  be  removed  more  easily  than  the  ammonia, 
ammonium  phosphate  and  ammonium  sulfate  produced  from  the 
unsuccessful  starting  compounds,  decomposition  of  the 
nitrate  and  acetate  would  be  expected  to  lead  to  relatively 
innocuous  anions  which  could  then  be  removed  as,  for 
example,  carbon  oxides  during  calcination. 

The  optimum  composition  of  the  catalyst  was  found  to  be 
3%  NiO  and  35%  Mo03,  although  Mo03  contents  of  10-15%  were 
only  slightly  less  active.  Below  about  2%  NiO,  the  activity 
of  the  catalyst  decreased  considerably;  this  was  explained 
as  being  due  to  the  formation  of  the  inert  NiMoO*. 
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2.5  Effect  of  Pretreatment  Methods 

Calcining,  prereduction,  and  presulfiding  are  all  known 
to  produce  higher  activity  and  longer  life  in  hydrodesulfur¬ 
ization  catalysts69.  In  addition,  there  may  be  a  relation¬ 
ship  between  coke  formation  on  the  catalyst  and  catalyst 
sulfur  content41.  Despite  this,  little  data  are  available  in 
the  literature  concerning  the  optimum  method  of  pretreatment 
of  these  catalysts. 

Calcining  has  been  the  least  studied  of  the  three  pre¬ 
treatment  steps,  however  the  available  data  on  the  effects 
of  varying  calcination  temperature  are  in  good  agreement. 

In  1977,  Ripperger  and  Saum20tested  the  effects  of 
calcining  temperature  on  the  activities  of  Co-Mo-y-Al 203  and 

i 

Ni-Mo-y-Al 203  catalysts.  The  results  of  this  study  showed 
catalyst  activity  to  be  significantly  impaired  at  calcining 
temperatures  above  550°C. 

This  result  was  supported  in  1979  by  Laine,  Pratt,  and 
Trimm68,  who  used  the  hydrodesulfurization  of  thiophene  over 
a  Ni-Mo-y-Al 203  catalyst  to  test  the  effects  of  various  pre¬ 
treatment  methods  on  catalyst  activity.  An  optimum  calcina¬ 
tion  temperature  was  discovered  at  about  500°C,  and  a  slow 
rate  of  heating  of  the  catalyst  to  the  maximum  calcination 
temperature  was  found  to  be  necessary  to  avoid  a  significant 
loss  of  activity.  Based  on  X-ray  diffraction  studies,  the 
authors  hypothesized  that  above  500°C,  loss  of  catalyst 
activity  through  the  formation  of  Al2(Mo04)3  was  the  most 
important  factor,  rather  than  sintering.  The  lower  activity 
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at  lower  calcination  temperatures  was  explained  as  due  to 
insufficient  reaction  between  the  various  inorganic  com¬ 
pounds  present  on  the  surface  of  the  catalyst.  For  example, 
heating  the  catalyst  too  quickly  could  result  in  vapor  phase 
transport  of  intermediates  between  the  uncalcined  and 
calcined  phases:  if  decomposition  of  nickel  nitrate  were 
incomplete,  it  would  be  expected  to  boil,  resulting  in  a 
loss  of  nickel  in  the  final  catalyst.  This  scenario  was 
verified  by  X-ray  fluorescence  analysis.  Adsorption  of 
ammonia  onto  acidic  sites  was  also  cited  as  a  possible 
factor  to  explain  the  differences  in  activity  produced  by 
various  calcining  procedures. 

Prereduction  has  been  studied  to  a  much  greater  degree 
than  calcining,  but  the  results  are  not  quite  so  conclusive. 
This  may  be  partly  due  to  the  more  complex  nature  of  the 
chemical  changes  associated  with  prereduction,  and  partly 
because  the  effectiveness  of  a  given  prereduction  tempera¬ 
ture  may  depend  more  on  the  conditions  seen  by  the  catalyst 
after  the  prereduction  has  finished,  i.e.  the  actual  reac¬ 
tion  conditions  and  whether  or  not  the  catalyst  is  presul¬ 
fided  prior  to  use. 

Ripperger  and  Saum20  found  that  no  appreciable  reduc¬ 
tion  occurred  below  250°C  for  Mo03-y-Al203  or  Co-Mo-y-Al 203 
catalysts,  and  that  above  500°C  excessive  oxygen  was  lost 
from  the  catalyst,  resulting  in  a  loss  of  activity.  It  was 
suggested  that  there  was  an  optimum  degree  of  reduction  for 
catalysts  used  in  the  sulfided  state:  excess  oxygen  should 
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be  removed  from  the  catalyst,  but  enough  should  remain  to 
allow  formation  of  oxysulfides  rather  than  the  cataly t ically 
less  active  metal  sulfides. 

Laine,  Pratt,  and  Trimm1 8  also  found  that  the  lowest 
prereduction  temperature  used  (300°C)  produced  the  most 
active  sulfided  catalyst,  and  that  a  high  activity  was 
obtained  even  though  no  prereduction  was  carried  out.  These 
authors  suggested  that  a  reduction  of  Mo03  to  Mo02  took 
place  at  above  about  350°C,  which  was  undesirable  because 
the  latter  was  more  difficult  to  sulfide. 

In  1977,  Massoth  and  Kibby52  used  a  gravimetric  tech¬ 
nique  to  test  the  effects  of  various  catalyst  pretreatment 
methods  on  the  rate  of  thiophene  hydrodesulfurization  over 
Mo-y-Al203.  These  authors  also  found  that  prereduced  cata¬ 
lysts  sulfided  to  a  lesser  extent  than  those  which  had  not 
been  prereduced,  but  no  correlation  with  catalyst  activity 
was  found.  It  was  also  found  in  this  study  that  for  a  given 
prereduction  temperature,  catalysts  used  in  the  reduced 
state  showed  a  higher  activity  than  sulfided  catalysts, 
although  little  difference  was  found  between  the  activities 
of  the  optimum  prereduced  catalysts  and  the  optimum  (non- 
prereduced)  sulfided  catalysts. 

Massoth  and  Kibby  hypothesized  that  sulfur  adsorbed 
onto  the  catalyst  in  two  ways  during  the  presulfiding  step: 
by  the  replacement  of  oxygen  atoms  with  sulfur  atoms, 
producing  extra  anion  vacancies  and  thus  more  active  sites, 
and  by  the  irreversible  adsorption  of  inert  H2S  onto  the 
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active  sites.  At  low  prereduction  temperatures,  sulfur  was 
added  mainly  by  exchange  with  oxygen,  while  at  higher  tem¬ 
peratures,  adsorption  of  H2S  was  more  important. 

In  1980,  Gissy,  Bartsch,  and  Tanielan27  used  the  con¬ 
version  of  benzothiophene  over  a  Co-Mo-y-Al 203  catalyst  to 
test  the  effect  of  pretreatment  temperature  on  activity. 
Sulfiding  was  found  to  markedly  increase  the  activity  of  all 
catalysts  studied.  A  mild  degree  of  prereduction  was  found 
to  be  beneficial;  250°C  was  found  to  be  inadequate  to 
produce  this  effect,  but  400°  was  an  effective  temperature 
for  prereduction  prior  to  presulfiding. 

Presulfiding  methods  have  been  studied  to  an  even 
greater  degree  than  those  for  prereduction.  The  results  are 
less  conclusive  still,  and  are  complicated  by  the  fact  that 
the  H2S  produced  in  the  test  reaction  tends  to  affect  the 
level  of  sulfur  on  the  catalyst  and  thus  obscure  the  effects 
of  the  pretreatment. 

Ternan  and  Whalley70,  using  a  Ni-Mo-y-Al 203  catalyst  in 
a  bottom-feed  reactor  operating  on  Athabasca  bitumen  at  high 
pressure,  presulfided  with  mixtures  of  hydrogen  and  either 
feedstock,  carbon  disulfide,  or  hydrogen  sulfide  at  various 
temperatures  and  pressures  to  determine  the  effects  of  pre¬ 
sulfiding  conditions  on  catalyst  performance.  These  authors 
found  that  after  8  hours,  the  sulfur  content  of  the  catalyst 
approached  an  asymptotic  value,  regardless  of  the  presulfid¬ 
ing  conditions.  Although  there  were  initial  differences  in 
catalyst  activity,  no  differences  were  found  after  the  8 
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hour  break-in  period  in  any  of  the  properties  of  the  liquid 
hydrocarbon  product.  The  fraction  of  coke  in  the  catalyst 
was  also  found  to  be  independent  of  presulfiding  conditions. 

In  1976,  De  Beer  et  al.71  tested  the  effects  of  various 
presulfiding  conditions  on  the  sulfur  contents  of  various 
catalysts,  including  Co-Mo-y-Al 203 .  At  a  presulfiding  tem¬ 
perature  of  400°C,  catalyst  sulfur  content  increased  with 
time,  levelling  off  after  about  2  hours.  Higher  presulfiding 
temperatures  produced  higher  sulfur  contents  in  the  range 
250-500°C,  and  in  general,  higher  sulfur  contents  correlated 
with  higher  hydrodesulfurization  activities,  based  on  con¬ 
version  of  thiophene. 

Composition  of  the  presulfiding  mixture  was  also 
tested;  an  increase  in  H2S:H2  ratio  from  0.04  to  1.0  pro¬ 
duced  only  a  10%  increase  in  sulfur  content. 

It  was  also  found  that  as  much  as  30%  of  the  sulfur 
present  after  presulfiding  could  be  removed  within  2  hours 
by  subsequent  reduction.  This  surprisingly  had  no  effect  on 
hydrodesulfurization  activity,  and  it  was  concluded  that  a 
significant  amount  of  "mobile  sulfur"  was  present  which  had 
no  involvement  in  the  desulfurization  activity  of  the  cata¬ 
lyst  . 

The  next  year,  Ripperger  and  Saum20,  working  with  both 
Co-Mo-y-Al 203  and  Ni-Mo-y-Al 203  catalysts  also  used  mixtures 
of  either  carbon  disulfide  or  hydrogen  sulfide  as  presulfid¬ 
ing  agents.  In  this  study,  it  was  found  that  presulfiding 
should  take  place  at  about  250°C  for  optimum  performance.  It 
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was  suggested  that  complete  reduction  of  the  oxides  before 
conversion  to  oxysulfides  was  the  main  reason  for  the  loss 
of  activity  when  prereduction  and  presulfiding  took  place  at 
higher  temperatures. 

Gissy,  Bartsch,  and  Tanielan27  in  1980  also  found  250°C 
to  be  the  optimum  temperature  for  presulfiding,  and  that  a 
high  catalyst  sulfur  content  at  steady  state  corresponded  to 
a  high  catalyst  activity.  These  authors  also  saw  presulfid¬ 
ing  as  an  important  step  in  protecting  the  catalyst  against 
coking . 

Laine,  Pratt,  and  Trimm1 8  also  tested  presulfiding 
methods  in  their  thiophene  hydrodesulfurization  studies. 
While  standard  presulfiding  techniques  involve  a  mixture  of 
10-15%  H2S  in  in  H2,  Laine  et  al.  found  the  most  effective 
presulfiding  mixture  to  be  pure  H2S,  and  the  most  effective 
temperature  (when  using  pure  H2S)  to  be  400°C. 

In  conclusion,  the  method  of  catalyst  pretreatment 
appears  to  have  a  significant  effect  on  catalyst  activity. 
Calcination  should  take  place  at  about  500°C,  followed  by 
prereduction  in  hydrogen  at  about  300°C  and  presulfiding, 
for  which  the  data  on  the  optimum  method  is  inconclusive. 


2.6  On-stream  Catalyst  Behavior 

Catalyst  behavior  over  long  periods  of  time  has  been 
one  of  the  least  studied  aspects  of  the  hydrodesulfurization 
It  has  been  found  that  the  activities  of  prereduced 
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catalysts  increase  from  near  zero  and  asymptotically 
approach  much  higher  values  during  the  first  few  hours  on 
stream,  while  presulfided  catalysts  start  off  much  more 
active,  and  level  off  at  a  much  lower  activity  after  the 
first  few  hours  on  stream18'27'71. 

De  Beer71  suggests  that  the  increase  in  activity  of 
prereduced  catalysts  appears  to  be  due  to  sulfiding  by  H2S 
produced  in  the  reaction,  while  the  decrease  in  activity  of 
sulfided  catalysts  is  due  to  the  establishment  of  a  new 
oxidation/reduction  equilibrium  on  the  catalyst.  Although  a 
significant  loss  of  sulfur  from  the  fresh  catalyst  and  a 
correlation  between  sulfur  content  and  activity  were 
observed,  De  Beer  hesitated  to  attribute  the  initial  loss  in 
activity  to  a  loss  of  sulfur.  It  was  suggested  that  the  high 
activities  found  in  catalysts  which  had  not  been  presulfided 
would  have  required  unreasonably  high  rates  of  on-stream 
sulf idat ion . 

After  a  catalyst  has  been  on  stream  for  4  or  5  hours, 
most  researchers  have  assumed  the  activity  to  be  constant, 
and  there  is  little  data  available  in  the  literature 
describing  catalyst  behavior  beyond  this  point. 

Rollman 2  9 ,  working  at  high  pressure,  recorded  the 
activities  of  his  catalysts  for  dibenzothiophene  hydrodesul¬ 
furization  over  a  period  of  several  days,  but  did  not  detect 
any  change  in  activity  during  this  time. 

Aitken  et  al.72,  working  with  a  Co-Mo-y-Al 203  catalyst 
in  a  trickle-bed  reactor  did  find  a  linear  drop  in  activity 
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over  a  period  of  hundreds  of  hours,  with  the  rate  increasing 
exponentially  with  a  decrease  in  hydrogen  partial  pressure; 
this  was  attributed  to  coking  via  dehydrogenation  reactions 
on  the  catalyst  surface. 

The  only  study  to  date  devoted  to  monitoring  the 
changes  in  hydrodesulfurization  catalysts  over  long  periods 
of  time  is  that  of  Broderick,  Schuit,  and  Gates73.  In  this 
study,  the  hydrodesulfurization  of  dibenzothiophene  over  a 
Co-Mo-y-Al 203  catalyst  was  studied  in  a  liquid  phase  flow 
reactor  at  10  MPa  and  300°C,  with  space  velocity  and  amount 
of  H2S  added  to  the  feed  the  independent  variables.  Typic¬ 
ally,  catalyst  activity  increased  by  about  10%  over  the 
first  10-20  hours  of  operation,  after  which  point  no  change 
was  observed  up  to  150  hours  of  operation. 

From  the  results  of  step  changes  in  the  reaction  condi¬ 
tions,  the  authors  concluded  that  at  low  H2S  concentrations 
in  the  reactor,  a  decrease  in  H2S  concentration  resulted  in 
an  irreversible  loss  of  catalytic  sites,  while  at  higher  H2S 
concentrations,  changes  in  this  parameter  affected  the  reac¬ 
tion  only  by  the  competitive  (reversible)  adsorptions  of 
dibenzothiophene  and  H2S. 

In  1980,  Kilanowski  and  Gates32  reported  a  loss  in  con¬ 
version  of  benzothiophene  with  time  in  their  flow  reactor 
operating  at  atmospheric  pressure  and  low  conversions.  These 
authors  agreed  with  the  hypothesis  of  Broderick  et  al.73, 
that  the  loss  in  activity  was  attributable  to  structural 
changes  in  the  catalyst  associated  with  sulfur  depletion.  To 
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counter  this  loss  in  activity,  approximately  2%  H2S  was 
added  to  their  feed;  this  proved  successful  in  stabilizing 
the  catalyst. 

Also  in  1980,  Gissy,  Bartsch,  and  Tanielan74  reported  a 
similar  loss  in  activity,  again  desulfurizing  benzothiophene 
over  Co-Mo-y-Al 203  in  a  flow  reactor  at  atmospheric  pres¬ 
sure.  It  was  also  found  that  the  presence  of  H2S  eliminated 
this  deactivation,  but  Gissy  et  al.  concluded  that  the  loss 
in  activity  was  due  to  coke  formation  from  the  cracking  of 
n-dodecane,  the  solvent  used  for  benzothiophene.  (No  solvent 
had  been  used  in  the  Kilanowski  and  Gates  study.)  It  was 
suggested  that  the  presence  of  H2S  improves  the  hydrogena¬ 
tion  activity  of  the  catalyst,  and  in  this  way  interfered 
with  the  coking  process.  Substitution  of  n-hexane  for 
n-dodecane  was  found  to  eliminate  the  deactivation;  on  the 
other  hand,  substitution  of  n-dodecene  markedly  increased 
the  rate  of  loss  of  activity. 
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3.  DESCRIPTION  OF  EQUIPMENT 


The  basic  design  of  the  experimental  apparatus  was  that 
of  a  fixed  bed  reactor  with  on-line  gas  sampling,  similar  to 
that  used  by  Satterfield65.  A  flow  diagram  of  the  equipment 
is  shown  in  Figure  3. 

The  liquid  feed  to  the  reactor  was  a  30%  (wt.)  solution 
of  benzothi ophene  (Aldrich  Chemical,  specified  95%  pure, 
although  a  chromatographic  analysis  showed  the  purity  to  be 
approximately  99.8%)  in  n-hexane  (Fisher  Chemicals,  99% 
pure).  Gases  fed  to  the  reactor  were  hydrogen  (99.9995%), 
nitrogen  (99.998%),  and  hydrogen  sulfide  (99.6%),  and  the 
carrier  gas  for  the  gas  chromatograph  was  helium  (99.9999%), 
all  supplied  by  Union  Carbide. 

Trace  oxygen  in  the  hydrogen  and  nitrogen  streams  was 
converted  to  water  by  a  Matheson  OR- 1 0  catalytic  converter. 
The  water  was  subsequently  removed  by  passing  the  gas  stream 
through  a  glass  tube  22  mm  in  diameter  packed  with  10  cm  of 
1/16"  pellets  of  Molecular  Sieve  Type  4A.  A  2.5  cm  layer  of 
Drierite  (indicating  calcium  sulfate)  was  located  above  and 
below  the  bed  to  indicate  when  the  molecular  sieve  needed 
replacing.  The  Drierite  and  molecular  sieve  were  held  in 
place  with  plugs  of  glass  wool.  Trace  water  in  the  helium 
was  removed  by  passing  the  gas  through  a  steel  tube  1/2"  in 
diameter  containing  a  35  cm  bed  of  molecular  sieve. 

Gas  flow  rates  were  monitored  with  Matheson  Type  600 
rotameters  and  controlled  by  needle  valves.  The  hydrogen  and 
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Figure  3  -  Equipment  Flow  Diagram 
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nitrogen  rotameters  were  calibrated  using  a  stopwatch  and  a 
soap  bubble  flowmeter  (the  calibration  data  are  given  in 
Appendix  1),  while  the  manufacturer's  supplied  calibration 
curve  was  used  for  the  hydrogen  sulfide  rotameter,  with  a 
correction  made  for  the  difference  in  molecular  weights 
between  the  calibration  fluid  (air)  and  hydrogen  sulfide. 

Liquid  feed  rate  was  controlled  with  a  Sage  Instruments 
Model  355  syringe  pump  fitted  with  a  50  ml  glass  syringe; 
calibration  data  for  the  syringe  pump  are  given  in  Appendix 
2.  The  liquid  was  vaporized  in  the  gas  stream  in  an  insu¬ 
lated  box  heated  to  300°C  by  a  strap  heater.  A  mercury- 
filled  barometer  used  to  measure  atmospheric  pressure  and  a 
manometer  filled  with  Meriam  fluid  (density  1.75)  located 
upstream  of  the  hot  box  were  used  to  measure  the  gauge  pres¬ 
sure  of  the  gas  stream. 

A  100  ml  cylinder  was  located  downstream  of  the  syringe 
pump  to  eliminate  any  concentration  gradients  in  the  feed  to 
the  reactor  which  might  have  been  caused  by  irregular  vapor¬ 
ization  of  the  liquid  feed.  Preliminary  tests  had  shown  that 
when  the  liquid  feed  rate  was  too  high,  droplets  were  formed 
on  the  end  of  the  syringe  tubing  which  flashed  quickly  when 
they  fell  onto  the  wall  of  the  flowing  gas  stream  tubing, 
rather  than  vaporizing  continuously  at  the  end  of  the 
syringe  tubing.  This  caused  oscillations  in  the  pressure  in 
the  tubing  with  a  magnitude  of  about  0.25  kPa,  measured 
using  the  manometer.  Because  no  oscillations  were  observed 
during  the  actual  experimentation,  it  was  concluded  that  any 
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concentration  gradients  in  the  gas  stream  entering  this 
100  ml  cylinder  were  very  small  and  would  have  been 
eliminated  in  this  tank,  although  the  mean  residence  time 
was  only  8  s  at  the  highest  flow  rate  used. 

The  reactor  consisted  of  a  0.1000  g  (±  0.0004  g)  bed  of 
catalyst  (depth  approximately  5  mm)  contained  in  1/2”  (nom.) 
stainless  steel  tubing.  The  bed  was  supported  by  a  100  mesh 
stainless  steel  screen,  which  was  in  turn  supported  by  a  bed 
of  1.0  mm  glass  beads.  The  bottom  of  the  tubing  was  plugged 
with  a  small  piece  of  glass  wool,  to  keep  the  beads  in  place 
while  the  reactor  was  operating.  The  catalyst  bed  was 
covered  by  a  5  mm  layer  of  -40+60  mesh  glass  beads,  to 
produce  a  more  constant  velocity  profile  of  the  reactant  gas 
stream  when  it  came  into  contact  with  the  catalyst.  This  was 
covered  by  a  second  5  mm  layer  of  1.0  mm  glass  beads  to  hold 
the  smaller  beads  in  place.  Both  the  glass  beads  and  the 
stainless  steel  screen  were  verified  to  be  cataly t ically 
inert  prior  to  experimentation. 

The  catalyst  used  was  a  commercial  catalyst,  Nalco 
NM  502,  containing  3%  Ni  and  14%  Mo  (wt.)  supported  on 
y-alumina,  ground  in  the  laboratory  and  sieved  to  -40+50 
mesh.  The  surface  area11  of  the  catalyst  was  270  m2/q,  and 
the  pore  volume  was  0.55  cm3/g.  The  pellet  density  was 
680-780  kg/m3. 

The  reactor  was  held  at  a  constant  temperature  in  a 
fluidized  bed  of  90  mesh  Alundum,  which  was  heated 
externally  with  nichrome  wire.  The  fluidized  bed  was 
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approximately  10  cm  in  diameter  and  30  cm  deep  when  fluid¬ 
ized;  the  catalyst  bed  was  located  approximately  10  cm  below 
the  surface  of  the  fluidized  bed.  The  top  of  the  fluidized 
bed  was  contained  by  approximately  6  cm  of  1/4"  glass 
Raschig  rings,  supported  on  a  coarse  wire  screen.  The  tem¬ 
perature  of  the  fluidized  bed  was  controlled  by  a  Techne 
TC4A  temperature  controller  to  within  2°C  and  monitored  by  a 
Type  J  ( iron-constantan )  thermocouple  located  close  to  the 
reactor  tubing  at  about  the  same  height  as  the  catalyst  bed. 
Temperature  in  the  reactor  was  measured  using  a  Type  J  ther¬ 
mocouple  inserted  into  the  catalyst  bed.  (The  temperature  of 
the  catalyst  bed  did  not  differ  significantly  from  that  of 
the  fluidized  bed  at  any  time.)  Pressure  in  the  reactor  was 
held  at  101  kPa  (±  3  kPa)  using  a  globe  valve  located  in  the 
hot  box  downstream  of  the  reactor.  The  pressure  in  the  reac¬ 
tor  was  assumed  not  to  differ  significantly  from  that  of  the 
gas  stream  as  measured  by  the  manometer  upstream  of  the 
syringe  feeder. 

The  products  were  sampled  with  a  Valeo  model  V-4-HT  gas 
sampling  valve  (internal  volume  3  fil)  located  in  an  insu¬ 
lated  steel  box  mounted  on  the  side  of  the  gas  chromatograph 
and  heated  to  approximately  250°C  by  a  strap  heater.  Analy¬ 
sis  was  performed  using  a  dual  column  F  &  M  720  gas  chroma¬ 
tograph;  the  conditions  under  which  this  analysis  was  per¬ 
formed  are  summarized  in  Table  1.  The  sampling  valve,  injec¬ 
tion  port,  and  detector  temperatures  shown  are  approximate, 
and  varied  by  about  ±5°  over  the  course  of  this  study.  The 
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Table  1  -  Summary  of  Chromatographic  Analysis  Conditions 


Column:  Length:  10T 

Diameter:  1/8" 

Packing:  Dexsil  300  on  80/100  Chromosorb  W/AW 

Carrier  gas:  Helium 

Flow  Rate:  40  cm3/nun 

Detector  Type:  Thermal  Conductivity 
Current:  150  ma 

Temperatures:  Gas  Sampling  Valve:  260°C 
Injection  Port:  200°C 
Detector:  260°C 


Temperature  Programming:  Rate:  10°/min 
Initial  Temperature:  40°C 
Final  Temperature:  190°C 
Total  Analysis  Time:  15  min 
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flow  of  carrier  gas  was  checked  regularly  with  a  soap  bubble 
flowmeter.  Output  from  the  gas  chromatograph  was  integrated 
and  recorded  using  an  HP  3380A  integrator. 

To  illustrate  the  degree  of  separation  of  the  major 
components  of  the  product  stream,  a  typical  chromatogram  is 
shown  in  Figure  4.  (Note  that  the  peak  heights  have  been 
plotted  on  a  logarithmic  scale,  and  as  a  result  the  relative 
peak  areas  may  be  deceptive.) 

Liquid  products  were  collected  in  a  glass  condenser 
held  at  room  temperature  (approximately  22cC)  to  obtain  a 
concentrated  sample  of  the  heavier  products  for  qualitative 
analysis.  Non-condensibles  were  vented  to  a  fume  hood.  No 
mass  balances  were  computed  for  either  the  liquid  products 
or  the  gas  stream  leaving  the  condenser. 

All  equipment  in  direct  contact  with  the  reactants  or 
products  was  constructed  of  Type  316  stainless  steel  unless 
noted  otherwise.  All  equipment  was  connected  with  1/4" 
tubing,  except  the  line  from  the  syringe  pump  to  the  flowing 
gas  stream  and  the  lines  connecting  the  gas  sampling  valve 
to  the  columns,  all  of  which  were  1/16".  All  lines  between 
the  reactor,  hot  box,  sampling  valve,  and  condenser  were 
electrically  traced  with  "Heats-by-the-Yard"  and  insulated. 
The  temperatures  of  these  lines  were  monitored  with  Type  J 
thermocouples  inserted  under  the  tracing.  The  voltages 
across  each  strap  heater  and  section  of  tracing  were 
controlled  manually  using  Variac  variable  voltage 
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RT 
1  .07 
1.80 
5.75 
6.64 
12.00 
13.17 


TYPE 

M 


AREA 

37373 

37192 

2201 

283 

14229 

195 


AREA  % 
40.86 
40.66 
2.406 
.3094 
15.56 
.213  2 


HP  3380A 
DLY  OFF 
MV/M  .03 


STOP  15 
ATTN  LOG 


REJECT  OFF 


Figure  4  -  Sample  Chromatogram 
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4.  EXPERIMENTAL  PROCEDURE 


4.1  Catalyst  Pretreatment 

No  tests  were  done  in  the  present  study  to  assess  the 
effects  of  various  pretreatment  methods.  Instead,  it  was 
decided  to  use  the  basic  method  of  Ripperger  and  Saum2 0 :  5 
hours  calcining  in  air  at  500°C,  followed  by  4  hours  prere¬ 
duction  at  250°,  followed  by  4  hours  presulfiding  at  250° 
with  a  mixture  of  10%  H2S  in  H2 . 

The  catalyst  was  ground,  sieved  and  graded  as  to  parti¬ 
cle  size,  and  calcined  in  air  at  500°C  for  5  hours.  Follow¬ 
ing  the  suggestions  of  Laine  et  al.68,  the  catalyst  was 
heated  to  the  maximum  calcining  temperature  as  slowly  as 
possible,  about  2°/min.  This  heatup  time  was  not  included  in 
the  5  hours  calcining  time.  After  calcining,  the  catalyst 
was  stored  in  a  desiccator  until  required. 

When  a  new  catalyst  bed  was  required,  the  fluidized  bed 
was  cooled  to  250°C  and  the  system  was  flushed  with  nitro¬ 
gen.  The  reactor  tube  was  cleaned  and  packed  with  the  glass 
beads  and  catalyst,  then  transferred  to  the  fluidized  bed 
and  connected  to  the  rest  of  the  equipment.  The  nitrogen 
flow  was  terminated,  the  hydrogen  flow  was  set  at  50 
cm3 /min ,  and  the  pressure  was  adjusted  to  101  kPa.  After  4 
hours,  the  hydrogen  sulfide  flow  was  set  at  5  cm3/min,  while 
the  hydrogen  flow  was  left  unchanged.  After  another  4  hours, 
the  hydrogen  sulfide  flow  was  terminated,  and  the  reactor 
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operating  conditions  were  set  at  the  desired  values. 

4.2  Reactor  Operation 

When  the  pretreatment  had  ended,  the  fluidized  bed  tem¬ 
perature  controller  was  reset  to  the  desired  reaction  tem¬ 
perature,  the  syringe  pump  was  started,  the  hydrogen  flow 
rate  was  reset,  and  the  pressure  was  adjusted.  This  was  con¬ 
sidered  the  start  of  reactor  operation. 

After  the  reactor  temperature  had  risen  to  the  desired 
value,  (the  temperature  of  the  fluidized  bed  could  be  raised 
by  about  10°/min)  and  a  stabilization  period  of  about  15 
minutes  had  passed,  gas  samples  from  the  product  stream  were 
taken  at  varying  intervals  and  analyzed  for  the  next  25  to 
30  hours,  after  which  the  run  was  terminated  and  the  cata¬ 
lyst  discarded. 

The  syringe  usually  had  to  be  refilled  at  least  once 
during  a  run.  When  this  became  necessary,  the  gas  flow  was 
terminated  and  the  syringe  was  disconnected  from  the  rest  of 
the  apparatus  and  refilled.  This  usually  took  about  5  min¬ 
utes,  during  which  time  the  connection  between  the  syringe 
and  the  gas  stream  was  kept  sealed.  A  period  of  at  least  15 
minutes  after  the  syringe  was  reconnected  and  the  gas  flow 
restarted  was  allowed  for  reactor  conditions  to  restabilize 
before  sampling  was  resumed. 

When  either  the  syringe  was  refilled  (i.e.  a  new  liquid 
feed  batch)  or  a  new  catalyst  bed  was  installed,  the 
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subsequent  gas  samples  were  considered  to  constitute  a  new 
series.  Approximately  halfway  through  each  series,  the  fol¬ 
lowing  data  on  the  operating  conditions  were  taken:  atmos¬ 
pheric  pressure,  reactor  gauge  pressure,  hydrogen  and  nitro¬ 
gen  rotameter  readings,  syringe  feeder  setting,  and  fluid¬ 
ized  bed  temperature  set  point. 

When  a  new  catalyst  bed  had  been  installed  and  the 
reactor  was  operating,  the  reactor  temperature  was  monitored 
for  a  period  of  about  15  minutes  using  the  thermocouple 
located  in  the  catalyst  bed  itself.  Since  no  significant 
differences  were  ever  found  between  these  readings  and  those 
from  the  fluidized  bed  temperature  controller,  no  record  was 
kept  of  these  data.  At  this  time,  the  temperatures  of  the 
hot  box,  gas  chromatograph  detector,  sampling  valve  and 
injection  ports  were  also  checked. 

Solutions  of  benzothiophene  in  n-hexane  were  made  up  as 
required.  Typically,  a  100  g  bottle  of  benzothiophene  was 
first  heated  in  a  steam  bath  to  liquify  the  contents.  (The 
melting  point  of  benzothiophene  is  32°C75.)  This  was  added 
to  a  500  ml  flask  containing  200  ml  of  n-hexane  plus  any 
solution  left  over  from  the  previous  batch.  Another  100  ml 
of  n-hexane  was  used  to  dissolve  any  residual  benzothiophene 
remaining  in  the  original  bottle;  this  was  also  added  to  the 
500  ml  flask.  The  flask  was  then  sealed  with  a  ground  glass 
stopper,  shaken  to  mix  the  contents  thoroughly  and  dissolve 
any  benzothiophene  which  might  have  crystallized  on  the 
sides  or  in  the  neck,  and  left  to  cool  to  room  temperature 
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in  a  fume  hood.  These  solutions  were  intended  to  have  an 
approximate  benzothiophene/n-hexane  weight  ratio  of  0.4. 

This  value  was  chosen  so  that  the  solution  would  be  concen¬ 
trated  enough  to  minimize  refilling  of  the  syringe  during 
reactor  operation,  but  still  be  dilute  enough  that  the  reac¬ 
tion  products  would  not  solidify  in  the  condenser. 

The  high  volatility  of  n-hexane  resulted  in  slight  dif¬ 
ferences  in  concentration  of  the  contents  of  the  syringe 
between  one  series  of  gas  samples  and  the  next  series  after 
it  had  been  refilled,  due  to  such  factors  as  evaporation  in 
the  storage  bottle  subsequent  to  the  previous  filling  of  the 
syringe,  evaporation  during  the  transfer  of  solution  from 
bottle  to  syringe,  and  differences  in  the  original  composi¬ 
tion  of  the  solution  as  it  was  made  up.  Cooling  the  solution 
was  impractical,  since  benzothi ophene  would  start  to  crys¬ 
tallize  out  of  solution  at  temperatures  only  slightly  below 
room  temperature.  It  was  impossible  to  sample  the  contents 
of  a  given  syringe  directly  in  many  cases  (such  as  when  the 
syringe  was  refilled  in  the  middle  of  a  run),  so  instead  it 
was  decided  to  calculate  the  composition  of  the  feed  solu¬ 
tion  from  a  molar  balance  based  on  an  analysis  of  the  prod¬ 
uct  stream.  The  details  of  this  calculation  are  given  in  the 
next  chapter. 

When  the  reactor  was  not  in  operation,  all  equipment 
was  kept  hot  (above  250°)  and  a  small  flow  of  nitrogen  was 
kept  running  through  all  lines. 
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The  products  of  typical  runs  were  condensed  at  room 
temperature  and  compared  to  a  sample  of  "condensed  feed", 
produced  by  injecting  the  liquid  feed  into  the  flowing  gas 
stream  and  recondensing  it  after  it  had  bypassed  the  reac¬ 
tor.  These  products  were  qualitatively  analyzed  by  a 
combination  of  gas  chromatography  and  mass  spectrometry. 


5.  EXPERIMENTAL 


A  summary  of  the  run  conditions  is  given  in  Table  2. 
Where  the  series  numbers  are  not  consecutive,  this  indicates 
that  the  accumulated  data  for  a  single  run  ("Data  Point")  is 
a  composite  of  data  from  two  separate  catalyst  beds. 

All  experiments  were  conducted  with  hydrogen  present  in 
great  excess.  The  hydrogen  and  benzothiophene  partial  pres¬ 
sures  in  the  feed  to  the  reactor  were  held  fairly  constant 
at  about  97  and  1.7  kPa,  respectively.  The  temperatures  used 
were  350,  400,  and  450°C,  and  7,  5,  and  4  data  points  were 
obtained  at  each  of  these  respective  temperatures.  At  each 
temperature,  runs  were  performed  with  the  space  time  ranging 
from  1x104  to  7x104  g  catalyst-s/mol  benzothiophene. 

To  verify  the  absence  of  internal  mass  transfer  limita¬ 
tions,  one  run  (Data  Point  6)  was  performed  using  a  signifi¬ 
cantly  smaller  catalyst  pellet  size  than  was  used  in  the 
other  runs  (-60+70  mesh,  as  opposed  to  -40+50).  It  was 
assumed  that  if  the  absence  of  internal  diffusion  effects 
could  be  shown ,  it  would  also  eliminate  the  possibility  of 
external  mass  transfer  limitations  on  the  reaction. 

Although  the  thermocouple  inserted  in  the  catalyst  bed 
gave  some  evidence  against  any  significant  thermal  gradients 
in  the  bulk  phase,  it  could  not  rule  out  completely  the  pos¬ 
sibility  of  internal  or  external  heat  transfer  limitations 
on  the  reaction  rate.  This  reaction  is  highly  exothermic, 
and  if  the  entire  heat  load  from  the  reaction  were  absorbed 
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-  Summary  of  Run  Conditions 


T 

(°C) 

W/F 
x10-  4 

p  ( BT ) 
kPa 

p  ( H  2  ) 
kPa 

Series 

400 

5.912 

1  .74 

102.60 

62,63 

400 

4.750 

1  .28 

97.32 

73,74 

400 

1  .  244 

1  .77 

97.24 

75,76,96 

400 

3.031 

1  .65 

96.97 

90,91 

400 

2.947 

1.87 

91.10 

92,93 

400 

4.011 

1  .73 

96.10 

94,95 

400 

1  .947 

1  .60 

101.35 

68,69,70,97 

450 

5.991 

1  .65 

96.15 

80 

450 

2.057 

1  .34 

96.73 

83,84,85 

450 

1  .227 

1 .79 

97.55 

86,87,88,89 

450 

4.057 

1  .63 

96.29 

79,98 

450 

3.082 

1  .53 

96.71 

81,82,99 

350 

5.500 

1  .68 

95.08 

113,114 

350 

3.001 

1  .54 

97.94 

102,103,104 

350 

3.636 

1  .90 

96.41 

115,116 

350 

6.C31 

1 .44 

97.63 

108 
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by  the  gas  stream,  it  is  estimated  that  the  rise  in  tempera¬ 
ture  would  be  approximately  50°C  (cf.  Appendix  F). 

To  test  this,  one  run  (Data  Point  5)  was  made  with 
about  9%  of  the  hydrogen  in  the  feed  replaced  with  nitrogen. 
If  the  mass  heat  capacities  of  the  two  gases  were  similar, 
this  would  increase  the  specific  heat  of  the  mixture  by 
about  25%.  However,  calculations  made  subsequent  to  the 
actual  experimentation  showed  that  these  mass  heat  capaci¬ 
ties  differed  greatly,  and  that  the  effect  of  substituting 
hydrogen  with  nitrogen  on  heat  capacity  was  negligible. 

Another  data  point  (Point  4)  was  run  using  a  more 
dilute  liquid  feed  solution  than  in  the  other  runs  (0.34 
rather  than  0.43  wt .  ratio  benzothiophene  to  n-hexane);  this 
was  estimated  to  have  increased  the  specific  heat  of  the 
mixture  by  about  6%.  Two  more  data  points  (Points  13  and  14) 
used  a  more  concentrated  solution  (0.51)  than  normal;  this 
was  estimated  to  have  decreased  the  specific  heat  of  the 
mixture  by  about  9%. 
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6.  DATA  ANALYSIS 


All  calculations  described  in  the  following  sections 
were  done  using  the  computer  programs  listed  in  Appendix  5, 
unless  noted  otherwise.  A  sample  output  from  the  main  data 
analysis  program  is  shown  in  Table  3.  This  series  was  chosen 
because  it  represents  both  an  entire  data  point  (Point  5) 
and  a  single  complete  liquid  feed  batch.  The  chromatogram 
shown  in  Figure  4  represents  Run  18  in  this  series. 


6.1  Calculation  of  Reactor  Operating  Conditions  and  Feed 
Rates 

Reactor  temperature  for  any  data  point  was  considered 
identical  to  the  set  point  of  the  fluidized  bed  temperature 
controller.  Total  pressure  in  the  reactor  was  calculated  as 
the  sum  of  atmospheric  pressure  and  the  gauge  pressure  meas¬ 
ured  by  the  manometer  upstream  of  the  syringe  pump.  (It  was 
assumed  that  the  pressure  losses  caused  by  the  1/4”  tubing 
and  the  glass  beads  above  the  catalyst  bed  and  the  bed 
itself  were  negligible.) 

Volumetric  flow  rates  of  hydrogen  and  nitrogen  were 
computed  by  fitting  third-  and  fourth-order  polynomials 
respectively  to  the  calibration  data  using  the  method  of 
least  squares  (cf.  Appendix  1).  Molar  flow  rates  were  calcu¬ 
lated  assuming  ideal  gas  behavior. 

Calculation  of  the  flow  rates  of  n-hexane  and 
benzothiophene  was  more  complex.  Calibration  data  for  the 


62 


v 


63 


Table  3  -  Sample  Output  from  Data  Analysis  Program 


SERIES  92 
800123 


TOTAL  PRESSURE  IN  REACTOR:  104.59  KPA 
TEMPERATURE  IN  REACTOR:  673.2  DEG.  K 
CATALYST:  NALCO  NM  502,  0.1002  G.,  -40+50  MESH 

CATALYST  BED  NUMBER:  23 
LIQUID  FEED  BATCH  NUMBER:  36 
FEED  RATES  TO  REACTOR: 

H2 :  0. 1672E-03  GMOL/SEC 

BT :  0 . 3420E-05  GMOL/SEC 

HEX:  0 . 6826E-05  GMOL/SEC 

N2 :  0. 1340E-04  GMOL/SEC 

PARTIAL  PRESSURES  OF  FEED  TO  REACTOR: 

H2 :  91.63  KPA 

BT:  1.87  KPA 

HEX:  3.74  KPA 

N2 :  7.34  KPA 


0 

.  2930E  +  05 

G-SEC/GMOL ( BT ) 

DATA  (AREA  %) 
RUN  N2+H2 

N-HEXANE 

ETHYL¬ 

STYRENE 

BENZO- 

DIHYDROBENZO 

1 

+H2S 

39.7765 

43.3200 

BENZENE 

2.8470 

0.0193 

THIOPHENE 

14.0100 

THIOPHENE 

0.0 

2 

40.0670 

42.4100 

2.5060 

0.0 

14.8300 

0.0 

3 

40.6509 

42.0400 

2.6660 

0.0 

14.6500 

0.0 

4 

43.7300 

38.0400 

2.7270 

0.0 

15.3000 

0.1859 

5 

42.7600 

40.1900 

2.5360 

0.0 

14.4600 

0.0 

6 

43.1632 

39.3900 

2.6610 

0.0534 

14.5200 

0.0 

7 

41.5336 

41.7900 

2.5520 

0.0 

14.0100 

0.0 

8 

41.6775 

40.9200 

2.4150 

0.0 

14.8800 

0.0361 

9 

42.3756 

39.8500 

2.4650 

0.0 

15.1500 

0.1322 

10 

39.8329 

43.9200 

2.3020 

0.0 

13.9100 

0.0344 

1 1 

40.3800 

42.7600 

2.4060 

0.1629 

14.1500 

0.0 

12 

44.0360 

38.5700 

2.6010 

0.0 

14.7600 

0.0 

13 

40.2400 

42.4600 

2.5820 

0.0305 

14.6900 

0.0 

14 

37.9236 

44.6100 

2.3330 

0.0655 

15.0600 

0.0 

15 

38.3500 

44.6800 

2.8410 

0.0860 

14.0500 

0.0 

16 

40.5959 

40.9100 

2.5580 

0.0 

15.6700 

0.0 

17 

38.7200 

43.8700 

2.3510 

0.0 

14.9000 

0.0647 

18 

40.8600 

40.6600 

2.4060 

0.3094 

15.5600 

0.2132 

19 

37.8500 

44.3700 

2.3240 

0.1304 

15.2600 

0.0 

20 

41.2700 

40.6100 

2.5790 

0.0 

15.2300 

0.1217 

21 

37.9300 

44.5800 

2.2850 

0.0 

14.8400 

0.0672 

22 

42.7300 

39.4100 

2.3960 

0.0 

15.2900 

0.1004 

23 

43.0700 

39.6100 

2.3040 

0.0572 

14.9200 

0.0 

24 

42.4700 

39.2200 

2.4740 

0.0 

15.8400 

0.0 

25 

38.9100 

43.0900 

2.3280 

0.0 

15.3900 

0.0 

26 

42.9600 

38.5500 

2.6930 

0.0 

15.6900 

0.0 

27 

38.9500 

42.8900 

2.6850 

0.0 

15.3800 

0.0641 

28 

40.9800 

41.0600 

2.4860 

0.0 

15.4300 

0.0476 

29 

41.2400 

41.0000 

2.1170 

0.0 

15.5300 

0.0 

30 

45.1200 

37.3800 

2.1170 

0.0 

15. 1700 

0.0 

31 

39.3600 

42.3400 

2.7230 

0.0 

15.4700 

0.0 

32 

39.3652 

43.1600 

2.1700 

0.0 

15.1200 

0.0 

33 

41.3000 

41.9200 

2.2900 

0.0 

14.4000 

0.0610 

34 

40.8700 

41.8100 

2.1260 

0.0 

15.1200 

0.0452 

35 

40.8100 

42.0900 

2.1830 

0.1189 

14.5900 

0.0522 

36 

40.9500 

41.6700 

2.1960 

0.1023 

14.9700 

0.0900 
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37 

38 

39 

42.4600 

39.8200 

44.0400 

40.4000 

42.4800 

37.7200 

2.2460 

2.1010 

2.1480 

0.0 

0.0314 

0.0 

14.8700 

15.3800 

16.0800 

0.0 

0.0 

0.0 

LIQUID 

FEED  COMPOSITION  BASED 

ON  MASS 

BALANCE 

(WT.  RATIO): 

RUN  BT/HEX 

1  0.415425 

2  0.433518 

3  0.437722 

4  0.508265 

5  0.448947 

6  0.465315 

7  0.421179 

8  0.448654 

9  0.471679 

10  0.392046 

11  0.415453 

12  0.477912 

13  0.432769 

14  0.414391 

15  0.405397 

16  0.472103 

17  0.417778 

18  0.482340 

19  0.422835 

20  0.468145 

21  0.407936 

22  0.477575 

23  0.461886 

24  0.494204 

25  0.434830 

26  0.505943 

27  0.448601 

28  0.463342 

29  0.454029 

30  0.488106 

31  0.456208 

32  0.423094 

33  0.423040 

34  0.436525 

35  0.426088 

36  0.440563 

37  0.448008 

38  0.435034 

39  0.509521 

FRACTIONAL  CONVERSION  OF  BENZOTHIOPHENE : 


RUN 

TIME 

TO  EB: 

TO  STY: 

TO  DHBT : 

EB+STY: 

1 

1.667 

0.198996 

0.001340 

0.0 

0.200336 

2 

2.000 

0. 171453 

0.0 

0.0 

0.171453 

3 

2.333 

0. 182237 

0.0 

0.0 

0.182237 

4 

2.917 

0.177416 

0.0 

0.009731 

0.177416 

5 

3.583 

0.176797 

0.0 

0.0 

0.176797 

6 

3.917 

0. 182621 

0.003640 

0.0 

0.186261 

7 

4.500 

0. 182381 

0.0 

0.0 

0. 182381 

8 

4.833 

0.165466 

0.0 

0.001990 

0.165466 

9 

5.250 

0. 164961 

0.0 

0.007118 

0.164961 

10 

5.583 

0.168168 

0.0 

0.002022 

0.168168 

11 

5.917 

0.170363 

0.011458 

0.0 

0.181821 

12 

6.333 

0.177493 

0.0 

0.0 

0.177493 
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13 

10.333 

0.176749 

0.002074 

0.0 

0.178823 

14 

10.583 

0. 158749 

0.004427 

0.0 

0.163176 

15 

10.917 

0.197295 

0.005933 

0.0 

0.203228 

16 

11.167 

0.166599 

0.0 

0.0 

0.166599 

17 

11.500 

0.161353 

0.0 

0.003573 

0.161353 

18 

11.833 

0.154317 

0.019713 

0.011003 

0.174030 

19 

12.167 

0.155817 

0.008685 

0.0 

0. 164502 

20 

12.583 

0.170638 

0.0 

0.006479 

0.170638 

21 

12.917 

0.158049 

0.0 

0.003740 

0.158049 

22 

13.250 

0.160131 

0.0 

0.005399 

0. 160131 

23 

13.667 

0.158409 

0.003907 

0.0 

0.162316 

24 

14.167 

0.160555 

0.0 

0.0 

0.160555 

25 

14.500 

0.156288 

0.0 

0.0 

0.156288 

26 

14.750 

0.173680 

0.0 

0.0 

0.173680 

27 

15.083 

0.175536 

0.0 

0.003372 

0.175536 

28 

15.417 

0.164369 

0.0 

0.002532 

0.164369 

29 

16.833 

0.143051 

0.0 

0.0 

0. 143051 

30 

17.167 

0.145950 

0.0 

0.0 

0.145950 

31 

17.500 

0.177326 

0.0 

0.0 

0.177326 

32 

17.833 

0.149479 

0.0 

0.0 

0.149479 

33 

18. 167 

0.162432 

0.0 

0.003481 

0. 162432 

34 

18.417 

0.146525 

0.0 

0.002507 

0.146525 

35 

18.750 

0.153113 

0.008284 

0.002946 

0. 161398 

36 

19.083 

0.150466 

0.006963 

0.004962  ' 

0.157429 

37 

19.500 

0.156092 

0.0 

0.0 

0.156092 

38 

19.750 

0.143007 

0.002123 

0.0 

0. 145130 

39 

20.083 

0.140585 

0.0 

0.0 

0.140585 

MEAN: 

0.164997 

0.002014 

0.001817 

0. 16701  1 

♦/- 

0.004562 

0.001336 

0.000928 

0.004635 

LIQUID  FEED  BATCH:  36 

COMPOSITE  MASS  BALANCE  BASED  ON  39  DATA  POINTS 
AVG.  FEED  RATIO  BT/HEX:  0.4484  +/-  0.0098 
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syringe  pump  yielded  a  linear  relationship  between  volumet¬ 
ric  flow  rate  and  flow  setting  (cf.  Appendix  5).  It  had  been 
verified  that  the  volumes  of  solid  benzothiophene  and  liquid 
n-hexane  were  additive  in  forming  solutions  in  the  range  of 
concentrations  of  interest  (cf.  Appendix  4),  and  so,  given 
the  relative  concentrations  of  benzothiophene  and  n-hexane 
in  the  solution,  the  molar  flow  rate  could  be  readily  calcu¬ 
lated  . 

The  concentration  of  the  liquid  feed  solution  was  esti¬ 
mated  from  an  analysis  of  the  product  stream.  Since  one  mole 
of  each  of  the  reaction  products,  ethylbenzene,  styrene,  and 
dihydrobenzothi ophene ,  represents  one  mole  of  benzothiophene 
in  the  original  feed,  and  since  these  three  products  com¬ 
prised  an  estimated  99.5%  of  the  total  reaction  products, 
the  mole  ratio  (and  thus  weight  ratio)  of  "equivalent  benzo¬ 
thiophene"  to  n-hexane  could  be  calculated  for  each  gas  sam¬ 
ple.  These  weight  ratios  were  averaged  to  give  a  value  for 
each  liquid  feed  batch.  The  values  of  liquid  feed  concentra¬ 
tion  obtained  from  this  method  were  precise  to  within  about 
3%  based  on  95%  confidence  limits;  this  precision  varied 
from  about  1%  to  5%  depending  mainly  on  the  number  of  mass 
balances  averaged  to  obtain  this  concentration. 

The  space  time  of  benzothiophene  for  each  series  of  gas 
samples  was  calculated  as  the  quotient  of  the  catalyst  bed 
mass  and  the  molar  feed  rate  of  benzothiophene  computed  from 
the  mean  batch  composition  for  that  series.  The  space  time 
for  each  data  point  was  calculated  as  a  weighted  mean  of  the 
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space  times  of  each  series  of  samples  used  to  generate  that 
point : 


where : 


rp 

=  E  (ns 

X  TS)  /  2  (ns 

►  • 

TS 

=  space 

time  for  a  gi 

rp 

=  space 

time  for  a  gi 

ns 

=  total 

number  of  gas 

The 

partial 

pressure  of 

(1) 


lated  for  each  series  as  the  product  of  the  mole  fraction  of 
that  component  in  the  feed  and  the  total  pressure. 

Because  space  time  for  a  given  data  point  was  calcu¬ 
lated  from  the  molar  feed  rate  of  benzothiophene  to  the 
reactor,  which  was  in  turn  based  on  an  analysis  of  the  prod¬ 
uct  stream,  a  relationship  is  introduced  between  the  meas¬ 
urement  errors  of  the  independent  variable,  space  time,  and 
the  dependent  variable,  rate  of  conversion  of  benzothio¬ 
phene.  Consequently,  any  subsequent  statistical  analyses 
which  were  based  on  this  assumption  (such  as  nonlinear 
regression)  are  technically  invalid.  The  relationship 
between  the  errors  of  these  two  variables  appeared  to  be 
fairly  weak,  however,  and  for  this  reason  it  was  assumed 
that  no  serious  errors  should  result  from  this  slight 
relaxation  of  statistical  rigor. 
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6.2  Calculation  of  Fractional  Conversions  of  Reactants  to 
Products 

For  each  gas  sample,  the  fractional  conversions  of 
benzothiophene  to  ethylbenzene,  styrene,  and  dihydrobenzo- 
thiophene  were  calculated  as  the  quotients  of  the  molar  con¬ 
centrations  of  each  product  in  the  sample  and  the  total  of 
the  molar  concentrations  of  these  four  components  in  the 
sample . 

It  was  discovered  early  in  the  course  of  this  study 
that  the  fractional  conversion  of  benzothiophene  to  ethyl¬ 
benzene  decreased  noticeably  with  time.  Typically,  the  cata¬ 
lyst  activity  would  decay  quickly  for  the  first  0  to  3 
hours,  after  which  time  the  rate  of  loss  of  activity  would 
decrease,  reaching  a  constant  value  after  about  5  to  10 
hours  and  continuing  for  as  long  as  it  was  monitored  (40 
hours  in  the  longest  case). 

Kilanowski  and  Gates32  had  independently  found  this 
same  phenomenon,  and  had  eliminated  it  by  adding  small 
quantities  of  H2S  to  their  feed.  Rather  than  introduce  any 
additional  variables  into  this  study,  (i.e.  the  competitive 
adsorption  between  H2S  and  benzothiophene),  it  was  decided 
to  eliminate  the  effects  of  a  decay  in  catalyst  activity  by 
applying  a  transformation  to  the  data.  To  obtain  a  consis¬ 
tent  value  for  catalyst  activity,  a  hypothetical  "zero-time 
activity"  was  defined,  representing  that  activity  which  the 
catalyst  would  have  had  if  it  were  in  equilibrium  with  its 
surroundings  at  the  moment  it  went  on  stream,  before  the 
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long  term  linear  decay  had  had  any  effect.  To  obtain  an 
estimate  of  this,  a  function  of  the  following  form  was 
fitted76  to  the  data: 

F  =  C,  +  C  2  * t  +  C  3  *exp(C  4  *t )  (2) 

where : 


F=  Fractional  conversion  to  ethylbenzene 
t=  Time 


,c 

2  f  C  3  ,  C  4 

=  Parameters  to 

be  fitted 

by 

the  method  o 

least 

squares 

llowing  constraints  were 

placed  on 

the 

parameters : 

0.0 

A 

o 

A 

1.0 

( 2a ) 

1  .0 

<  c2  < 

0.0 

(2b) 

o 

• 

o 

<  c3  < 

1.0 

( 2c ) 

1  .0 

A 

O 

*» 

A 

-0.3 

(2d) 

The  linear  portion  of  this  function  was  extrapolated  back  to 
the  ordinate,  and  this  initial  value  of  fractional  conver¬ 
sion  (i.e.  C, )  was  the  value  subsequently  used.  Two  further 
advantages  of  this  treatment  of  the  data  were  that  it 
enabled  all  of  the  data  from  a  given  catalyst  bed  to  be 
used,  and  it  avoided  the  necessity  of  making  any  arbitrary 
decisions  concerning  the  length  of  time  required  for  a  cata¬ 
lyst  bed  to  reach  steady  state. 
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While  the  constraints  on  the  parameters  may  seem  consi¬ 
derable,  their  removal  had  virtually  no  effect  on  the  resul¬ 
tant  values  of  C,  in  any  of  the  cases  examined.  Further, 
there  was  in  most  cases  some  physical  justification  for 
their  use. 

The  reasons  for  the  upper  and  lower  limits  on  C,  are 
self-evident:  fractional  conversion  must  by  definition  fall 
between  zero  and  unity. 

The  upper  limit  on  C2  precluded  a  steady  state  increase 
in  activity;  while  this  boundary  was  encountered  in  one 
instance,  there  was  no  statistical  justification  for  its 
elimination.  The  lower  limit  on  C2  was  arbitrary;  in  any 
case  it  was  never  encountered. 

Forcing  C3  to  be  positive  precluded  the  possibility  of 
catalyst  activity  increasing  initially  and  later  decreasing. 
Although  this  physical  situation  was  never  encountered,  the 
mathematical  constraint  frequently  was.  If  the  decay  in  cat¬ 
alyst  activity  was  linear  by  the  time  that  sampling  began 
(as  it  occasionally  was)  and  if  the  residual  error  of  the 
first  sample  was  negative,  then  the  curve  fitting  program 
would  attempt  to  force  a  negative  C3  (and  a  positive  C4)  in 
order  to  better  accomodate  this  point.  The  lower  boundary  on 
C 3  was  designed  to  counteract  this.  The  upper  boundary  on  C3 
theoretically  should  be  (1-C, )  and  not  1,  since  (C,+C3) 
represents  the  extrapolation  of  the  experimental  data  to 
time  zero,  and  it  is  this  value  which  must  be  kept  below 
unity.  Neither  value  for  the  upper  limit  of  C3  was 
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approached  in  this  study,  however. 

The  lower  limit  on  C4  was  again  chosen  to  counteract 
the  problem  of  the  residual  of  the  first  gas  sample  affect¬ 
ing  the  shape  of  the  exponential  portion  of  the  fitted 
curve.  A  value  of  1.0  seemed  to  keep  the  problem  reasonably 
well  under  control  while  still  leaving  enough  flexibility  to 
allow  a  good  fit  for  the  cases  in  which  there  actually  was  a 
large  initial  rate  of  decay.  The  upper  limit  on  C4  was  cho¬ 
sen  to  force  the  fitted  curve  to  become  linear  after  a  maxi¬ 
mum  of  about  15  to  20  hours. 

In  Table  3,  it  can  be  seen  that  there  is  a  very  large 
amount  of  scatter  in  the  measurements  of  the  concentrations 
of  styrene  and  dihydrobenzothiophene ,  and  that  the  amounts 
of  these  compounds  in  the  gas  sample  are  very  close  to  the 
minimum  that  could  be  detected  by  the  gas  chromatograph. 
Because  of  this,  no  conclusions  could  be  drawn  concerning  a 
relationship  between  the  fractional  conversions  to  these 
compounds  and  the  time  that  the  catalyst  had  been  on  stream, 
and  it  was  decided  instead  to  calculate  the  fractional  con¬ 
version  for  a  given  data  point  as  the  mean  of  the  values  for 
each  individual  gas  sample.  It  should  be  noted  that  the 
errors  in  the  fractional  conversions  to  these  two  compounds 
is  probably  very  large  -  both  the  fact  that  catalyst  deacti¬ 
vation  was  not  taken  into  account  and  the  fact  that  these 
data  were  taken  so  close  to  the  limits  of  detection  of  the 
apparatus  would  tend  to  yield  values  of  conversion  lower 
than  the  actual  ones. 
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7.  RESULTS  AND  DISCUSSION 


7 . 1  General 

The  three  main  products  from  this  reaction  (if  H2S  is 
not  considered)  were  ethylbenzene,  styrene,  and  dihydro- 
benzothi ophene .  Ethylbenzene  was  the  dominant  product,  and 
comprised  over  95%  of  the  organic  carbon  in  the  products  in 
all  but  two  of  the  cases  studied.  Thus,  many  features  of  the 
overall  reaction  can  be  approximated  by  a  study  of  this  one 
component . 

The  time-dependence  of  the  fractional  conversion  of 
benzothiophene  to  ethylbenzene  for  each  run  is  shown  in  Fig¬ 
ures  5a  through  5p.  Where  the  use  of  a  second  catalyst  bed 
was  necessary  to  supplement  the  data  for  a  particular  set  of 
reaction  conditions,  the  data  from  this  second  bed  have  been 
shown  as  +'s  rather  than  x’s.  The  constants  of  the  exponen¬ 
tial  curves  which  were  fitted  to  these  data  are  summarized 
in  Table  4.  In  all  cases,  the  curves  appeared  to  fit  the 
data  extremely  well. 

The  shapes  of  these  curves  were  similar  for  all  runs: 
either  the  rate  of  loss  of  activity  was  rapid  at  first, 
decreasing  asymptotically  to  a  constant,  or  else  the  rate  of 
deactivation  was  essentially  constant  from  the  time  that 
sampling  was  started.  This  is  in  agreement  with  the  results 
of  most  desulfurization  studies  reported  in  the  literature. 

The  only  exception  to  this  behavior  was  for  data  point 
12.  The  low  conversions  obtained  during  the  first  5  hours  of 
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Data  Point  1 
T=  400  *C 

U/F=  5.912x10*  g-s/mol 
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Figure  5a  -  Time-dependence  of  Fractional  Conversion  to  Ethylbenzene 


Data  Point  2 
T=  400  *C 

U/F=  4. 7  50x10*  g-s/mol 
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Figure  5b  -  Time-dependence  of  Fractional  Conversion  to  Ethylbenzene 


Data  Point  3 
T=  400  °C 

W/F=  1.244x10+  g-s/mol 
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Figure  5c  -  Time-dependence  of  Fractional  Conversion  to  Ethylbenzene 


Data  Point  4 
T=  400  *C 

W/F=  3.031x1Q+  g-s/mol 
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Figure  5d  -  Time-dependence  of  Fractional  Conversion  to  Ethylbenzene 


Data  Point  5 
T=  400  *C 

W/F=  2.947x10+  g-s/moi 
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Figure  5e  -  Time-dependence  of  Fractional  Conversion  to  Ethylbenzene 


* 

Data  Point  6 
T=  400  *C 

U/F=  4.01  lxlO4  g-s/moi 
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Figure  5f  -  Time-dependence  of  Fractionai  Conversion  to  Ethyibenzene 


Data  Point  7 
T=  400  °C 

U/F=  1.947x104  g-s/mol 
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Figure  5g  -  Time-dependence  of  Fractional  Conversion  to  Ethylbenzene 


Data  Point  8 
T=  450  *C 

U/F=  5 . 99  1  x10+  g-s/mol 
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Figure  5h  -  Time-dependence  of  Fractional  Conversion  to  Ethylbenzene 


Data  Point  9 
T=  450  *C 

U/F=  2 .  057  x10+  g-s/mol 
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Figure  51  -  Time-dependence  of  Fractional  Conversion  to  Ethylbenzene 


Data  Point  10 
T=  450  *C 

U/F=  1  .227x10+  g-s/mol 
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SAMPL I  NO  TIME,  HRS 

Figure  5j  -  Time-dependence  of  Fractional  Conversion  to  Ethylbenzene 


Data  Point  1 1 
T=  450  *C 

W/F=  4.Q57x1Q+  g-s/mol 
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Figure  5k  -  Time-dependence  of  Fractional  Conversion  to  Ethylbenzene 


Data  Point  1 2 
T=  450  °C 

U/F=  3  .  Q82x1Q+  g-s/moi 
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Figure  51  -  Time-dependence  of  Fractional  Conversion  to  Ethylbenzene 
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Data  Point  13 
T=  350  *C 

W/F=  5.500x104  g-s/mol 
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SRMPLING  TIME.  HRS 

Figure  5m  -  Time-dependence  of  Fractional  Conversion  to  Ethyibenzene 


Data  Point  1 4 
T=  350  T 

U/F=  3.Q01x1Q+  g-s/mol 
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Figure  5n  -  Time-dependence  of  Fractional  Conversion  to  Ethylbenzene 


Data  Point  1 5 
T=  350  X 

W/F=  3.636x104  g-s/mol 
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Figure  5o  -  Time-dependence  of  Fractional  Conversion  to  Ethylbenzene 


Data  Point  1 6 
T=  350  *C 

U/F=  6.031x104  g-s/moi 
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Figure  5p  -  Time-dependence  of  Fractional  Conversion  to  Ethylbenzene 
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Table  4  -  Constants  of  Catalyst  Deactivation  Curve 


Data 

Point 

c, 

c2 

x103 

c3 

c4 

c2/c, 

1 

0.3651 

1.9996 

0.0667 

0.0667 

0.5476 

2 

0.2873 

3.0126 

0.0996 

0.0996 

1 .0486 

3 

0.0850 

1 .4343 

0.0330 

0.0330 

1 .6874 

4 

0.1983 

1.5598 

0.0610 

0.0610 

0.7866 

5 

0.1815 

1 .4702 

0.0487 

0.0487 

0.8100 

6 

0.2556 

1 .8228 

0.1033 

0.1033 

0.7131 

7 

0.1155 

0.7339 

0.0911 

0.0911 

0.6354 

8 

0.3753 

2.6228 

0.0377 

0.0377 

0.6989 

9 

0.1244 

2.2231 

0.0392 

0.0392 

1 .7822 

10 

0.0777 

1.4813 

0.1125 

0.1125 

1.9052 

1  1 

0.2546 

2.7376 

0.0075 

0.0075 

1 .0752 

12 

0.1898 

2.7186 

0.0000 

0.0000 

1 .4323 

13 

0.2182 

0.0000 

0.0876 

0.0876 

0.0000 

14 

0.1045 

0.8593 

0.0512 

0.0512 

0.8223 

15 

0.1421 

1.5221 

0.0719 

0.0719 

1.0922 

16 

0.1912 

0.2048 

0.0836 

0.0836 

0.1071 

- 


90 


operation  for  this  point  were  due  to  a  different  catalyst 
bed  having  been  used  in  this  period,  for  which  the  operating 
conditions  may  have  differed  slightly.  Rather  than  discard 
these  data,  it  was  decided  that  a  composite  of  the  data  from 
the  two  beds  would  give  a  better  estimate  of  the  initial 
activity  of  the  catalyst  than  would  either  bed  alone,  and 
for  this  reason,  this  bed  was  treated  identically  to  the 
other  beds  used  to  obtain  the  rest  of  the  data. 

The  conversion  for  Point  5  was  found  to  be  slightly 
low;  this  may  be  explained  by  the  fact  that  the  hydrogen 
partial  pressure  was  a  few  percent  lower  than  those  in  the 
other  runs. 


7.2  Reaction  Kinetics 

The  fractional  conversions  of  benzothiophene  to  ethyl¬ 
benzene,  styrene,  and  dihydrobenzothiophene  are  presented  in 
Table  5  and  plotted  as  a  function  of  space  time  in  Figures 
6-8,  respectively. 

From  Figure  6,  it  can  be  seen  that  over  the  range  of 
reaction  conditions  used,  the  reaction  rate,  obtained  from 
the  slope  of  the  curve,  is  independent  of  the  conversion  of 
benzothiophene  up  to  fractional  conversions  as  high  as  0.4. 
This  may  be  interpreted  as  evidence  for  intrinsic  zero-order 
kinetics,  however  if  the  space  time  is  increased  much  beyond 
the  range  used  here,  the  conversions  do  begin  to  fall  below 
the  linear  relationships  shown  in  this  figure.  This 
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Table  5  -  Fractional  Conversion  Data 


Data 

T 

W/F 

Point 

(°C) 

x10-  4 

1 

400 

5.912 

2 

400 

4.750 

3 

400 

1  .244 

4 

400 

3.031 

5 

400 

7.947 

6 

400 

4.011 

7 

400 

1  .947 

8 

450 

5.991 

9 

450 

2.057 

10 

450 

1  .227 

1  1 

450 

4.057 

12 

450 

3.082 

13 

350 

5.500 

14 

350 

3.001 

15 

350 

3.636 

16 

350 

6.031 

To 

To 

To 

EB 

STY 

DHBT 

.  3670 

0.0018 

0.0005 

.2870 

0.0011 

0.0002 

.  0984 

0.0035 

0.0007 

.  1949 

0.0021 

0.0019 

.  1832 

0.0024 

0.0017 

.2590 

0.0021 

0.0012 

.1219 

0.0027 

0.0005 

.4007 

0.0080 

0.0002 

.  1268 

0.0100 

0.0004 

.1815 

0.0073 

0.0006 

.  270  1 

0.0143 

0.0014 

.  1985 

0.0152 

0.0011 

.2195 

0.0013 

0.0043 

.  1047 

0.0001 

0.0012 

.  1420 

0.0009 

0.0063 

.  1950 

0.0009 

0.0027 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


350  °C 
400  °C 
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Figure  6  -  Effect  of  Space  Time  on  Fractional  Conversion  to  Ethylbenzene 
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Figure  7  -  Effect  of  Space  Time  on  Fractional  Conversion  to  Styrene 


0.020i 


94 


u  o  u 

«  «  o 

o  o  o 

LD  O  LD 
CO 


<]  H  O 


LO 

tH 

o 

0 

o 


NO  I S^BANQD  lUNOIlOd^J 


Figure  8  -  Effect  of  Space  Time  on  Fractional  Conversion  to  Dihydrobenzothiophene 
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indicates  that  the  intrinsic  rate  is  not  zero-order. 

Much  more  scatter  can  be  seen  in  the  data  taken  at 
350°C  than  at  the  other  temperatures.  This  was  probably 
caused  by  a  much  slower  equilibration  of  the  catalyst  with 
its  surroundings  at  this  temperature. 

The  absence  of  effects  from  heat  and  mass  transfer  lim¬ 
itations  can  also  be  seen  in  Figure  6:  neither  a  decrease  in 
catalyst  particle  size  (Point  6)  nor  an  increase  (Point  4) 
or  decrease  (Points  13,  14)  in  the  specific  heat  of  the  gas 
phase  appeared  to  produce  any  significant  effect  on  frac¬ 
tional  conversion. 

The  observed  reaction  rates  at  the  three  temperatures 
studied  are  shown  in  Table  6  and  the  logarithms  of  these 
rates  are  plotted  as  a  function  of  reciprocal  temperature  in. 
Figure  9. 

The  dependence  of  the  reaction  rate  on  temperature  is 
striking:  although  the  overall  rate  increased  significantly 
from  350  to  400°C,  no  increase  was  observed  from  400  to 
450°C,  implying  the  existence  of  a  maximum  or  a  point  of 
inflection  in  this  relationship.  Possible  causes  of  this 
behavior  will  be  discussed  below. 

The  conversions  to  styrene  and  dihydrobenzothiophene 
are  shown  in  Figures  7  and  8,  respectively.  In  each  case, 
the  conversion  appears  to  rise  to  a  maximum  with  increasing 
space  time  and  then  to  decrease,  although  the  data  for  dihy¬ 
drobenzothiophene  are  too  scattered  to  be  conclusive.  This 
would  support  the  conclusion  of  Furimsky  and  Amberg14  that 
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Table  6 


-  Observed  Reaction  Rates 


Temperature 

Rate 

°C 

mol/g  cat-s 

350 

0.03611 

400 

0.06243 

450 

0.06225 
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both  are  intermediates  in  the  overall  reaction. 

The  conversion  to  styrene  increases  with  increasing 
temperature,  while  the  conversion  to  dihydrobenzothiophene 
decreases.  If  these  intermediates  are  seen  as  representing 
respective  hydrogenolys i s-hydrogenat i on  and  hydrogenat i on- 
hydrogenolys i s  paths,  and  if  direct  conversion  from  dihydro¬ 
benzothiophene  to  styrene  is  considered  to  be  small,  then  it 
is  implied  that  the  activation  energies  for  hydrogenolysi s 
are  higher  than  for  hydrogenation  in  both  cases.  This  is 
what  would  be  expected  from  kinetic  theory. 

This  mechanism  is  in  itself  insufficient  to  explain  the 
lack  of  increase  in  the  overall  reaction  rate  between  400 
and  450°C,  even  if  conversion  of  dihydrobenzothiophene  to 
styrene  is  considered.  At  least  two  alternative  hypotheses 
could  reasonably  explain  this  result. 

One  possibility  is  that  the  lack  of  change  in  reaction 
rate  is  a  consequence  of  a  shift  in  the  equilibrium  between 
benzothiophene  and  dihydrobenzothiophene.  Givens  and 
Venuto16,  Furimsky  and  Amberg14  and  Kilanowski  et  al.31  all 
found  that  the  rate  of  dehydrogenation  of  dihydrobenzothio¬ 
phene  was  significant  even  at  the  relatively  high  hydrogen 
partial  pressures  used  in  this  study  (3-4  orders  of  magni¬ 
tude  higher  than  those  of  dihydrobenzothiophene).  Furimsky 
and  Amberg  also  found  a  higher  activation  energy  for  dehy¬ 
drogenation  than  for  desulfurization,  and  from  this  pre¬ 
dicted  that  dehydrogenation  would  increase  in  importance 
with  respect  to  desulfurization  as  a  rate-limiting  step  as 


J 


■ 


99 


the  temperature  increased. 

If  this  reverse  reaction  is  in  fact  rate-limiting,  it 
is  important  to  estimate  its  relevance  for  the  operation  of 
commercial  hydrotreating  reactors.  Since  these  reactors 
operate  at  hydrogen  partial  pressures  approximately  two 
orders  of  magnitude  higher  than  in  the  present  study,  it  is 
expected  that  the  temperature  at  which  the  rate  of  dehydro¬ 
genation  becomes  critical  should  be  higher  than  found  here. 
If  the  forward  reaction  is  first  order  in  hydrogen  (as  has 
been  found  for  the  overall  desulfurization  reaction  by 
Rollman29  and  Bartsch  and  Tanielan15),  and  the  heat  of  reac¬ 
tion  is  assumed  to  be  30  kcal/mole  (typical  for  hydrogena¬ 
tions),  then  an  increase  in  hydrogen  partial  pressure  by  a 
factor  of  100  would  raise  this  "critical  temperature"  by 
about  175°C.  This  is  in  agreement  with  the  prediction  of 
Furimsky66  that  the  dehydrogenation  reaction  should  not  be 
an  important  factor  at  high  pressures.  This  estimate  of 
"critical  temperature"  may  be  misleading,  however,  in  that 
it  does  not  allow  for  any  effects  of  H2S,  which  has  been 
found  to  severely  retard  some  hydrogenations5.  Consequently, 
the  actual  temperature  at  which  this  dehydrogenation  rate 
should  become  important  in  an  industrial  hydrotreater  with 
high  H2S  partial  pressures  may  be  much  lower  (or  even  in  the 
opposite  direction)  than  the  predicted  value. 

A  second  possible  cause  of  the  unusual  character  of  the 
temperature  dependence  of  the  overall  rate  is  that  the  rate 
of  adsorption  of  one  of  the  reactants  may  have  been  a  rate- 
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controlling  step.  If,  for  example,  the  reaction  took  place 
on  two  different  catalyst  sites,  one  more  active  for  initial 
hydrogenolys i s  of  benzothiophene  to  styrene  and  one  more 
active  for  initial  hydrogenation  to  dihydrobenzothiophene , 
and  if  the  exothermic  adsorption  of  either  benzothiophene  or 
hydrogen  onto  the  latter  site  were  stronger  than  onto  the 
former,  then  (at  some  point),  an  increase  in  temperature 
could  also  result  in  no  change  in  the  overall  rate  as  well 
as  in  an  increase  in  styrene  and  a  decrease  in  dihydrobenzo¬ 
thiophene  in  the  reaction  products.  Consequently,  a  reaction 
model  incorporating  adsorption  limitations  on  the  reaction 
path  through  dihydrobenzothiophene  could  also  predict 
kinetic  data  consistent  with  the  results  of  this  study. 

Most  studies  relevant  to  the  evaluation  of  this  hypoth¬ 
esis  have  been  conducted  at  too  low  a  temperature  to  permit 
any  definite  conclusions  to  be  drawn  with  respect  to  the 
present  data.  Lipsch  and  Schuit49  found  a  maximum  in  the 
adsorption  of  hydrogen  onto  supported  cobalt-molybdenum  cat¬ 
alysts  at  about  550°C,  and  Owens  and  Amberg77  found  that  the 
activation  energy  for  chemisorption  of  hydrogen  onto  chromia 
had  the  right  order  of  magnitude  to  have  been  able  to  affect 
the  rate  of  thiophene  desulfurization  at  250°C.  Owens  and 
Amberg  also  found  that  the  adsorption  of  thiophene  did  not 
appear  to  be  a  rate  controlling  step  in  this  reaction  at 
temperatures  up  to  about  370°C,  but  allowed  that  this  may 
not  be  true  at  higher  temperatures. 
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Desikan  and  Amberg78,  working  with  methyl-substituted 
thiophenes  at  temperatures  up  to  350°C,  concluded  that 
adsorption  of  either  thiophene  or  hydrogen  onto  Co-Mo- 
y-Al  2 0 3  could  not  be  the  only  rate-determining  step,  and 
Kawaguchi  et  al. 1 1 ,  studying  thiophene  desulfurization  in 
the  same  temperature  range,  concluded  that  adsorption  of 
hydrogen  onto  Ni-Mo-y-Al 203  was  not  significant  in  determin¬ 
ing  the  overall  reaction  rate.  Frye  and  Mosby7’,  working  at 
high  pressures  and  330°C,  also  suggested  that  hydrogen  chem¬ 
isorption  was  probably  not  significant. 

If  the  reaction  path  is  highly  temperature-dependent, 
then  one  may  ask  on  what  basis  have  previous  researchers 
calculated  apparent  activation  energies  and  pseudo-first 
order  rate  constants  for  this  reaction.  From  the  experi¬ 
mental  data  for  350  and  400°C,  two  points  alone,  a  value  of 
9.1  kcal/mol  was  estimated  as  a  minimum  value  for  the  appar¬ 
ent  activation  energy  of  the  reaction  path  through  dihydro- 
benzothiophene .  From  this,  a  pseudo-zero  order  preexponen¬ 
tial  factor  of  57.33  mol/g  cat-s  was  calculated.  This  value 
of  activation  energy  is  shown  in  Table  7,  along  with  the 
previously  reported  literature  values.  The  details  of  this 
calculation  are  given  in  Appendix  G. 

The  experimentally  obtained  minimum  value  for  the  acti¬ 
vation  energy  falls  within  the  range  of  values  reported  by 
other  researchers,  however  the  error  in  this  value  is  sig¬ 
nificant,  and  it  is  likely  that  the  true  apparent  activation 
energy  for  this  path  may  be  somewhat  higher.  The  use  of  only 
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Table  7  -  Experimental  and  Literature  Values  of 
Activation  Energy 


Reference 

Catalyst 

Ea 

kcal/mol 

Temp.  Range 
°C 

Bartsch  and  Tanielan15 

Co-Mo- y-Al 2  0  3 

4.9 

300-400 

Bartsch  and  Tanielan28 

Fe  203  -y-Al 2 0 3 

7.9-14.0 

300-400 

Furimsky  and  Amberg14 

Co-Mo-y-Al 203 

12.0 

250-400 

Daly  1  7 

Co-Mo-y-Al 2  0  3 

16.  1 

200-400 

Kilanowski  and  Gates32 

Co-Mo-y-Al 20  3 

20±3 

252-332 

Present  Study 

N i -Mo-y -Al 2  0  3 

>9.  1 

350-400 

* 
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two  points  in  the  estimation  of  this  activation  energy  cer¬ 
tainly  introduces  some  error  into  the  calculations,  and  the 
fact  that  the  minimum  temperature  used  in  this  calculation 
was  somewhat  higher  than  those  used  in  previous  studies 
meant  that  the  present  value  is  probably  less  reliable  as  an 
estimate  of  the  true  activation  energy  of  the  low  tempera¬ 
ture  reaction  path  than  are  the  other  reported  values.  Since 
the  apparent  activation  energy  accounts  for  the  kinetics  of 
the  adsorption  and  desorption  of  the  reactants  and  products 
as  well  as  for  the  kinetics  of  the  actual  reaction,  it  is 
possible  that  the  use  of  a  nickel-promoted  catalyst  in  the 
present  study  rather  than  a  cobalt  catalyst  may  also  have 
been  responsible  for  some  of  the  discrepancy  between  the 
observed  minimum  value  of  apparent  activation  energy  and 
previously  reported  values  of  this  parameter. 

7.3  Trace  Reaction  Products 

A  chromatogram  of  the  condensate  from  a  typical  set  of 
reaction  conditions  is  shown  in  Figure  10  along  with  the 
molecular  weights  associated  with  each  peak  by  mass  spec¬ 
trometry.  These  results  are  only  qualitative  (and  at  best 
incomplete),  since  most  of  the  lighter  products  would  have 
been  lost  from  the  condenser  in  the  flowing  gas  stream,  but 
it  is  estimated  (using  benzothiophene  as  a  basis)  that  none 
of  these  trace  products  was  present  in  excess  of  0.1%  of  the 
total  products,  and  that  most  were  present  at  a  level  of 
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Figure  10 


Chromatogram  of  Condensed  Products 
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0.01%  or  lower.  Another  factor  affecting  these  results  is 
the  relatively  large  amounts  of  the  four  main  components  of 
the  product  stream  with  respect  to  the  trace  products:  it  is 
almost  certain  that  a  small  peak  with  a  retention  time  close 
to  that  of  one  of  the  main  products  would  have  been  lost  in 
this  larger  peak,  especially  if  it  had  the  same  molecular 
weight  as  the  large  peak. 

The  only  peaks  present  in  an  analysis  of  the  liquid 
feed  were  ethylbenzene,  styrene,  benzothiophene ,  dihydro- 
benzothi ophene ,  the  peak  with  molecular  weight  146  eluted 
just  prior  to  benzothiophene,  and  the  peaks  with  masses  176 
and  226  eluted  just  after  dihydrobenzothiophene .  It  appeared 
that  the  relative  areas  of  all  these  peaks  in  the  product 
analysis  were  significantly  greater  than  in  the  feed,  and  as 
a  result  it  was  concluded  that  all  were  legitimate  products 
of  the  reaction,  rather  than  just  impurities. 

A  summary  of  the  molecular  weights  of  all  the  compo¬ 
nents  of  the  product  stream,  along  with  possible  identities 
of  some  of  the  trace  products  is  given  in  Table  8.  These 
identities  could  not  be  checked  by  an  analysis  of  retention 
times  because  the  pure  compounds  were  unavailable. 

From  this  table,  it  can  be  inferred  that  a  small  degree 
of  cracking  of  the  reactants  and/or  products  occurs,  along 
with  some  recombination  of  the  hydrocarbon  fragments.  The 
presence  of  some  species  implies  that  their  isomers  should 
also  be  present;  in  some  cases  peaks  corresponding  to 
possible  structural  isomers  were  found,  in  some  cases  they 
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Table  8  -  Composition  of  Product  Stream 
Molecular  Number 

Weight  of  Peaks  Possible  Identity 


78  1 


92 


1 


104 


106 


120 


.SH 


+  isomers 


132 


,CH2  SH 


+  isomers 


134 


107 


Table  8  ( cont . ) 


Molecular  Number 

Weight  of  Peaks  Possible  Identity 


136 


142 


1 


146 


CH2  CH2  SH 


148 


CH2  CH2  SH 


176 


210  1 


226  1 


108 


Table  8  (cont . ) 


Molecular 

Weight 


Number 
of  Peaks 


Possible  Identity 


346 
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were  not,  possibly  due  to  inadequate  separation  or  sensitiv¬ 
ity  of  the  apparatus. 

A  peak  at  molecular  weight  210  was  identified  as  disty¬ 
rene  and  a  second  peak  at  M.W.  346  could  represent  the  prod¬ 
uct  of  a  reaction  between  distyrene  and  benzothiophene .  This 
indicates  that  even  at  400°C  some  degree  of  styrene  polymer¬ 
ization  was  occurring. 

Furimsky  and  Amberg  had  predicted  that  2-mercapto- 
phenylacetylene  (M.W.  134)  and  2-mercaptostyrene  (M.W.  136) 
might  be  found  in  the  reaction  products  in  trace  quantities. 
In  the  present  analysis,  only  one  unidentified  peak  corre¬ 
sponding  to  either  of  these  masses  was  found,  at  M.W.  134. 

Daly17  found  separate  peaks  which  were  attributed  to  1- 
and  2-phenylethanethiol .  The  reaction  conditions  used  in  the 
present  study  (low  pressure  and  low  conversion)  would  not  be 
expected  to  favor  a  back-reaction  between  styrene  and  H2S, 
Daly’s  major  path  for  the  formation  of  these  products,  but 
it  could  allow  for  2-phenylethanethiol  to  be  formed  from 
dihydrobenzothiophene  through  2-phenylethenethiol .  No  evi¬ 
dence  was  found  for  the  presence  of  either  of  these  com¬ 
pounds  in  the  reaction  products,  although  it  is  quite  possi¬ 
ble  that  a  peak  representing  the  intermediate  may  have  been 
lost  in  the  much  larger  dihydrobenzothiophene  peak. 
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7.4  Catalyst  Deactivation 

Only  one  researcher  (Kilanowski  et  al.32)  has  reported 
a  loss  in  catalyst  activity  comparable  to  that  found  in  the 
present  study.  Kilanowski  concluded  that  the  decrease  in 
activity  was  due  to  a  loss  of  sulfur  from  the  catalyst  sur¬ 
face. 

To  test  whether  this  explanation  is  consistent  with  the 
results  of  the  present  study,  or  whether  another  explanation 
such  as  coking  on  the  catalyst  might  be  required,  the  rate 
of  deactivation  was  normalized  for  each  data  point  to 
account  for  differences  in  fractional  conversion  between  the 
runs.  This  normalized  rate  of  decay  (expressed  as  the  frac¬ 
tion  of  initial  activity  lost  per  hour)  was  calculated  as 
the  rate  of  decrease  of  fractional  conversion  divided  by  the 
initial  fractional  conversion  (i.e.  C2/C, ) .  (Although  a  mean 
value  of  fractional  conversion  might  have  been  chosen  rather 
than  the  initial  value,  this  would  not  have  affected  the 
normalized  rate  significantly.) 

The  normalized  deactivation  rate  is  shown  in  Table  4, 
and  is  plotted  against  fractional  conversion  to  ethylbenzene 
in  Figure  11.  (Since  fractional  conversion  to  ethylbenzene 
was  found  to  be  directly  proportional  to  space  time,  the 
shapes  of  these  curves  will  be  similar  regardless  of  which 
of  these  variables  is  used  as  the  abscissa.)  The  curves  in 
this  figure  are  exponentials  which  were  fitted  using  a  least 
squares  algorithm  with  no  constraints  on  either  of  the 
parameters,  and  were  included  only  to  show  the  general 
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relationships  between  the  variables  rather  than  to  suggest 
any  particular  mathematical  dependence  or  to  provide  any 
basis  for  extrapolation  of  these  data. 

Although  the  errors  in  these  deactivation  rates  are 
large,  some  general  conclusions  can  still  be  drawn  from  the 
data . 

If  the  deactivation  rate  were  related  to  a  chemical 
change  on  the  catalyst  surface,  a  deactivation  rate  with  a 
negative  order  with  respect  to  one  of  the  reaction  products 
may  be  more  plausible.  One  hypothesis  which  would  predict 
this  type  of  behaviour  is  that  the  loss  in  activity  is  the 
result  of  a  loss  of  sulfur  from  the  catalyst  surface,  as 
suggested  by  Broderick  et  al.73  and  Kilanowski  et  al.32  In 
this  case,  it  would  be  expected  that  an  increase  in  frac¬ 
tional  conversion  of  benzothiophene  (and  hence  an  increase 
in  H2S  partial  pressure)  would  lead  to  a  decrease  in  deacti¬ 
vation  rate. 

Another  explanation  for  this  observed  long  term  deacti¬ 
vation  might  be  coking  reactions  on  the  catalyst  surface.  If 
the  overall  rate  of  coke  formation  is  zero  order  in  benzo¬ 
thiophene  (as  has  been  found  for  desulfurization),  then  this 
deactivation  rate  will  be  independent  of  space  time.  If,  in 
addition,  the  rate  of  coking  is  a  positive  order  with 
respect  to  any  (or  all)  of  the  products  of  the  desulfuriza¬ 
tion  reaction,  then  an  increase  in  space  time  would  result 
in  an  increase  in  the  partial  pressures  of  these  components 
in  the  reactor,  and  thus  the  deactivation  rate  should 
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increase  with  increasing  space  time.  It  is  only  when  the 
rate  of  coking  has  a  negative  order  with  respect  to  one  of 
the  reaction  products  (perhaps  due  to  blocking  of  one  of  the 
catalyst  sites  by  a  harmless  component)  that  a  negative 
reaction  order  may  be  possible. 

The  nature  of  the  trace  products  of  the  desulfurization 
reaction  would  seem  to  indicate  that  coke  formation  could 
result  from  reactions  between  styrene,  benzothiophene , 
hydrogen  sulfide,  and  thiophenic  fragments  formed  by  crack¬ 
ing  of  the  benzene  ring  in  benzothiophene.  If  this  is  in 
fact  occurring,  then  the  rate  of  coking  should  have  a 
positive  order  with  respect  to  these  compounds.  A  negative 
order  with  respect  to  dihydrobenzothiophene  is  also 
unlikely,  since  it  appears  to  be  in  equilibrium  with  benzo¬ 
thiophene,  and  consequently  an  increase  in  dihydrobenzothio¬ 
phene  partial  pressure  would  also  be  associated  with  an 
increase  in  that  of  benzothiophene.  Since  an  increase  in 
ethylbenzene  pressure  can  only  be  the  result  of  increased 
amounts  of  styrene  and/or  dihydrobenzothiophene,  it  is  also 
unlikely  to  be  an  inhibitor  of  coking. 

Hydrogen  sulfide  was  found  by  Gissy  et  al.74  to  inhibit 
coking  of  n-dodecane  on  Co-Mo-y-Al 203  in  a  similar  reactor 
to  that  used  in  the  present  study.  It  was  suggested  by  Gissy 
that  H2S  protected  the  catalyst  by  improving  its  hydrogena¬ 
tion  function,  and  thus  limiting  the  concentration  of  ole- 
finic  coke  precursors  in  the  reactor.  If  this  hypothesis  is 
applicable  to  the  present  case,  then  the  rate  of 
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deactivation  should  decrease  with  increasing  space  time  - 
the  same  result  as  would  be  predicted  by  the  "loss  of  cata¬ 
lyst  sulfur"  hypothesis. 

The  data  in  Figure  11  support  both  the  hypothesis  that 
deactivation  is  the  result  of  a  loss  of  sulfur  from  the  cat¬ 
alyst  and  the  hypothesis  that  H2S  inhibits  coking  by 
improving  the  hydrogenation  function  of  the  catalyst:  even 
very  small  H2S  partial  pressures  (on  the  order  of  0.5  kPa) 
are  associated  with  a  significant  decrease  in  the  rate  of 
loss  of  activity. 

Other  data  in  the  present  study  would  also  seem  to 
favor  the  former  explanation,  however.  Not  only  were  heavy 
aliphatic  hydrocarbons  such  as  n-dodecane  absent  from  the 
present  study,  but  an  analysis  of  the  reaction  products 
would  seem  to  indicate  a  greater  effect  of  increasing  tem¬ 
perature  on  potential  for  coking  than  was  observed.  It 
appears  that  the  change  in  deactivation  rate  from  350  to 
450°C  is  fairly  small  for  a  fixed  conversion  of  benzothio- 
phene,  increasing  by  a  factor  of  between  two  and  ten,  but 
that  the  ratio  of  the  amounts  of  styrene  (presumably  a  main 
coke  precursor)  produced  by  the  reaction  at  these  two  tem¬ 
peratures  may  be  much  greater,  perhaps  between  20:1  and 
50:1,  based  on  the  relative  amounts  of  styrene  and  dihydro- 
benzoth i ophene  in  the  product  stream.  In  addition,  the  data 
for  all  three  temperatures  show  both  a  low  rate  of  deactiva¬ 
tion  and  a  relatively  high  degree  of  unsaturation  in  the 
products  stream  at  intermediate  to  high  values  of  space 
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time.  Unfortunately,  these  data  are  not  sufficiently 
accurate  to  completely  eliminate  the  possibility  of  coking. 
Regardless  of  which  of  these  two  explanations  is  correct  (if 
either),  if  the  protection  of  the  catalyst  is  in  fact  due  to 
H2S,  this  type  of  deactivation  would  not  be  expected  to  be  a 
significant  factor  in  high-pressure  processes,  and  could 
explain  why  this  effect  has  not  been  found  in  similar 
studies  conducted  at  high  pressure. 

Other  hypotheses  could  also  be  consistent  with  the  cat¬ 
alyst  deactivation  data  of  the  present  study.  For  example, 
the  presence  of  hydrogen  sulfide  may  protect  the  catalyst 
surface  in  a  different  way  (perhaps  by  hindering  migration 
of  metal  ions  below  the  catalyst  surface),  or  that  this 
protection  may  be  mainly  due  to  the  presence  of  a  different 
product,  or  that  deactivation  may  be  the  result  of  several 
concurrent  processes.  If  desulfurization  occurs  on  two  dif¬ 
ferent  sites,  it  is  also  possible  that  one  may  be  affected 
and  not  the  other,  or  that  one  type  of  site  may  be  slowly 
transformed  into  the  other.  These  possibilities  were  not 
investigated  further,  due  to  the  large  amount  of  scatter  in 
the  data  and  the  consequent  difficulty  of  obtaining  meaning¬ 
ful  results  from  subsequent  tests. 
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8.  CONCLUSIONS 


1.  A  fixed-bed  microreactor  with  on-line  gas  sampling 

provides  an  effective  method  for  studying  hydrodesul¬ 
furization  reactions,  and  has  the  advantages  over  an 
analysis  of  condensed  products  that  the  analysis  is 
quantitative  and  that  small  changes  in  catalyst  activ¬ 
ity  with  time  on  stream  can  be  detected  more  easily. 

2.  The  pretreatment  method  of  Ripperger  and  Saum20  is  effec¬ 

tive  in  producing  a  catalyst  with  consistent  activity. 

3.  Benzothiophene  is  desulfurized  via  two  main  reaction 

paths.  At  lower  temperatures,  benzothiophene  is  first 
hydrogenated  to  dihydrobenzothiophene ,  which  is  then 
desulfurized  to  ethylbenzene.  At  higher  temperatures, 
benzothiophene  is  first  desulfurized  to  styrene,  which 
is  then  hydrogenated  to  produce  ethylbenzene.  It  is 
suggested  that  this  could  have  been  due  either  to  a 
shift  in  the  equilibrium  between  benzothiophene  and 
dihydrobenzothiophene  or  to  adsorption  limitations  on 
one  of  the  reactants.  Because  of  this,  a  rate  expres¬ 
sion  which  does  not  allow  for  a  change  in  the  main 
reaction  path  with  a  change  in  temperature  may  be 
applicable  only  over  a  very  limited  range  of  tempera¬ 
tures. 
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4.  A  minimum  value  of  9.1  kcal/mol  was  calculated  for  the 

apparent  activation  energy  of  the  reaction  path  through 
dihydrobenzothiophene .  This  is  within  the  range  of 
other  reported  literature  values,  however  it  is  sug¬ 
gested  that  the  true  value  may  be  somewhat  higher. 
Possible  reasons  for  this  discrepancy  include  the  error 
caused  by  the  fact  that  the  rates  at  only  two  tempera¬ 
tures  were  used  to  calculate  this  quantity,  the  wider 
temperature  range  used  by  other  researchers  which  may 
be  more  representative  of  the  reaction  path,  and  the 
differences  between  the  nickel-molybdenum  catalyst  used 
in  this  study  and  the  cobalt-molybdenum  and  iron  cata¬ 
lysts  used  by  other  researchers. 

5.  Other  reactions  taking  place  to  a  lesser  degree  in  the 

reactor  include  cracking  of  the  benzene  ring,  dealkyla¬ 
tion  of  ethylbenzene,  and  reaction  of  the  resultant 
hydrocarbon  fragments  with  each  other,  hydrogen 
sulfide,  and  the  other  reactants  and  products.  It  is 
estimated  that  these  reactions  involve  only  about  0.5% 
of  the  total  reaction  products.  In  addition  to  this, 
evidence  was  found  for  a  slight  degree  of  styrene 
polymerization  in  the  reactor. 

6.  Catalyst  activity  was  found  to  decay  linearly  at  a  rate 

on  the  order  of  0.1%/hr.  Although  this  necessitated 
defining  a  hypothetical  "zero-time"  activity  for  the 
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catalyst,  this  was  not  considered  to  have  affected  the 
other  results  in  any  way.  Because  this  deactivation 
rate  varied  inversely  with  space  time  (and  with  hydro¬ 
gen  sulfide  partial  pressure),  it  was  concluded  that 
coking  was  not  the  cause  of  this  deactivation.  It  is 
suggested  instead  that  a  loss  of  sulfur  from  the  cata¬ 
lyst  may  be  a  more  reasonable  explanation,  suggesting 
that  the  loss  in  activity  could  be  significantly 
reduced  by  a  slight  increase  in  the  hydrogen  sulfide 
partial  pressure  in  the  reactor. 

Because  the  equilibria  between  benzothiophene  and  dihy- 
drobenzothiophene  and  between  H2S  in  the  gas  phase  and 
sulfur  on  the  catalyst  surface  are  strongly  dependent 
on  pressure,  data  from  atmospheric  studies  should  not 
be  extrapolated  to  the  high-pressure  trickle-bed 
reactors  used  in  industry. 


9.  RECOMMENDATIONS 


1.  Analysis  of  the  intermediates  styrene  and  dihydrobenzo- 

thiophene  would  be  improved  using  a  gas  sampling  valve 
with  a  larger  internal  volume  than  was  used  in  the  pre¬ 
sent  study,  preferably  at  least  10  /il,  and/or  a  more 
sensitive  gas  chromatograph. 

2.  The  limitations  introduced  by  inconsistent  liquid  feed 

concentrations  could  be  lessened  if  the  apparatus  were 
redesigned  to  allow  refilling  of  the  syringe  without 
exposing  the  solution  to  the  atmosphere.  One  possible 
modification  to  the  apparatus  may  be  to  include  a 
three-way  valve  in  the  line  between  the  syringe  and  the 
hot  box,  the  third  arm  of  which  is  connected  to  a  large 
reservoir  of  feed  solution  kept  under  pressure,  and 
from  which  the  syringe  could  be  refilled  once  it  was 
empty.  This  modification  would  have  the  additional 
advantage  that  a  more  dilute  feed  solution  could  be 
used.  As  a  result,  the  heat  capacity  of  the  gas  stream 
in  the  reactor  would  be  increased,  and  this  would 
lessen  any  possible  effects  of  heat  transfer  limita¬ 
tions  on  the  reaction  rate. 

3.  The  reaction  mechanism  for  the  hydrodesulfurization  of 

benzothi ophene  should  be  studied  in  more  detail.  If 
dihydrobenzothiophene  and  styrene  are  included  in  the 
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feed,  this  would  assist  in  obtaining  more  data  on  the 
rates  of  the  individual  reactions  making  up  the  overall 
network . 

4.  The  cause  of  the  long  term  catalyst  deactivation  should 
also  be  investigated  further.  The  addition  of  small 
amounts  of  hydrogen  sulfide  to  the  gas  stream  should 
show  whether  or  not  this  compound  is  involved  in  the 
deactivation,  although  spectroscopic  measurements  of 
the  catalyst  itself  would  probably  be  necessary  to  draw 
any  conclusions  concerning  the  actual  mechanism  of  the 
deactivation . 
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APPENDIX  A  -  ROTAMETER  CALIBRATION  DATA 


The  hydrogen  and  nitrogen  rotameters  were  calibrated 
using  a  stopwatch  (accurate  to  0.1  s)  and  a  soap-bubble 
flowmeter  manufactured  from  a  50  ml  burette.  The  procedure 
consisted  of  setting  a  flow  rate  on  the  rotameter,  adjusting 
the  pressure  to  101  kPa  using  a  globe  valve  located  down¬ 
stream  of  the  rotameter,  and  measuring  the  length  of  time  it 
took  for  a  given  volume  of  gas  to  pass  through  the  flow¬ 
meter.  This  volume  was  usually  chosen  so  that  the  time  meas¬ 
ured  would  be  on  the  order  of  30  s,  although  when  the  flow 
rate  was  very  high  and  the  full  volume  of  the  flowmeter  was 
used,  the  measured  times  were  as  low  as  5  s. 

The  calibration  data  are  presented  in  Table  A-1.  The 
volumetric  flow  rates  shown  were  measured  at  24°C  and  92.6 
kPa.  The  following  polynomials  were  fitted  to  these  data 
using  a  least  squares  algorithm*0: 


Qhg  =  1 . 6258x 1  0 "  1  +  1.3138x10-3R  +  1 . 1 1 97x 1 0 - 4R2 
-  2. 1589x10- 7R3 


( A- la ) 


Qhs  =  7.2424  +  3 . 70 1 7x 1 0" 2R  +  1 . 8662x 1 0 ' 4R2 
+  4. 1 18x10" 7R3 


(A- 1b) 


Qng  =  8.51 57x 1 0 " 2  +  6.4351x10‘4R  +  8 . 2782x 1 0 ‘ 5R2 
-  5.7607x10' 7R3  +  1 .8400x10- ’R4 


( A- 1c ) 


Qns  =  -3.0917  +  1 . 4885x 1 0" 1 R  -  2 . 1 1 37x 1 0 " 3Ra 
+  1 .4609x10' 6R3  -  3.5941x10' *R4 


(A- Id) 


where 


Q  =  Volumetric  flow  rate,  cm3/s 
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R  =  Rotameter  reading 
h  =  Hydrogen  Rotameter 
n  =  Nitrogen  Rotameter 
g  =  Glass  Float 
s  =  Steel  Float 

The  order  of  the  polynomial  was  the  minimum  such  that 
the  addition  of  another  term  would  not  give  a  statistically 
better  fit  to  the  data.  Both  the  data  and  the  fitted  polyno¬ 
mials  are  shown  graphically  in  Figure  A- 1 . 
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Table  A-la  -  Calibration  Data  for  Hydrogen  Rotameter 

(glass  ball) 


Rotameter 

Flow  Rate 

Reading 

cm3  /s 

25.0 

0.2503 

25.5 

0.2519 

26.0 

0.2407 

25.5 

0.2519 

29.0 

0.3035 

29.0 

0.2869 

28.5 

0.2813 

28.5 

0.2703 

33.3 

0.3231 

32.0 

0.3211 

32.0 

0.3247 

30.0 

0.3284 

40.5 

0.3752 

40.0 

0.3953 

40.0 

0.4024 

40.0 

0.3839 

51.3 

0.4890 

51.3 

0.5195 

51.7 

0.5168 

52.5 

0.4975 

52.5 

0.5102 

59.5 

0.5634 

59.7 

0.5848 

59.8 

0.5747 

59.6 

0.5682 

67.  1 

0.6515 

66.4 

0.6601 

66.0 

0.6515 

66.2 

0.6410 

65.8 

0.6515 

74.2 

0.7605 

73.4 

0.7519 

73.4 

0.7463 

72.8 

0.7605 

73.2 

0.7547 

72.0 

0.7463 

71.5 

0.7519 

80.0 

0.8389 

79.2 

0.8451 

78.2 

0.8446 

77.2 

0.8361 

84.5 

0.9146 

83.6 

0.9202 

82.4 

0.9615 

82.5 

0.9375 

91.3 

1.0169 

Rotameter 

Flow  Rate 

Reading 

cm 3  /s 

91.8 

1 .0000 

91.3 

1 .0067 

90.8 

1.0135 

90.5 

1.0239 

90.6 

1 .0239 

90.0 

0.9967 

90.5 

0.9772 

101.0 

1.1364 

100.8 

1.1173 

101.2 

1.1236 

100.5 

1.1205 

100.7 

1.1364 

107.0 

1  .  2270 

106.0 

1 .2346 

105.8 

1.2461 

105.5 

1 .2422 

111.8 

1.3333 

110.5 

1.3513 

111.0 

1.3468 

110.6 

1 .3559 

110.0 

1 .3559 

110.0 

1 .3559 

120.8 

1 .4970 

120.0 

1.5198 

119.5 

1.5106 

119.5 

1.5060 

119.5 

1.5106 

118.2 

1 .4620 

118.4 

1 .4620 

118.4 

1 .4663 

118.6 

1 .4620 

124.3 

1 .6026 

124.3 

1.5873 

124.2 

1  .  5974 

124.2 

1.5974 

124.2 

1 .5823 

124.2 

1  .5924 

133.0 

1 .6949 

133.0 

1.7065 

133.0 

1 .7007 

132.7 

1.7065 

132.7 

1 .7065 

132.7 

1.7123 

140.8 

1.9011 

140.8 

1.9011 

140.8 

1.8797 

. 
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Table  A- la  (cont . ) 


Rotameter 

Flow  Rate 

Rotameter 

Flow  Rate 

Reading 

cm3/s 

Reading 

cm3/s 

140.8 

1.9011 

135.5 

1 .7483 

140.8 

1.9084 

135.5 

1 .7422 

150.0 

2.0492 

114.8 

1 .4205 

150.0 

2.0243 

114.8 

1 .4045 

150.0 

2.0243 

114.8 

1 .4045 

150.0 

2.0492 

114.8 

1  .4006 

150.0 

2.0325 

95.  1 

1.1390 

145.6 

1 .9084 

94.9 

1.1359 

145.6 

1 .9084 

95.0 

1.1013 

145.6 

1 .8727 

95.0 

1.1312 

135.5 

1 .7606 

. 
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Table  A-lb  Calibration  Data  for  Hydrogen  Rotameter 

(steel  ball) 


Rotameter 

Flow  Rate 

Reading 

cm3/s 

28.0 

1  .  9305 

22.5 

1  .  6340 

23.0 

1  .  6287 

23.0 

1  .  6234 

24.0 

1  .  5924 

22.0 

1.6129 

30.0 

2.0243 

30.0 

2.0243 

29.5 

2.0000 

29.2 

1.9608 

28.5 

2.0000 

52.5 

3.2051 

52.5 

3.1250 

53.0 

3.1646 

52.5 

3.1646 

53.0 

3.1250 

37.5 

2.3474 

36.0 

2.3697 

36.0 

2.3364 

34.5 

2.3364 

33.5 

2.2624 

34.0 

2.3148 

71.5 

4.3860 

71.3 

4.3860 

71.3 

4.3860 

71.6 

4.3478 

71  .4 

4.3103 

71.2 

4.3478 

71.0 

4.3860 

79.8 

5.1546 

79.5 

5.0505 

79.3 

4.8544 

78.0 

5.0505 

77.5 

5.1020 

76.0 

5.0505 

76.0 

5.0505 

86.2 

5.6180 

86.3 

5.5556 

86.0 

5.5556 

85.5 

5.4945 

84.5 

5.5556 

81.0 

5.6180 

83.0 

5.5556 

94.0 

6.0976 

93.0 

6.0976 

92.5 

6.0976 

Rotameter 

Flow  Rate 

Reading 

cm3/s 

93.0 

6.1728 

93.5 

6.1728 

93.0 

6.1728 

107.5 

7.3529 

107.3 

7.3529 

107.3 

7.2464 

107.2 

7.2464 

106.5 

7.3529 

106.0 

7.2464 

120.5 

8.6207 

120.0 

8.6207 

120.0 

8.7719 

120.0 

8.4746 

119.0 

8.4746 

119.0 

8.6207 

118.5 

8.4746 

130.8 

9.8039 

130.8 

9.8039 

130.8 

9.8039 

129.5 

9.6154 

129.6 

9.8039 

130.0 

9.8039 

150.0 

1  1  .6279 

150.0 

1  1  .9048 

150.0 

1  1 .9048 

150.0 

1 1 .9048 

145.0 

1 1 .6279 

145.0 

1 1 .6279 

145.0 

1 1 .6279 

139.5 

10.8696 

139.5 

10.2041 

139.5 

10.2041 

139.5 

10.6383 

139.0 

10.4167 

135.0 

10.2041 

134.8 

9.8039 

134.8 

10.0000 

134.7 

9.8039 

134.2 

9.8039 

126.5 

9.2593 

126.3 

8.9286 

125.5 

9.6154 

125.5 

9.4340 

125.0 

9.0909 

125.0 

9.0909 

114.5 

7.9365 

• 
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Table  A- 1 b  ( cont . ) 


Rotameter 

Flow  Rate 

Rotameter 

Flow  Rate 

Reading 

cm 3  /s 

Reading 

cm3/s 

114.5 

7.6923 

57.5 

3.4722 

114.5 

7.8125 

57.5 

3.4722 

114.0 

7.8125 

57.0 

3.4965 

113.8 

7.8125 

58.0 

3.4965 

111.2 

7.8125 

58.0 

3.4965 

111.2 

7.6923 

47.3 

2.8249 

111.2 

7.6923 

47.8 

2.8249 

111.2 

7.6923 

47.5 

2.8409 

111.2 

7.8125 

47.3 

2.8249 

102.2 

6.7568 

47.5 

2.8409 

101.8 

6.7568 

39.8 

2.5381 

101.8 

6.7568 

40.5 

2.5510 

101.5 

6.6667 

38.8 

2.5381 

101.5 

6.7568 

44.0 

2.8090 

89.0 

5.8139 

43.0 

2.7933 

88.5 

5.6818 

44.0 

2.8090 

88.8 

5.8823 

34.0 

2.2831 

88.2 

5.6818 

35.0 

2.4390 

63.5 

4.0323 

36.0 

2.4390 

63.0 

4.0000 

26.0 

1 .8050 

64.5 

4.0650 

26.5 

1 .8248 

65.0 

4.0650 

26.0 

1.8315 

64.5 

4.0650 

26.0 

1 .8382 

57.5 

3.4722 

• 
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Table  A-lc  -  Calibration  Data  for  Nitrogen  Rotameter 

(glass  ball) 


Rotameter 

Flow  Rate 

Reading 

cm3/s 

19.6 

0.1247 

19.6 

0. 1289 

19.6 

0.1012 

19.6 

0.1572 

19.6 

0.1359 

24.2 

0. 1433 

24.2 

0.1401 

24.2 

0.1335 

24.2 

0.1292 

24.2 

0.1276 

36.8 

0. 1776 

37.0 

0.1799 

36.2 

0.1789 

40.8 

0.2347 

40.7 

0.2257 

40.0 

0.2262 

40.2 

0.2262 

40.2 

0.2247 

40.0 

0.2222 

39.5 

0.2217 

43.0 

0.2242 

43.4 

0.2247 

42.5 

0.2247 

50.0 

0.2632 

50.0 

0.2611 

49.5 

0.2611 

49.3 

0.2691 

49.0 

0.2558 

48.7 

0.2545 

48.2 

0.2545 

58.5 

0.3077 

58.2 

0.3058 

58.2 

0.3040 

57.9 

0.3049 

57.5 

0.3049 

57.2 

0.3040 

63.8 

0.3448 

63.3 

0.3460 

62.8 

0.3413 

62.5 

0.3378 

62.2 

0.3367 

70.0 

0.3759 

69.4 

0.3745 

69.3 

0.3731 

69.0 

0.3717 

68.5 

0.3717 

Rotameter 

Flow  Rate 

Reading 

cm3/s 

68.2 

0.3717 

94.0 

0.5348 

93.3 

0.5333 

92.9 

0.5333 

92.6 

0.5333 

92.2 

0.5333 

100.0 

0.5848 

100.0 

0.5848 

99.8 

0.5848 

99.8 

0.5848 

99.5 

0.5865 

104.8 

0.6192 

104.4 

0.6173 

104.2 

0.6173 

104.2 

0.6173 

104.0 

0.6173 

111.0 

0.6667 

110.6 

0.6667 

110.2 

0.6667 

110.1 

0.6667 

110.0 

0.6667 

115.4 

0.7067 

115.3 

0.7067 

115.0 

0.7042 

115.0 

0.7017 

1  14.8 

0.7042 

119.4 

0.7407 

119.8 

0.7353 

123.1 

0.7576 

122.9 

0.7576 

122.6 

0.7605 

122.5 

0.7634 

122.2 

0.7634 

130.0 

0.8163 

129.9 

0.8231 

129.6 

0.8197 

129.5 

0.8231 

129.6 

0.8197 

134.6 

0.8645 

134.0 

0.8571 

134.0 

0.8596 

133.8 

0.8571 

133.6 

0.8596 

140.3 

0.9231 

139.2 

0.9174 

138.9 

0.9174 
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Table  A- 1c  (cont . ) 


Rotameter 

Flow  Rate 

Rotameter 

Flow  Rate 

Reading 

cm3/s 

Reading 

cm 3  /s 

138.8 

0.9293 

141.7 

0.9464 

138.8 

0.9464 

149.2 

1  .0274 

142.2 

0.9464 

149.0 

1.0135 

142.0 

0.9434 

148.6 

1 .0204 

141.9 

0.9434 

148.7 

1.0169 

141.8 

0.9464 

148.8 

1.0135 

141.7 

0.9464 
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Table  A-ld  Calibration  Data  for  Nitrogen  Rotameter 

(steel  ball) 


Rotameter 

Flow  Rate 

Reading 

cm3/s 

59.8 

0.8850 

59.5 

0.8876 

59.0 

0.8902 

59.7 

0.8902 

58.2 

0.8902 

64.0 

0.9901 

63.3 

0.9868 

63.  1 

0.9836 

62.7 

0.9868 

61.8 

0.9836 

68.5 

1.0714 

67.8 

1 .0753 

67.0 

1.0791 

66.3 

1 .0753 

65.5 

1 .0753 

74.2 

1.1494 

73.7 

1.1494 

72.9 

1.1538 

72.3 

1.1493 

71.8 

1.1538 

80.4 

1 .3072 

79.5 

1 .3029 

78.8 

1 .2987 

78.6 

1  .  3029 

78.3 

1  .  3029 

84.7 

1 .3605 

84.2 

1 .3559 

83.3 

1 .3605 

82.6 

1 .3605 

82.0 

1 .3605 

91.8 

1.4981 

91.7 

1.4981 

91.4 

1.4981 

91  .3 

1.4925 

91.1 

1 .4925 

96.9 

1 .6026 

96.7 

1 .5924 

96.6 

1.5974 

96.3 

1  .5974 

96.2 

1 .6026 

99.5 

1 .6502 

99.5 

1  .6447 

99.4 

1 .6556 

99.5 

1 .6556 

99.3 

1.6556 

105.6 

1.8182 

Rotameter 

Flow  Rate 

Reading 

cm3/s 

105.7 

1.8116 

105.8 

1.8182 

105.8 

1.8182 

105.9 

1 .7986 

109.8 

1 .8868 

109.9 

1 .8797 

109.9 

1  .8657 

109.9 

1.8797 

109.9 

1  .8727 

114.6 

1 .9920 

114.5 

1 .9685 

114.3 

1.9685 

114.3 

1.9608 

114.1 

1 .9608 

120.5 

2.1368 

120.3 

2.1097 

120.1 

2.1097 

120.1 

2.1186 

120.0 

2.0920 

124.0 

2.2624 

124.1 

2.2422 

124.0 

2.2522 

124.2 

2.2624 

124.3 

2.2422 

129.4 

2.3697 

129.1 

2.3809 

129.6 

2.3585 

129.5 

2.3148 

129.5 

2.4510 

129.7 

2.3697 

134.5 

2.5641 

145.9 

2.6455 

145.8 

2.6738 

145.8 

2.6455 

145.7 

2.6455 

145.7 

2.6455 

149.9 

2.8249 

149.9 

2.8409 

149.9 

2.8409 

149.9 

2.8249 

149.9 

2.8409 

140.6 

2.5510 

140.5 

2.5126 

137.5 

2.5510 

137.5 

2.5641 

137.4 

2.5510 
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Table  A- Id  (cont . ) 


Rotameter  Flow  Rate 
Reading  cm3/s 


Rotameter  Flow  Rate 
Reading  cm3/s 


137.5  2.5907  137.7  2.5641 

137.7  2.5641 
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Figure  A- la  -  Calibration  Curve  for  Hydrogen  Rotameter  (glass  ball) 
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Figure  R— lb  -  Calibration  Curve  for  Hydrogen  Rotameter  (steel  ball) 
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Figure  R-ld  -  Calibration  Curve  for  Nitrogen  Rotameter  (steel  ball) 


APPENDIX  B  -  SYRINGE  PUMP  CALIBRATION  DATA 


The  procedure  for  calibrating  the  syringe  pump  was  as 
follows:  The  syringe  was  filled  with  water  and  connected  to 
the  pump.  The  section  of  1/16"  tubing  which  was  to  lead  from 
the  syringe  to  the  flowing  gas  stream  during  the  actual 
experimentation  was  attached  to  the  syringe.  A  flow  rate  was 
set  and  the  pump  was  allowed  to  run  for  a  short  while  for 
the  flow  rate  to  stabilize.  A  weighed  flask  was  set  under 
the  end  of  the  tubing  and  water  was  collected  for  a  given 
period  of  time,  measured  with  a  stopwatch.  The  length  of  the 
time  period  was  adjusted  so  that  the  mass  of  the  liquid  col¬ 
lected  was  on  the  order  of  a  few  grams.  Since  the  flow  of 
water  fell  into  the  flask  as  discrete  drops  rather  than  as  a 
continuous  flow,  the  time  period  started  and  finished  just 
after  a  drop  was  released  from  the  end  of  the  tubing. 

The  flow  rates  measured  fell  between  10%  x  1/1000  and 
100%  x  1/100  of  the  maximum  flow  rate  of  the  pump;  these 
data  are  reported  in  Table  B-1  and  plotted  in  Figure  B-1. 
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Table  B-1  _  Syringe  Pump  Calibration  Data 


Flow 

Setting 

Flow  Rate 

%  Max 

.  Rate 

ml/min . 

10% 

X 

1/1000 

0.0088 

20% 

X 

1/1000 

0.0172 

30% 

X 

1/1000 

0.0262 

40% 

X 

1/1000 

0.0335 

50% 

X 

1/1000 

0.0443 

60% 

X 

1/1000 

0.0520 

70% 

X 

1/1000 

0.0596 

80% 

X 

1/1000 

0.0690 

90% 

X 

1/1000 

0.0768 

100% 

X 

1/1000 

0.0881 

10% 

X 

1/100 

0.0894 

20% 

X 

1/100 

0.1798 

30% 

X 

1/100 

0.2639 

40% 

X 

1/100 

0.3513 

50% 

X 

1/100 

0.4390 

60% 

X 

1/100 

0.5182 

70% 

X 

1/100 

0.6103 

80% 

X 

1/100 

0.6998 

90% 

X 

1/100 

0.8026 

15% 

X 

1/100 

0.1337 

25% 

X 

1/100 

0.2201 

35% 

X 

1/100 

0.3057 

45% 

X 

1/100 

0.3954 

55% 

X 

1/100 

0.4851 

65% 

X 

1/100 

0.5833 

75% 

X 

1/100 

0.6574 

85% 

X 

1/100 

0.7622 

95% 

X 

1/100 

0.8291 

100% 

X 

1/100 

0.8758 

0.90^ 
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Figure  B-1  -  Syringe  Pump  Calibration  data 


APPENDIX  C  -  GAS  CHROMATOGRAPH  DETECTOR  CALIBRATION  DATA 

Prior  to  analysis,  the  gas  chromatograph  detector  was 
calibrated  for  the  three  major  components  in  the  product 
stream,  benzothiophene ,  ethylbenzene,  and  styrene,  using 
n-hexane  as  an  internal  standard. 

The  procedure  was  identical  for  all  three  components. 
First,  a  small  volume  (approximately  50  ml)  of  a  solution  of 
known  composition  of  the  component  in  n-hexane  was  made. 

Five  samples  of  10  fil  each  were  injected  into  the  gas  chro¬ 
matograph,  the  response  ratio  of  the  component  to  n-hexane 
was  calculated  for  each  sample,  and  the  mean  of  the  five 
samples  was  computed.  This  procedure  was  repeated  until  six 
or  seven  data  points  had  been  collected. 

For  all  three  components,  the  response  ratio  varied 
linearly  with  the  weight  ratio.  The  data  are  summarized  in 
Table  C~ 1  and  plotted  in  Figure  C-1.  The  relative  response 
factors  are  summarized  in  Table  C-2. 

The  only  relative  response  factor  available  in  the  lit¬ 
erature  was  that  for  ethylbenzene81,  1.17.  No  explanation 
can  be  offered  for  the  discrepancy  between  this  value  and 
the  value  obtained  experimentally  (0.99). 

The  analysis  was  performed  under  the  same  conditions  as 
described  in  Chapter  3,  except  that  the  temperature  program¬ 
ming  rate  was  15°/min.  for  ethylbenzene  and  styrene  and 
20°/min.  for  benzothiophene,  and  that  the  run  was  terminated 
as  soon  as  the  component  of  interest  had  been  fully  eluted. 
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Table  C-1  -  Detector  Calibration  Data 


Component 

Styrene 

Ethylbenzene 

Styrene 


Weight  Ratio 
Component/HEX 


0.0543 

0.0188 

0.0067 

0.0857 

0.0261 

0.1092 


0.1814 

0.0208 

0.0134 

0.0088 

0.1396 

0.0704 

0.0511 


0.0374 

0.0181 

0.0050 

0.0728 

0.1063 

0.1637 


Response  Ratio 
Component/HEX 


0.0562 

0.0194 

0.0067 

0.0888 

0.0254 

0.1098 


0.1823 

0.0222 

0.0139 

0.0089 

0.1401 

0.0706 

0.0521 


0.0364 

0.0174 

0.0045 

0.0711 

0.1040 

0.1648 
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WEIGHT  RATIO,  STY/HEX 

Figure  C-la  -  Calibration  Data  for  Styrene 
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0.15  0.20 

WEIGHT  RATIO,  BT/HEX 

Figure  C-lb  -  Calibration  Data  for  Benzothiophene 


0.20-1 
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Figure  C-lc  -  Calibration  Data  for  Ethylbenzene 
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Table  C-2  -  Relative  Response  Factors 


Component 


Relative 
Response  Factor 


Benzothiophene/n-Hexane 
Ethy 1 benzene /n -Hexane 
Sty rene/n -Hexane 


1  .0272 
0.9942 
0.9697 


APPENDIX  D  -  DENSITY-COMPOSITION  DATA  FOR  LIQUID  FEED 

SOLUTIONS 

This  appendix  describes  work  done  to  verify  that  the 
specific  gravity  of  solutions  of  benzothiophene  in  n-hexane 
could  be  predicted  by  simply  adding  the  products  of  the 
weight  fractions  and  specific  gravities  of  each  pure  compo¬ 
nent  . 

A  total  of  four  data  points  were  collected.  In  each 
case  a  small  volume  (about  50  ml)  of  solution  was  mixed  on  a 
balance  and  immediately  used  to  fill  a  25  ml  specific  grav¬ 
ity  bottle.  The  bottle  was  weighed,  and  the  specific  gravity 
calculated.  The  weight  ratios  of  benzothiophene  to  n-hexane 
in  the  test  solutions  and  the  experimental  and  theoretical 
specific  gravities  are  shown  in  Table  D-1.  The  experimental 
points  and  the  theoretical  curve  are  shown  in  Figure  D-1. 

The  specific  gravity  of  pure  n-hexane  is  0.6603*  2  ;  the 
specific  gravity  of  pure  solid  benzothiophene  is  1.14975. 

The  experimental  conditions  were  20  °C  and  93.3  kPa. 
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Table  D- 1  -  Density-Composition  Data  for 

Benzothiophene/n-Hexane  Solutions 


Wt .  Ratio 
BT/HEX 

0.7334 

0.6206 

0.5521 

0.3974 


Experimental 

S.G. 

0.8648 

0.8524 

0.8302 

0.8007 


Theoretical 

S.G. 

0.8671 

0.8474 

0.8341 

0.7993 


0.90H 
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APPENDIX  E  -  LISTING  OF  COMPUTER  PROGRAMS 


The  main  data  analysis  program  used  in  the  preparation 
of  this  thesis  was  entitled  THEUNIVERSE,  listed  beginning  on 
the  next  page.  It  was  used  to  calculate  the  fractional  con¬ 
versions  of  benzothiophene  to  each  of  the  major  products  and 
to  compute  an  estimate  of  the  composition  of  the  liquid 
feed.  A  typical  output  from  this  program  is  given  in  Table 
3. 

The  program  used  to  fit  a  least  squares  exponential  to 
the  fractional  conversion  data  making  up  each  individual 
data  point  was  called  REGPT,  and  is  listed  after 
THEUNIVERSE.  The  constants  returned  by  this  program  are 
given  in  Table  4,  and  the  curves  themselves  are  plotted  in 
Figures  5a  through  5p.  The  actual  output  from  this  program 
can  be  obtained  from  the  author  on  request. 

The  library  subroutine  ("AR")  called  by  this  program 
has  been  documented  elsewhere74.  The  function  to  which  the 
data  are  fitted  (Equation  2)  is  passed  to  the  program  as  a 
single  line  of  FORTRAN  in  a  separate  file. 
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C  THEUNI VERSE 
C 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

C  The  input  to  the  program  for  each  series  of  gas 
C  samples  is  as  follows: 

C  NSER:  Identification  number  for  the  series 

C  NBED :  Identification  number  for  the  catalyst  bed 

C  NFEED :  Identification  number  for  the  liquid  feed 

C  batch 

C  NDATE :  Date  on  which  the  data  were  obtained 

C  NPTS:  Number  of  gas  samples  in  the  series 

C  BTCONC :  Wt .  ratio  bt/hex  in  the  liquid  feed 

C  (Although  this  program  will  calculate  a  value 

C  based  on  a  composite  mass  balance  of  all  gas 

C  samples  in  the  data  file,  this  input  value 

C  will  be  the  one  used  to  calculate  bt  and  hex 

C  flow  rates.  In  other  words,  this  program  must 

C  be  run  twice  for  any  given  set  of  data  -  once 

C  to  find  out  the  value  of  BTCONC,  and  a 

C  second  time  to  use  it  in  the  calculations.) 

C  ATM:  Average  atmospheric  pressure  for  the  series, 

C  mm.  Hg 

C  P:  Average  reactor  pressure  for  the  series,  inches 

C  Meriamfluid(s.g.  1.75) 

C  TC:  Average  reactor  temperature,  deg.  C 


This  is  the  main  data  analysis  program  used  in  the 
present  thesis.  It  takes  all  the  data  recorded 
during  the  actual  experimentation  as  input,  and 
outputs  the  following  for  each  series  of  gas  samples: 
Series  identification  number 
Date 

Reactor  pressure 
Reactor  temperature 
Catalyst  bed  mass 
Catalyst  pellet  mesh  size 
Catalyst  bed  identification  number 
Liquid  feed  batch  identification  number 
Molar  flow  ratesof  all  feed  components 
Partial  pressures  of  all  feed  components 
Reciprocal  space  velocity 
For  each  individual  gas  sample,  the  following  are 
printed  out: 

Percent  of  total  peak  area  for  each  component 
in  the  product  stream 
A  hypothetical  feed  concentration  based 
on  a  mass  balance  for  that  gas  sample 
Fractional  conversion  of  bt  to  each  of  the 
reaction  products 

Finally,  the  program  prints  the  means  and  95% 
confidence  intervals  of  the  fractional  conversions 
over  the  series  of  gas  samples,  and  gives  an 
estimate  of  the  original  liquid  feed  composition 
based  on  a  composite  mass  balance  over  the  entire 
liquid  feed  batch. 
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C  CATWT:  Catalyst  bed  weight,  g. 

C  MESH1,MESH2:  Catalyst  pellet  size,  standard  mesh 

C  (e.g.  -40+60) 

C  H2 :  Hydrogen  rotameter  reading 

C  NHROT:  Hydrogen  rotameter  float  (1=glass,  2=steel) 

C  SP:  Syringe  pump  setting,  %  x  1/1000  of  full  rate 

C  EN2:  Nitrogen  rotameter  reading 

C  NNROT:  Nitrogen  rotameter  float  (as  above) 

C  In  addition,  the  following  data  are  inputted  for 
C  each  gas  sample: 

C  THR,TMIN:  Time  elapsed  between  reactor  going 

C  on  stream  and  sampling  time,  in  hours  and 

C  minutes 

C  GAS:  Sum  of  peaks  for  hydrogen,  nitrogen,  and 

C  hydrogen  sulfide,  %  of  total  area  of  all 

C  peaks  (not  used  in  any  calculations) 

C  HEX:  N-hexane  peak  size,  %  of  total  area 

C  EB:  Ethylbenzene  peak  size,  %  of  total  area 

C  STY:  Styrene  peak  size,  %  of  total  area 

C  BT:  benzothiophene  peak  size,  %  of  total  area 

C  DH:  Dihydrobenzothiophene  peak  size,  %  of  total 

C  area 

C 

DIMENSION  MESH ( 2 ) ,GAS(80) , HEX (80) ,EB(80) , STY (80) , 

<5>  CONV (80)  ,  BTTOEB  (  80  )  ,BTTOST(80)  ,T(70)  ,BALNCE(80)  , 

@  BTTODH ( 80 ) ,THR(80) ,TMIN(80) , BT ( 80 ) , DH ( 80 ) 

Initialize  variables 

DATA  T/ 12. 71, 4. 30, 3. 18, 2. 78, 2. 57, 2. 45, 2. 36, 2. 31, 

<S>  2.26,2.23,2.2,2.  18,2.16,2.  14,2.13,2.  12,2.  1  1  , 

<6>  2.  1  ,2*2.09,2.08,2*2.07,3*2.06,2*2.05,2*2.04, 

@  5*2.03,8*2.02,12*2.01,16*2./ 

TOTCOM=  0 . 

TOTSQ=  0 . 

NCOMP=  0 
NNROT= 1 

Input  data  for  a  new  series  of  runs  and  re-initialize 
variables 

1 0 0  READ ( 5 , 1 0 1 ) NSER , NBED , NFEED , NDATE , NPTS , BTCONC , ATM , P , 

$  TC, CATWT, MESH ( 1 ) ,MESH(2) , H2 , NHROT , SP , EN2 , NNROT 

IF(NSER.EQ.O)  GO  TO  200 
CONVST=0 . 0 
CONVDH=0 . 0 
BTDHT=0 . 

BTDHT2=0 . 

BTEBT=0 . 

BTEBT2=0 . 

BTSTT=0 . 

BTSTT2=0 . 

BTCT=0 . 

BTCT2=0 . 
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DO  110  1=1, NPTS 

READ ( 5 , 1 1 1 ) THR ( I ) ,TMIN(I ) ,GAS(I ) ,HEX(I ) , EB ( I ) , 

§  STY ( I ) , BT ( I ) , DH ( I ) 

THR ( I ) =THR ( I ) +TMI N ( I ) /6  0 . 

1 10  CONTINUE 

Calculate  reactor  pressure  (kPa)  and  temperature 
(deg.  K) 

PR=P*0 . 4358935+ATM* . 1333224 
TK=TC+273 . 16 

Calculate  volumetric  (Q)  and  molar  (F)  flow  rates  for 
gas  (cm3/s, mol/s)  and  liquid  (ml/min,  mol/s)  feeds; 
Calculate  reciprocal  space  velocity 

I F ( NHROT . EQ . 1 )  QH2=1 . 6258 1 56 1 5E- 1 + 1 . 3 1 3837578E-3*H2 
<6>  +1  .  1  19669541E-4*H2**2-2.518931661E-7*H2**3 

IF (NHROT. EQ. 2)  QH2=7 . 242407802E- 1 +3 . 70 1 69792 1E-2*H2 
<P>  +1  ,866222346E-4H2**2+4. 1  1 888704 6E-7*H2** 3 

I F ( NNROT . EQ . 1 )  QN2=8 . 5 1 565 360 7E- 2+6 . 43509775 1E-4*EN2 
@  +8.2781 9538 1E-5*EN2** 2-5. 76040392E-7*EN2**3 

@  +1 .83998824E-9*EN2**4 

IF (NNROT. EQ. 2)  QN2=-3 . 09 1 740296+ 1 . 488 5233 68E- 1 *EN2 
(P>  -2.  1  1  37  1  6022E-3*EN2**2  +  1  . 46093048E~5*EN2**3 

<5>  -3 . 594074268E-8*EN2**4 

FH2=3.7476E-5*QH2 
FN2=3 . 7476E~5*QN2 
I F ( EN2 . EQ . 0 )  FN2=0 . 

BTFRAC=BTCONC/( 1 . +BTCONC ) 

QLIQ=9.678E-4*SP 
FBT=BTFRAC*QLIQ* 1 .42697E-4 
FHEX= ( 1 . -BTFRAC ) *QLIQ* 1 .2769E-4 
WF=CATWT/FBT 

Calculate  total  molar  flow  rate  and  partial  pressures 
of  each  feed  component  (kPa) 

F=FH2+FBT+FHEX+FN2 

PH2=PR*FH2/F 

PBT=PR*FBT/F 

phex=pr*fhex/f 

PN2=PR*FN2/F 

Calculate  fractional  conversion  of  bt  to  eb,  sty, 
dhbt  for  each  run 

First,  a  weight  ratio  and  then  a  mole  ratio  of  each 
component  to  n-hexane  is  calculated 

DO  120  1=1 , NPTS 

WRBTHX=BT ( I )*  1  . 027  1 93/HEX(  I  ) 

WREBHX=EB ( I ) *0.99421 1/HEX(I ) 

WRSTHX=STY ( I ) *0 . 9 6 97 3 9/HEX ( I ) 
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WRDHHX=DH ( I ) *  1  .  0272/HEX(l ) 

BTMOL=WRBTHX*86. 18/134.2 
EBMOL=WREBHX*86 . 18/106.07 
STYM0L=WRSTHX*8 6 . 18/104.15 
DHMOL=WRDHHX*86 . 18/136.2 
TOTMOL = BTM0L+ E  BM0L+ STYMOL + DHMOL 
BALNCE ( I ) =T0TM0L* 134.2/86. 18 

Next,  a  benzene  ring  balance  on  the  product  stream 
is  computed.  TOTMOL  is  the  mole  ratio  of  total 
aromatics  to  n-hexane  in  the  product 
stream,  which  should  be  the  same  as  the  mole  ratio 
of  bt/hex  in  the  feed,  (the  weight  ratio  of  bt/hex 
in  the  original  feed  is  given  by  BALNCE.)  Based  on 
this  estimate  of  bt/hex  mole  ratio,  the  fractional 
conversions  of  benzothiophene  to  ethylbenzene, 
styrene,  and  dihydrobenzothiophene  are  calculated. 

BTTOEB ( I ) =EBMOL/TOTMOL 
BTTOST ( I ) =  STYMOL /TOTMOL 
BTTODH ( I ) = DHMOL /TOTMOL 
CONV ( I ) = BTTOEB ( I ) +BTTOST ( I ) 

The  rest  of  the  loop  is  used  to  store  intermediate 
values  of  variables  used  to  compute  means 
and  confidence  intervals  in  the  next  section. 

BTEBT=BTEBT+BTTOEB ( I ) 

BTSTT=BTSTT+BTTOST ( I ) 

BTCT=BTCT+CONV ( I ) 

BTDHT=BTDHT+ BTTODH ( I ) 

BTEBT2=BTEBT2+BTTOEB ( I ) **2 . 

BTSTT2=BTSTT2+BTTOST ( I ) **2 . 

BTCT2=BTCT2+CONV ( I ) **2 . 

BTDHT2=BTDHT2+BTTODH ( I ) **2 
NCOMP=NCOMP+ 1 
TOTCOM=TOTCOM+BALNCE ( I ) 

TOTSQ=TOTSQ+ BALNCE ( I ) **2 
120  CONTINUE 

Compute  mean  and  confidence  interval  at  .95  level 
of  significance 

CONVEB=BTEBT/NPTS 

convbt=btct/npts 

convst=btstt/npts 

CONVDH= BTDHT/NPTS 

SIGEBC= (BTEBT2/NPTS-CONVEB*CONVEB) **0 . 5 
SIGCON= (BTCT2/NPTS~CONVBT*CONVBT) **0 . 5 
SIGSTC= (BTSTT2/NPTS-CONVST*CONVST) **0 . 5 
SIGDHC= (BTDHT2/NPTS-CONVDH*CONVDH) ** . 5 
NU=NPTS-1 

CIEB=T(NU) *SIGEBC/NU**0 . 5 
Cl STY=T(NU) *SIGSTC/NU**0 . 5 


non  non 
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CITOT=T (NU ) *SIGCON/NU**0 . 5 
CIDH=T(NU) *SIGDHC/NU**0 . 5 


Output  everything 

WRITE (6,131)  NSER , NDATE 
WRITE (  6 , 132)  PR , TK 

WRI TE ( 6 , 1 3  3 )  CATWT , MESH ( 1 ) , MESH ( 2 ) 

WRI TE ( 6 , 139)  NBED , NFEED 
WRI TE (  6 , 134)  FH2 , FBT , FHEX , FN2 
WRITE (6, 135)  PH2 , PBT , PHEX , PN2 
WRI TE ( 6 , 138)  WF 
WRI TE ( 6 , 136) 

DO  130  1=1 , NPTS 

WRI TE ( 6 , 137)  I , GAS ( I ) ,HEX(I ) ,EB(l ) ,STY(I ) , BT ( I ) , 
<3>  DH  ( I  ) 

130  CONTINUE 

WRI TE ( 6 , 1 5 1 ) 

DO  150  1=1 , NPTS 

WRI TE (  6  , 152)  I , BALNCE ( I ) 

150  CONTINUE 

WRI TE ( 6 , 141 ) 

DO  140  1=1 , NPTS 

WRI TE ( 6 , 142)  I , THR ( I ) , BTTOEB ( I ) ,BTTOST(I ) , 

@  BTTODH ( I ) ,CONV(I ) 

WRI TE ( 7 , 1 4  2 )  I, THR ( I ) , BTTOEB ( I ) , BTTOST ( I ) , 

<6>  BTTODH  (I  )  ,  CONV  ( I  ) 

140  CONTINUE 

WRI TE ( 6 , 1 4  3 )  CONVEB , CONVST , CONVDH , CONVBT ,  CIEB , 

<3>  Cl  STY, Cl DH, Cl TOT 
NFEEDX=NFEED 
GO  TO  100 


Calculate  avg.  feed  composition  based  on  all  runs 

200  AVCOMP=TOTCOM/NCOMP 

S I GCOM= ( TOTSQ/NCOMP-AVCOMP* *  2 ) *  * . 5 
NUNU=NCOMP- 1 

Cl  COMP =T ( NUNU ) *S I GCOM/NUNU*  *0.5 

WRITE (6, 201 )  NFEEDX , NCOMP , AVCOMP , C I COMP 

STOP 

101  FORMAT (31 3, 1  7, 1  3, F7 . 4, F6. 1 ,F5 . 1 ,F6 . 1 , F8 . 5 , 21 3 ,F6 . 1 , 

@  I  2 , 2F6 .1,12) 

111  FORMAT (2F2.0,6F8. 4) 

131  FORMAT! ’ 1 ’,' SERIES  f ,13,/, 17,//) 

132  FORMAT! ’  ',’ TOTAL  PRESSURE  IN  REACTOR:  ’,F7.2,'  KPA ’ 

(5>  ,/,’  TEMPERATURE  IN  REACTOR:  ’  ,F6.1,’  DEG.  K') 

133  FORMAT!  ’  ’/CATALYST:  NALCO  NM  502,  ’  ,F7.4,’  G.  ,  ’, 

<3>  13, ’  +  ’,12,’  MESH’  ) 

139  FORMAT (’  ’,' CATALYST  BED  NUMBER:  ’,12,/, 

<5>  ’LIQUID  FEED  BATCH  NUMBER:  ’,12) 

134  FORMAT! ’  ’ , ’ FEED  RATES  TO  REACTOR:  ’,/,5X,’H2: 

<5>  El  1.6,’  GMOL/SEC  ’  ,  / ,  5X ,  ’  BT :  ’,E11.4,’  GMOL/SEC  ' 


•  - 


43 


■ 
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(5)  , /, 5X, ’  HEX:  '  ,  E 1  1  .  4  ,  '  GMOL/SEC '  ,  / ,  5X r  ’  N2 : 

<P>  El  1  .4,  ’  GMOL/SEC’  ) 

135  FORMAT ( ’  'PARTIAL  PRESSURES  OF  FEED  TO  REACTOR:’, 

<6>  / ,  5X ,  ’  H2  :  '  ,  F6 . 2  ,  ’  KPA  ’  ,  / ,  5X ,  ’  BT :  ’,F6.2, 

<3>  ’  KPA’  ,/,5X, ’HEX:  ’,F6.2,’  KPA '  ,  / ,  5X ,  '  N2  : 

<5>  F6.2, '  KPA’ ) 

1 38  FORMAT ( ’  ’,’W/F=  ’ , E 1 1 . 4 , ’  G~SEC/GMOL ( BT ) ' ,/) 

136  FORMAT ('  ’ , ’ RAW  DATA  (AREA  %):’,/, 5X RUN ’, 3X , 

@  ’  N2+H2  ’  ,2X,  'N-HEXANE'  ,4X,  ’ETHYL-’  ,3X,  ’STYRENE’  , 

@  4X ,  ’  BENZO-  DIHYDROBENZO'  ,/,  12X,  '+H2S*  ,  13X, 

<5>  ’  BENZENE  ’  ,  1  1 X ,  ’  THIOPHENE  '  ,  4X ,  '  THI OPHENE  ’  ) 

137  FORMAT ( ’  ’  , 4X , I  2 , 2X, 6 (F7 . 4 , 3X) ) 

151  FORMAT (’  ’,/, ’Ll  QUID  FEED  COMPOSITION  BASED  ON  MASS’, 

(P>  ’  BALANCE  (WT.  RATIO)  :’,/,  5X,  '  RUN ',  3X,  ’  BT/HEX’  ) 

152  FORMAT ('  ’  , 4X, I  2 , 3X, F8 . 6 ) 

141  FORMAT ( '  ',/, ’FRACTIONAL  CONVERSION  OF  BENZOTHIO’ , 

@  ’PHENE: ’ ,/,5X, ’RUN’ ,3X, ’TIME’ ,6X, ’TO  EB:’,4X, 

§  ’ TO  STY : ’ , 3X , ’ TO  DHBT : ’ , 4X , ’ EB+STY : ’ ) 

142  FORMAT ( '  ’  , 4X, I  2 , 3X , F6 . 3 , 4 ( 3X, F8 . 6 ) ) 

143  FORMAT (’  ’,3X,’MEAN:  ’ , 9X , 4 ( F8 . 6 , 3X ) , / , 5X , ’ +/- ’ , 1 1 X , 

<5>  4 ( F8 . 6 , 3X ) ) 

201  FORMAT( ’  1 ’, ’LIQUID  FEED  BATCH:  ’  , I  2 , // , ’ COMPOSI TE  ’, 

<5>  ’MASS  BALANCE  BASED  ON  ’,12,'  DATA  POINTS’,/, 

<6>  ’AVG.  FEED  RATIO  BT/HEX:  ’,F7.4,’  +/-  ’,F7.4) 

END 
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/PROB 

/INPUT 


/VAR 

/REGR 

/PARAM 

/END 


TITLE  IS  ?  REGPT ’  . 

VARIABLES  ARE  2. 

FORMAT  IS  ' (9X,F6.3,3X,F8.6) ' . 
CASES  ARE  54. 

MTSFILE  IS  ’ POINT  1 2 '  . 

NAMES  ARE  THR , CT . 

DEPEND  IS  CT. 

INDEP  IS  THR. 

PARAM  ARE  4. 

ITER  ARE  200. 

HALV  ARE  200. 

TITLE  IS  ’ POINT  12  '  . 

INITIAL  ARE  0.15,-0.002,0.1,-0.5. 
MAX  ARE  1.0, 0.0, 1.0, -0.3. 

MIN  ARE  0.0, -0.1, 0.0, -1.0. 


' 

4 

■ 


'' 


APPENDIX  F  -  CALCULATION  OF  HEAT  LOAD  IN  REACTOR 


Literature  values  of  the  heat  capacities  of  the 
reactants  and  products  are  shown  in  Table  F-1.  The  relevant 
bond  strengths  at  25°C  are  summarized  in  Table  F-2. 

The  heat  of  reaction  for  the  hydrodesulfurization  of 
dihydrobenzothiophene  was  estimated  as  62.4  kcal  per  mole  of 
dihydrobenzothiophene  reacted,  based  on  the  net  change  in 
bond  energy  between  the  reactants  and  products,  assuming 
that  the  C-S  bond  strengths  do  not  differ  greatly  from  the 
value  shown  in  this  table  for  an  aliphatic  chain.  If  the 
heat  of  reaction  for  the  hydrogenation  of  benzothiophene  is 
assumed  to  be  30  kcal/mol,  then  the  net  heat  of  reaction  for 
the  overall  desulfurization  of  benzothiophene  can  be  calcu¬ 
lated  as  about  92.4  kcal/mol. 

Assuming  that  benzothiophene  is  converted  into  only 
ethylbenzene  and  hydrogen  sulfide,  the  heat  capacities  of 
the  reactants  and  products  at  25  and  400°C  were  used  to  cal¬ 
culate  an  estimate  of  the  heat  of  reaction  at  the  higher 
temperature  of  98.4  kcal/mol. 

The  feed  to  the  reactor  under  normal  operating  condi¬ 
tions  consisted  of  hydrogen,  n-hexane,  and  benzothiophene  in 
the  approximate  molar  proportions  of  57:4:1.  Assuming  only 
40%  of  the  benzothiophene  is  desulfurized,  the  average  heat 
capacity  of  the  product  stream  was  calculated  as  about  11.5 
cal/mol-K  at  400°C  and  the  net  temperature  rise  as  53°. 
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Table  F-1  -  Heat  Capacities  of  Reactants  and  Products 


Compound  Heat  Capacity  (cal/mol-K) 


298  °K 

800°K 

67  3  °K 

Hydrogen 

6 . 64  8  3 

7  .  1  6  8  3 

Nitrogen 

6 . 80  8  3 

7  .  1  7  8  3 

Hydrogen  Sulfide 

8 . 28  8  3 

9 . 62  8  3 

n-Hexane 

34 . 20  8  4 

70.36s  4 

61.18* 

Benzothiophene 

31 . 6*  * 

64.0** 

52.8* 

Ethylbenzene 

30 . 70  8  4 

67  .  1  58  4 

57.89* 

*  estimated  by  linear  interpolation 

**  estimated  assuming  the  same  mass  heat  capacity  as 
methyl thiophene  8  5 
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Table  F-2  -  Bond  Strengths  of  Reactants  and  Products 


Bond 


Strength  (kcal/mol) 


H-H 

C-H  (aromatic) 
C-H  (aliphatic) 
S-H  (in  H2S) 

C~S  (aliphatic) 


104.2 

112.0 

98.0 

90.0 

60.0 


All  bond  strengths 8 6are  at  25°C 
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A  change  in  the  liquid  feed  concentration  from  a  mole 
ratio  of  n-hexane  to  benzothiophene  of  4:1  to  4.6:1  results 
in  an  average  heat  capacity  of  12.1  cal/mol-K,  and  a  change 
to  3.2:1  results  in  an  average  heat  capacity  of  10.6 
cal/mol-K . 


. 


APPENDIX  G  -  CALCULATION  OF  KINETIC  CONSTANTS 


The  minimum  apparent  activation  energy  for  the  reaction 
path  through  dihydrobenzothiophene  was  calculated  from  the 
reaction  rates  (cf.  Table  6)  at  350°  and  400°C.  It  was 
assumed  in  this  calculation  that  the  rate  could  be  described 
by  a  zero-order  Arrhenius-type  equation  of  the  form: 

Rate  =  A  exp(-E/RT) 

where 

A  =  preexponential  factor,  mol/g  cat-s 
E  =  apparent  activation  energy,  cal/mol 
R  =  gas  constant,  1.987  cal/mol-°K 
T  =  temperature,  °K 

Thus , 

E  =  -1.987  ln(0. 06243/0. 036107)  =  9123 

673-’  -  623 ' 1 

From  this  value  of  E  and  the  measured  rate  at  400°C,  the 
preexponential  factor  was  calculated: 

A  =  0.06243  exp(+9123/673x1 .987)  =  57.33 


164 


. 


